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FOREWORD 

The pronounced deviations which occur in the aerodynamic properties of 
airfoil profiles below a certain Reynolds number became known to me early, 
and I had therefore planned a systematic study of this "peacetime matter" for 
myself at some quieter time in the future. I am therefore happy that Mr. Schmitz 
has become involved so successfully, even in this time of war, with this task in 
the interest of working groups set up at the schools and of model-airplane flying 
in general. This book not only represents a very remarkable advance in model- 
airplane flying which will undoubtedly be an inspiration in this field to new 
research by scholars, but it also will give very instructive conclusions to the 
fluid dynamics experts through what are in part new and surprising results. 

Goettingen, July 1942. 


Professor L. Prandtl 
Director of the Kaiser Wilhelm Institute 
for Flow Research 
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PREFACE 


The Government Minister for Science and Education made two prizes 
available to the Presidency of the Lilienthal-Gesellschaft fuer Luftfahrtforschung 
( Lilienthal Society for Aeronautical Research) , at their general meeting in 1937 
in Munich, to be awarded every year on the anniversary of Lilienthal's death. 

As explained in detail in the proclamation of 13 April 1938, these prizes are to 
furnish the new generations with incentive for the study of the technology and 
science of aviation. First is the Lilienthal Prize, which gives to qualified 
Gorman secondary-school graduates the material prerequisite for the study of 
aeronautical sciences, and second is the Ludwig Prandtl Prize for Promotion of 
Aeronautical Physics in Connection with Airplane Model Building, awarded for 
the best accomplishments eaoh year in the described field in the secondary 
schools. 

According to the eligibility rules for tho Ludwig Prandtl Prize, the work 
can be done by teachers, by pupils, or as a joint effort. These efforts of pupils 
or joint efforts are directed toward promotion of the selection of future engineers 
and scientists for aviation, since the pupils learn to apply interesting aeronauti- 
cal fundamentals, so that they can be won for aeronautical study with a clearly 
directed vocational goal. 

To supply the schools with the fundamentals constituting a connection 
between aviation physics, model-airplane building, and model-airplane flying, 
the model airplane — as an entity in itself — is subjected in my paper to an 
aerodynamic study. Since I had observed certain contradictions between practi- 
cal data of model flight and die flow laws of large airplanes, I received from the 
National Ministry of Education in reply to my recommendations the mission of 
clarifying the relationships of airfoil measurements at the Reynolds numbers of 
model flight. The measurements were finished in June 1939, and the text was 
finished in June 1940. Publication was delayed by the war. 

The surprising results, which can serve as the first systematic founda- 
tion for scientifically directed model-airplane building, establish the necessity 
for an "aerodynamics of model flight" as a special field of work to complete the 
total physical picture. The first group of my test results to date and their 
explanations are hereby submitted to the German schools and to model-airplane 
building for application. 

Berlin, March 1942. 


F. W. Schmitz 


EXPLANATION OF SYMBOLS 


). Geometric Quantities 


a) Profile 


Symbol 

Unit of 

Measurement 

Explanation 

t 

m 

Chord length, projection of the profile 
outline on the lower-surface tangent 
( "t" is usually the airfoil thickness 
and "c" the chord length in the U. S. ) 

rt 

m 

Maximum profile thickness 

f 

m 

Camber rise of the profile mean lino 
above the theoretical chord of the 
profile mean lino 

r 

m 

Nose circle radius 

d/t, f/t, r/t 

- 

Profile parameters 

X * V y u 

- 

Profile outline coordinates 

Of 

deg 

"Geometric” angle of attack , measured 
between the lower-surface tangent 
and the direction of airflow 

cr 

deg 

Chord angle, angle between the lower- 
surface tangent and the theoretical 
chord of the profile mean line 

a - a + a 
s 

deg 

Angle of attack of the mean line chord 
( theoretical chord) 


deg 

"True" angle of atack of the lower- 
surface tangent to a wing of infinite 
span 

“s-o =a oo + * 

deg 

"True" angle of attack of the mean line 
chord of a wing of infinite span 

a 

deg 

"True" angle of attack, measured from 

0 


the zero-lift direction of airflow 
(c = O') 

p 

deg 

V a ) 

Zero-lift angle of attack of the lower- 
surface tangent 

B = B + <r 
8 

deg 

Zero-lift angle of attack of the chord 

“k 

deg 

Critical angle of attack at breakaway 
of the flow 
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b) Wing 




Symbol 

Unit of 
Measurement 

Explanation 



b 

m 

Wingspan 

F 

m 

Wing area, projection of the wing plan 
on the lower-surface tangent 

A = b/t - b 2 /F 

- 

Aspect ratio 

1/A « t/b = F/b 2 

- 

Reciprocal of aspect ratio 

a 

deg 

Angle of attack of the lower-surface 

IP. 


tangent of the wing model in the wind 
tunnel measurement 

<*5 

deg 

Angle of attack of the lower-surface 
tangent of the rectangular wing with 
A = 5 after taking into account the 
influence of finite airflow diameter 


2. Aerodynamic Quantities 
a) Air 


General subscripts: o on the upper surface of the profile 

u on the lower surface of the profile 
00 in the undisturbed flow 


Symbol 

Unit of 
Measurement 

Explanation 

Y 

kg/m 3 

Air density, weight of unit volume 

g = 9.81 

m/sec 2 

Acceleration of gravity 

P = y/g 

kg * sec 2 /m 4 

Air density, specific mass 

v or v 

00 

m/sec 

Velocity of undisturbed flow 

q = p . v 2 /2 

kg/m 2 

Dynamic pressure 

P 

kg/m 2 

Pressure 

p/q 

- 

Coefficient of pressure 

c = 340 

m/sec 

Velocity of sound 

v/c 

- 

Mach number 

V 

kg • sec/m 2 

Absolute viscosity 

v = rj/p 

m 2 /sec 

Kinematic viscosity 

8 

mm 

Thickness of the boundary layer 
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b) Wing 



Unit of 



Symbol 

Measurement 

Explanation 

Re = v • t/v 


Reynolds number 

Re K = V K ’ ^ V 

- 

Critical Reynolds number 

Re. = v • i/v 

- 

Local Reynolds number at the free edge 

O 


of the boundary layer; i ~ length of 
friction distance 

TP a 

- 

Turbulence factor 

Re Kai/ Re K tunnel 
Re eff = 




Re . TF 

— 

Effective Reynolds number 

meas 




A 

kg 

Lift | v 

— OO 


W 

kg 

Drag 1 1 


M 

M 0.25 

kg • m 
kg • m 

Moment about the front projected point 
on the lower-surface tangent 
Moment about a center of rotation at 

0. 25 t ( quarter chord) behind the 
lording edge and on the mean line 





chord 


M h 

kg • m 

Moment about the center of gravity 

N, T 

kg 

Normal force and tangential force 



components of the resultant air force 

e/t 

- 

Distance of the center of pressure 



from the leading edge 

•&! t 

- 

Distance of the center of gravity from 



the leading edge 

c = A/q ♦ F 

- 

Lift coefficient ("c_ " is generally used 

a 


in the U, S,) 


c = W/q • F 

- 

Drag coefficient ( "c " is generally used 

w 


in the U.S.) 


c = M/ q ♦ F • t; 



referred to the forward 

m 



center of rotation 

C m0. 25 = 


Moment 

referred to the center 

Mo. 25 / 9 • F • t; 


coefficients 

of rotation at t/4 

°mh = M h /q ' P '‘ 



referred to the center of 
gravity 


xiii 



Symbols 


Unit of 
Measurement 


Explanation 


c 

wm 

C wk 

c _ W C ~ C 

w5 wm wkfi 

a m k5 

c 

r 

c , c, 
n t 

C , c 

ns ts 
C f = W/q • 0 

deg 

Drag coefficient for the uncorrected 
measurement 

Correction value for the influence of 
finite stream diameter 
Drag coefficient of a rectangular wing 
with A = 5 , with finite airflow diam- 
eter taken into account 
The angle of attack corresponding to 

C w5 

Coefficient for resultant aerodynamic 
force 

Coefficients of normal force or tangen- 
tial force; center of pressure lying 
on the lower-surface tangent 
Coefficients for normal force and 
tangential force; center of pressure 
lying on the mean line chord 
Coefficient of surface friction (0=2 1) 


c) 

Profile 

c . - c 2 • F/tt ♦ b 2 
wi a 

- 

Coefficient of induced drag 

C a C _ - C , 

w«° w5 wi 

— 

Profile drag coefficient 

a. 23 57.3 • c 
i a 

• F/tt • b 2 

deg 

"Induced" angle of attack 

ot 90 = a s -a i 

deg 

"True" angle of attack of the lower- 
surface tangent, corresponding to 
c 

w®° 

a 8 oo = «co + a 

deg 

True angle of attack of file chord 

dc /da 
a 

- 

Slope of the lift curve 

c a^ c w°° 

- 

Lift-drag ratio ( for infinite span 

c 3 /c 2 
a w 00 

- 

Ceiling factor 




xiv 




d) Propeller 


Symbol 


Unit of 
Measurement 


se 

se 

3f 


m/sec 


Explanation 


Radius 

Blade thickness 

Air velocity in the propeller circle 
Flight velocity 
Peripheral velocity 
Velocity of a blade element on the 
propeller path 

Reynolds number of the blade element 
at radius 0 . 7 r 


e) Model Airplane 


V 

m/sec 

Gliding velocity 

V 

y 

V 
X 

m/sec 

Sinking velocity 

m/sec 

Horizontal velocity 

<p 

deg 

Gliding angle 

*’V T *’ V c a 

- 

Drag-lift ratio 

= W/A = H/B3 


H = height, B = base of the glide 

* tan (p 


triangle 

c 2 /c 8 
w a 

- 

Ceiling factor reciprocal 

G 

kg* g 

Gross weight 

G/F 

kg/m 2 ; 

g/dm 2 

Wing loading 








I. INTRODUCTION 


The measurements of aerodynamic force on airfoil models in the wind 
tunnel by means of the three-component balance are the most important basis of 
aeronautics. The airfoil measurements [3]*) published by the Aerodynamische 
Versuchsanstalt Goettingen ( Goettingen Aerodynamics Research Laboratory) 
should be mentioned as the most widely used. They were made at a Reynolds 
number Re = 420,000. The Reynolds number gives the order of magnitude of 
the flow conditions. The Goettingen measurements correspond to the order of 
magnitude of the smallest manned sailplane. For model flight, where the 
Reynolds number range is below 200,000, these measurements have no more 
validity than for the airfoil flow, for example, of the Messer schmidt record- 
holding airplane, which at v = 210 m/sec and an average chord length of about 
3 m, corresponds to Re = 44 • 10®. For model flight, special airfoil measure- 
ments had to be made in the Reynolds number range mentioned. 

The range of Reynolds numbers for model aviation covered by the follow- 
ing airfoil measurements extends from Re = 21,000 to 168,000. Five profiles 
were measured, chosen so that their results showed the greatest possible con- 
trast, to delimit the problem of profile properties in this range of Reynolds 
numbers. 

In the range of Reynolds numbers below 150, 000, each profile passes 

through a critical range in which the boundary layer becomes turbulent, and 

the previously laminar separated flow attaches itself so that the lift coefficient 

c becomes greater, the drag coefficient c smaller, and thus the drag-lift 
a w 

ratio c /c formed from the two becomes suddenly larger, so that the model 

w a 

airplane first becomes able to fly. 

From the comparison of test results for the five profiles there is found 
for the model airplane: 

1. Rules for profile choice in regard to the parameters: profile thick- 
ness d/t, mean camber f/t, and leading-edge radius r/t. 

2. Rules for the choice of wing plan. 

3. Rules for the choice of the model airplane propeller with respect 
to the blade profile and blade plan. 


*) The numbers in brackets refer to the bibliography at the end of the 


book. 


4. Comparison of the measurements on profiles N 60 and N 60R shows 
to what extend the S-shapcd curvature of the mean line affects the 
moment coefficient or aids in achieving constancy of the center of 
pressure at low Reynolds numbers. 

5. The results turn model -airplane building again to the study of the 
biological example, bird flight, in the sense of Lilienthal, because 
bird and model airplanes fly in the same range of Reynolds numbers. 

6. Measures for artificial turbulence by wall roughness, turbulence- 
producing edges, or wires stretched in front of the wing on the 
model airplane, in analogy to the classic sphere research with and 
without wire rings (by Prandtl and Wieselsberger) prove to have a 
favorable effect, in contrast to their harmful influence on the large 
plane. 

7. An attempt is made, using the study of the individual phenomena of 
flow around an airfoil, to explain the flow causes in the critical 
region. 

The measurements first became possible after the airflow could be made 
largely free of turbulence by modification of the wind tunnel used. 

Because these measurements were to be used on one hand for physics 
instruction and on the other hand for model airplane building, that is, to direct 
the interest of youth, through the teachers, to aviation research, it appeared 
advantageous to formulate the basic laws of boundary layer flow and flow around 
an airfoil necessary for explanation of the test results. 

The basic concepts of aeronautical physics required in the syllabus of 
physics of the secondary schools were presupposed for understanding of this. 

An example is the book Einfuehrung in die Physik des Fliegens ( Introduction 
to the Physics of Flight) by Schuett, Verlag C. J. E. Volckmann, Nachf. E. 
Wette, Berlin-Charlottenburg ) . In addition, knowledge of the quantitative 
application of these basic ideas - corresponding, say to the Flugphsikalischen 
Arbeitsbuch ( Workbook in Aeronautical Physics) by Kisse, Verlag G. Freytag, 
Berlin, and B. G. Teubner, Leipzig - is presupposed. The author published 
a short discussion of the basic problem in the No. 3 issue of the year 1942 of 
Unterrichtsblaetter fuer Mathematik und Naturwissenschaften ( Instruction 
Sheets for Mathematics and Natural Sciences) , Verlag Otto Salle, Frankfurt 
a.M. To allow the teacher to explain the interesting relationships to model 
fliers less well-educated mathematically, the January 1942 issue of the 
magazine Luftfahrt und Schule (Aviation and School) , Verlag C. J.E. Volckmann 
Nachf. E. Wette, Berlin-Charlottenburg ) , began a series by the author under 


the title "Fluid Dynamics Study of the Model Airplane," with simple examples 
of calculation, and flight and school experiments. 


II. FUNDAMENTALS AND FORMULATION OF THE PROBLEM 

1. Comparison: Sailplane, Model Sailplane, Bird 

a) The Similarity Law ( Re » constant) 

The bird is the natural prototype of the model airplane. The contempo- 
rary model airplane builder, however, in contrast to Lilienthal, thinks little 
about bird flight. His prototype is the large airplane, which he seeks to copy 
as closely as possible in the model. This tendency in development led to 
contradictions between the facts of experience of model flight and the aerody- 
namics of the large airplane. Although there is an aerodynamic relationship 
between bird flight and model flight, since both occur in the ranges of similar 
Reynolds numbers, fulfillment of the similarity law of Reynolds is not possible 
between the small free-flying model airplane and the large airplane, as a com- 
parison will show. 

A model airplane, which represents the geometrically similar imitation 
of a large sailplane, will give a poorer drag-lift ratio in comparable flight 
(unaccelerated gliding flight in still air) the smaller it is. If the sailplane has, 
for example, a drag-lift ratio of: 

( C w /C a) F = 1/20 * 

then experience shows that with a model airplane on a 1:10 scale the gliding dis- 
tance is about half, and thus the drag-lift ratio expected would be: 

(°v/ C a) M = 1/l0 ‘ 

( The subscript F = airplane and M = model. ) 

The deterioration in drag-lift ratio is caused, as is well known, by the lack of 
similarity of the streamlines, because the inertia forces of the air p v 2 /t, 
effective along the wing chord t, decrease more than do the viscosity forces of 
the air, r) v/t 2 , when the model is reduced. Only if it is possible to hold constant 
the ratio of the two, that is, the Reynolds number: 

inertia force P v 2 / 1 P vt yt 

e ~ viscosity force tj v/t 2 1 ? v 
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(p = air density, v = velocity, t = wing chord, ri => absolute viscosity, 
r/p = v => kinematic viscosity) 


will this similarity law of flow be fulfilled, allowing the model airplane to achieve 
the same flight performance as the large airplane. Therefore if the model air- 
plane had a wing chord t of Vio the wing chord of the large airplane, it would have 
to glide at 10 times the speed. Ten times the speed in unaccelerated gliding 
flight is not possible for free flying model airplanes, however. The reason for 
this is as follows. From: 


gliding speed v » 



there follows for v^ = 10 v^, the equation: 



and from this the wing loading of the model airplane: 

(G/F) m = 100(G/F) p 


Because a 1:10 linear reduction means an area reduction of = F^/lOO, the 

model airplane weight G required for 100 times the wing loading is calculated 
from 



100 °F . °M 100 °F 

— - f^ckT -iTT 


The fulfillment of this similarity law thus leads to the practically impossible 
requirement that the model airplane and the airplane must have the same weight. 


* ) The root under the fraction stroke gives the coefficient of the resultant 

air force c through use of the Pythagorean theorem for the triangle short sides 
r 

c & and c^. For flat gliding flight, c^ can be taken as approximately c & . 

4 



To mako model flight possible, model airplane building is compelled to 
go in the opposite direction, using lightweight construction to make the wing 
loading and thus the sinking speed as low as possible, to what extent the "gust 
sensitivity" of model flight permits. Fluid dynamics therefore could learn 
little in quantitative fundamentals from measurements on free-flight small model 
airplanes (see [5 J , vol. IV, part 2, A. Betz, p. 211) because it is much more 
suitable to hang the model airplane in the wind tunnel, allow the air to flow by 
at say ten times the speed, and to measure the aerodynamic forces acting on the 
model airplane by means of the wire suspension attaching it to the balances. 

The essential experimental fundamentals for aerodynamic calculation of the 
airplane were found through this model testing in the wind tunnel. There was 
no practical economic interest in measurements at the low speeds of free-flying * 
model airplanes, so that model airplane building - aerodynamically considered - 
was left to the hazard of testing. If, however, as described in this work, there 
are measurements both at the Reynolds number of the large airplane as well as 
at the Reynolds number of the model airplane, then it is possible to calculate the 
flight performance and flight characteristics of the model airplane and to test 
the correctness by flight measurements. 

« 

In the model airplanes common today, about five major classes can be 
differentiated, and the Reynolds number ranges of these are given in List 1. To 
show the aerodynamic relationship, three types of soaring birds in the same 
range of Reynolds numbers are added, and in addition, as a contrast, the 
Reynolds numbers of various airplanes. 

Note on the calculation of Reynolds number: 


The kinematic viscosity ( also called viscosity modulus) v has the 
dimension L 2 / T . 

For air at 13° and 760 mm Hg, v = 0.143 cm 2 /see = 0.0000143 m 2 /sec. 
For water at 20° and 760 mm Hg, v = 0. 01 cm 2 / sec = 0. 000001 m 2 / sec. 

The kinematic viscosity of air is about 14 times as great as that of 
water, so that in water V14 of the speed is required to achieve the same 
Reynolds number as in air. 

The numerical values for air density p and for v as a function of tempera- 
ture and air pressure can be taken from data in the literature ( (3) I, pp. 136 
and 136, or from [5] vol. IV, part 2, p. 113) . The Reynolds number change 
caused by atmospheric fluctuations at a given reading of the dynamic indicator 
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ia, for example, about +10 percent from B = 740 mm Hg and 20° C to 780 mm 
and 0°C. 

For the upper range of measurement of the following airfoil measure- 
ments there resulted at a maximum velocity of v = 27 m/ sec with a chord length 
t = 90 mm: 


_ 2700 cm/s • 9 cm ^ 27 m/s ■ 0. 09m 
e max ~ 0, 143 cm J /s ” 0.0000143 m 2 /s 


170,000 


List 1 

Ranges of Validity of Reynolds Numbers*) 



Mean Wing 
Chord t ( mm) 

Velocity 
v ( m/s) 

Re = v • t • 70 

Zanonia plate 

40 

1 

2800 

Gliding butterfly 

50 

2 

7000 

Indoor model airplane 

60 

2 

8400 


150 

4 

42,000 

Swift 

30 

min 6. 1 

12,800 



max 39* *) 

82,000 

Small model airplane 

100 

3 

21,000 


150 

6 

63,000 

Medium size model airplane 

150 

4 

42,000 


200 

10 

140, 000 

Herring gull 

140 

10 

100, 000 

Large model airplane 

200 

4 

56,000 


300 

10 

210,000 

Albatross 

200 

16 

224,000 

High-performance sailplane 

900 

v_ = 10 

Lt 

630, 000 



v =60 

3,800,000 



max 



* ) The Reynolds numbers are calculated from the mean wing chord. 
For a tapered wing, the Reynolds number at the wing tip is correspondingly 
smaller, and larger at the wing root. In this, = landing velocity, 
v = flight velocity. 

**) According to R. Schmidt ( see [18] , p. 84) , the swift does not soar, 
but like the swallow, makes only short gliding flights. 


* 
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If the comparison between model airplane and airplane is made in air of 
the same kinematic viscosity v , then the similarity law will be satisfied if the 
characteristic E = v • t is kept constant. This simplified but incomplete form 
of the similarity law has not established itself and was not used here. If v is in 
units of m/sec and t in mm, then Re is roughly 70 times the value of E; for 
example: 


Re = 27 m/s • 90 mm . 70 = 170,000 

It should be mentioned in this connection that changes in flow dependent 
upon Reynolds number are called "Reynolds number influence”; rounded bodies 
are "Reynolds number-sensitive"; sharp-edged bodies are "Reynolds 
number-insensitive. " 


b) Wing Loading and Drag- Lift Ratio 

For model airplane building there is technically no difficulty in reducing 
the wing loading by lightweight construction to 

( G/ F) = 4-0.1 kg/m 2 = 40-1 g/dm 2 , 

M 
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the lower value corresponding to the Indoor model airplane covered with thin 
film. This is about V 6 to V 2 oo of the wing loading of the sailplane, which is 
between 20 and 10 kg/m 2 . The wing loading of birds varies between the 16 kg/m 2 
of the albatross and the 1.7 kg/m 2 of the swallow; that of butterflies is between 
0.17 and 0.1. Here model airplane building follows the prototype of nature. 

The low wing loading considerably reduces the sinking velocity: 



and in spite of a low drag-lift ratio, an equal or longer flight time is achieved 

compared with that of the sailplane. The necessary performance required for 

horizontal floating in a thermal or by means of an engine is found, as is well 

known, at the minimum of v or / o 8 /c 2 \ in the polar plot. Flight at this 

y y a w /max r 1 

angle of attack therefore corresponds to the longest gliding time, while the some- 
what smaller angle of attack corresponding to c ft / c w max gives the longest glid- 
ing range. These two values are shown for two examples in Figure 1. For the 
same wing loadings ( G/ 1) M = ( G/ F) F » the sinking velocities are related as 

the roots of the ceiling factor reciprocals c 2 / c 3 . Applied to the example 
given, there result: w h 



The model airplane in this case will have a sinking velocity n/ 1T = 2.8 
times as great. To bring the model airplane to the same sinking velocity, its 
wing loading must be reduced to Vg; then it achieves the same gliding duration 
and the same climbing ability in a thermal as the large airplane. The ratio of 
gliding velocities is then: 


V M _ V + V M „ 1 

v f n 1 2 

■W + V F 

and the ratio of the Reynolds numbers is: 
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Figure 1. Comparison of two polar plots made at Re -- 42,000. 
This Reynolds number corresponds to the order of magnitude of, 
for example, an indoor airplane model, which with a chord length 
t - 150 mm glides with a velocity of v = 4 m/sec. Comparison 
shows the characteristic superiority of the curved plate at low 
Reynolds numbers ( supercritical flow conditions) compared with 
an airfoil profile of 12 percent thickness ( subcritical separated 
flow) . The polar line gives at the point of t&ngency the necessary 
angle of attack for the best drag-lift ratio (that is, for the 
greatest gliding range) . On the other hand, achievement of the 
maximum gliding flight time, that is, the lowest sinking velocity, 
always requires a somewhat larger angle of attack, which corre- 
sponds to the maximum of the ceiling factor c 3 /c 2 . 

& w 


(Re) 

(Ri) 


M 

F 


(v • t) 


M 


(v ♦ t) 


_ 1 _ 
20 * 


The simdng and gliding velocities are given in Figure 2 in dimensionless numbers 
as a function of the drag-lift coefficient for various wing loadings. The relation- 
ships can be easily seen in this diagram for transfer of tests on the airplane to 
the model airplane and vice versa, because the points of constant wing loading 
plotted in a log -log form, are straight lines. The comparison is applied to 
practical examples in List 2. In addition, the numerical values for some soar- 
ing bird types are given there. The aerodynamical superiority of the bird 


9 









indicates that essential improvement possibilities exist for the model airplane 
if it is made to resemble the bird wing to achieve better drag-lift ratio. The 
mechanical imitation of the bird's wing will achieve only one flight condition, 
however, while the bird in soaring flight can change the wing profile, the wing 
plan, the wing twist, the dihedral and sweepback of the wing, the size of the 
lifting surface, the position of the center of gravity, the position of the center of 
pressure, and the tail feather setting, all for suitable utilization of any possible 
flight situation. 


Figure 2. Chart for change of wing loading G/F, drag-lift 

ratio c /c , sinking velocity v or horizontal velocity v . 
w 8» y 

The left lower corner corresponds approximately to the butter- 
fly, the upper right to the high-speed or very large airplane. 
Example: If the drag-lift ratio of a model airplane ^c^c^ 2 

is poor in comparison with a large airplane ^c c a ^i> f° r 

example, if 2 : ( c w /° a )i = Vio* V20 * 2 ( cf » horizontal 

scale at the upper left) , then proceeding from the center of 
the chart outwardly toward the upper left upon intersection of 
the line v Jv with the perpendicular under the 2 there is 

y 3 

an increase in sinking velocity v = 2. 8 v . , and at the 

yz yi 

intersection point underneath with the line v ^/v there is an 


increase in horizontal velocity v^ = 1.4 v^. If however, die 

model airplane has the same sinking velocity as the large plane 
(a point perpendicularly under, on the 1.0 horizontal) then the 
wing loading of the model airplane (G/F) 2 = V 2 3 * Vs of the 
wing loading of the large plane (as seen by going along the 
new v line from the point mentioned on the 1. 0 horizontal to 

y 

the lower right at the intersection with the scale of the wing 
loading) . The horizontal velocity for this is then v^ =0.5 v^ , 
that is, half as great as in the large airplane. 


The scales correspond exactly to the 12-cm slide rule, so that 
Intermediate values can be read off on the slide-rule slide 
because the line grid gives directly only the changes in wing 
loading which correspond to the geometric series of the num- 
bers 2 and 10 (dashed lines) . 
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++) The fuselage width is not assumed to be load-carrying for the mean chord length. 



List 2 ( Concluded) 















2. Test Conditions and Test Installation 


The measurements were made in the wind tunnel of the State 
Engineering School in Cologne. *) 

The following difficulties are found for airfoil measurement in the wind 
tunnel in the Reynolds number range of model flight: 

1 . The aerodynamic forces are so small at low angles of attack that it 

is necessary to specify a measuring accuracy of tenths of a gram 

for the balances. For the balance used, the lower measuring limit 

was therefore Re . » 20,000, 

min 

2. At the chosen chord length t - 90 mm, a dynamic pressure reading 

q = 0.75 mm water column corresponds to this Re . . The "Debro" 

min 

miniscope of the de Bruyn firm was used; it permitted a reading 
accuracy of 0. 01 mm water column. 

3 . The key point in the work is in the study of the effect of the charge 
of laminar boundary layer flow on the upper side of the airfoil into 
a turbulent flow. Because the flow state of the free atmosphere can 
be considered to be practically turbulence-free, these effects in the 
wind tunnel appear in undistorted form only when the airstream is 
laminar. This condition is difficult to fulfill because the turbulence 
of a wind tunnel is probably always greater than that of free atmo- 
sphere. The wind tunnel used for the tests in Cologne appeared to 
be too turbulent in comparison with a Goettingen measurement in 
the first check measurement. As a rough check on turbulence, the 
flame probe showed a long eddying flame over the whole stream. 

The measurement [9J of the critical Reynolds number of the sphere 
Re^ = v ♦ d/v, commonly used as an indication of wind tunnel 

turbulence, gave for a polished wooden ball of 18 cm diameter: 

p/q = 0 at Re, = 2. 07 • 10 s . 

00 k 

Pressure p was measured at the rear support tube of the sphere, and 
q^ = the dynamic pressure in the undisturbed flow. 


* ) The wind tunnel was made available by Professor Grunewald, the 
principal, and Engr. D. Eck for the test and for the necessary modification. 
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After the basic changes in the wind tunnel, taking into account the guide- 
lines developed by Prandtl on "Setting up Satisfactory Airstreams" (see [5], 
vol. IV, part 2, p. 65) and following the proposal of Max Kramer, Eng. D. , 

Adlershof, p/q = 0 was established at Re =■ 3.8 • 10 B , the dynamic pressure 
00 K 

distribution in the horizontal diameter D of the stream, at distance D/2 ahead of 
the nozzle was: 


A q/q^ " 0*03 (0.10) , 

and the static pressure in the stream center at distance D/2 was: 

p/qm « 0.004 ( 0.015) . 

The earlier values are placed in parentheses. 

Flight tests of the DVL [ Deutsche Versuchsanstalt fuer Luftfahrt 
( German Institute for Aeronautics Research) 1 made by S. Hoerner ( 9 ) with 
spheres gave Re in still air as 4.05 • 10 5 . This number as the zero point of 

wind tunnel turbulence gave a turbulence factor [10] of the wind tunnel of: 


Re, in still air 
k 

Re, in the wind tunnel 
k 


= 1.06 (1.99) . 


In connection with turbulence factor, an effective Reynolds number is 
defined in NACA Report No. 558 as 


Re eff ^measured 


TF , 


in the assumption that a process observed in the turbulent wind tunnel would be 
seen in a turbulence-free flow at a Reynolds number depending on the turbulence 
factor. The NACA profiles N 60 and N 60R used in the following measurements 
were measured at Re « 3, 000, 000' in the high-pressure wind tunnel (VDT) 
known to be strongly turbulent. At a sphere characteristic of Re^ = 1.55 • 10 s 

(see Figure 11) the results for the NACA measurement were TF = 2.62 and 

Re „ = 3,000,000 * 2.62 « 8,000,000. This Reynolds number correction is 
eff 

not used in the following measurements, because a change of Reynolds number 
by 6 percent would not essentially alter the overall picture. 



3. Subcritical and Supercritical Flight State for the Model Airplane 

In contrast to flow around the airfoil on large airplanes, which is 
always supercritical, ’'subcritical" flow, in addition to supercritical, can 
appear in flow around the model airplane airfoil in the Reynolds number range 
under 100,000; at the subcritical condition where lift A is small, drag W is 
disproportionately large, so that the drag-lift ratio is very poor and is the 
cause of the failure of many model flight experiments. The following measure- 
ments now show why a model airplane can fly only in the supercritical flight 
state and how through a suitable choice of the airfoil profile the supercritical 
flight state can always be realized. Figure 3 gives the results of a measure- 
ment on an airfoil for an angle of attack of a. ~ 6° as an example. Instead of 
A and W, the dimensionless coefficients, lift coefficient c , drag coefficient 

c w » and the lift-drag ratio c a /° w » art plotted as a function of Reynolds number 

Re. The practical application of the measurement will be shown for an example. 
Below the critical Reynolds number Re^ the lift-drag ratio is 4 to 5, and in the 

supercritical region it is 10 to 11, that is, for 1-m drop the airplane model 
glides only 4 to 5 m subcritieally but more than twice as far supercritically. 

In the use of this airfoil orofile the model airplane with an assumed velocity of, 
for example, v = 6 m/sec would have to have a chord of at least 200 mm. The 
model airplane would then fly (in still air) at Re = 6 ♦ 200 ♦ 70 = 84,000 and 
thus in the neighborhood of the critical Reynolds number. Velocity loss from 
6 to 5 m/sec would suffice to bring the model airplane into the subcritical flight 
stage at Re = 5 ♦ 200 ♦ 70 => 70, 000; it would pancake or go into a spin. 

These phenomena in flow around the airfoil are controlled by the flow 
state of the thin friction layer near the wall, called the "boundary layer” by 
Prandtl, who first explained its effect in 1914 [1], Because of the viscosity of 
the air, the air molecules adhere to the body wall, act as a brake on the air 
particles moving above diem, anu thus generate a friction resistance. The 
friction loss in the boundary layer, the boundary layer thickness, the velocity 
distribution in it, and above all the ability to draw energy from the external 
flow around the foil can be basically different, depending on whether the bound- 
ary layer is laminar ( from lamina = layer) or turbulent ( eddying) . In the sub- 
critical flow state, the boundary layer flow is laminar on the airfoil upper side 
and in the supercritical flow state it is turbulent. Both flow forms have inter- 
changeably a favorable and an unfavorable characteristic. 





a) The Laminar Boundary Layer 

The laminar boundary layer ( Figure 4) , which is peculiar to the range of 
low Reynolds numbers, is very thin. It causes little frictional drag; wall rough- 
ness which is not excessive is smoothly covered by it, so that the nature of the 
surface on a small airplane model plays a different role from that on the large 
airplane. The laminar layer does not have the ability, however, to take energy 
from the outer flow, so that in spite of low friction its kinetic energy is soon 
used up. Consequently, it adheres poorly to curved surfaces in a divergent flow, 
such as at the rearward curvature of the airfoil upper side ( suction side) , 
because there, deceleration and pressure increase always prevail, through 
which it cannot penetrate. From the region of higher pressure, and therefore 
from the airfoil trailing edge, there follows shortly after beginning of flow a 
retrograde movement in the boundary layer which, like a wedge, moves under 
the decelerated boundary layer material, separates the flow from the wall, and 
through eddy formation creates higher drag behind the wing ( Figure 5) . The 
point of this reverse flow wedge is called the "separation point, " which at 
"laminar separation" lies at the point of the beginning pressure increase and 
therefore approximately at the highest point of curvature, moving forward with 
increasing angle of attack. If laminar separation exists in the whole range of 
angles of attack (which for the given airfoil profiles can be only for Reynolds 
numbers less than 100,000, and therefore only on model airplane airfoils) this 
flow state is called "subcritical. " The pertinence of the origin of laminar flow 
at low Reynolds numbers is understandable through the physical significance of 
the Reynolds number, which gives the ratio of inertia force to viscosity force; 
that is, at low Reynolds numbers the viscosity forces of the boundary layer 
conditions prevail and at higher Reynolds numbers the inertia forces. 


Figure 3. Example of an aerodynamic force measurement 
on an airfoil model at a given angle of attack and increasing 
Reynolds number. The transition from subcritical (laminar 
and separated) to supercritical ( turbulent and attached) flow 
state of the upper surface boundary layer occurs at the criti- 
cal Reynolds number, which is Re^ = 63,000 here. The lift 

coefficient c increases suddenly there, the drag coefficient 

a v 

c decreases and the lift-drag ratio c /c becomes about 
w aw 

three times as favorable (from 4 to 12) • 


TRANSITION POINT 



Figure 4. Variation in the boundary layer thickness 6, along 
a flat plate in laminar approach flow. The figure at the lower 
left is the velocity profile of the laminar boundary layer; at 
the lower right is the velocity profile of the turbulent boundary 
layer. 


b) The Turbulent Boundary Layer 

As the Reynolds number increases (that is, with increase in speed v, 
when the chord t and the kinematic viscosity of the air are constant) the bound- 
ary layer after initial laminar flow begins to be turbulent at a critical Reynolds 
number Re^, and the lift suddenly increases, the drag decreases, and thus the 

drag -lift ratio becomes better by several times, and then improves slightly 
further as the Reynolds increases. This state of affairs, which explains the 
cause of the sudden changes shown in Figure 3, gives as a first basic conclusion: 
the model airplane will give optimum performance in its size class only if its 
flight state is supercritical. The following measurements answer the questions: 
"On what is the critical transition dependent, and how can it be artificially 
influenced?" It may be said beforehand that the most important result is that the 
attainment of supercritical state is always possible by 

a) a suitable choice of airfoil profile, 

b) by measures for aritificial turbulence on any arbitrary profile. 
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SEPARATION POINT 


Figure 5, Velocity profiles of the boundary layer on the upper 
side of the airfoil before and after the separation point. The 
dashed line limits the reverse flow wedge. 


The superiority of the supercritical flight state is the result of the 
peculiar ability of the turbulent boundary layor to take energy from the outer flow 
and transport it to the wall as a result of mixing motions of the small vortices. 
Through impulse exchange on a "mixing length" ( see ( 2 ] , p. 94) perpendieular 
to the principal direction of flow, the boundary layer is able largely to overcome 
the pressure increase at the rearward curvature of the airfoil upper side, so 
that the point of separation is brought closer to the trailing edge, the higher is 
the energy content of the outer flow or, expressed in a form easier to understand, 
the higher becomes the "drag effect of the outer flow" as Reynolds number 
increases. 

The previously laminar separated boundary layer attaches itself to the 
upper side because of turbulence , the vortex street behind the airfoil becomes 
narrower, and thereby the pressure (eddy) drag (also called form drag) becomes 
suddenly smaller. In spite of this, in the total drag, which as is known consists 
of 1. the pressure drag, 2. friction drag at the body surface, turbulent friction 
on smooth walls can be three times as large as laminar friction, and on rough 
walls even more, as has been found from measurement of friction drag c^ on 

thin flat plates ( Figure 6). The "vorticity" of the friction layer causes, in 
addition to increased friction loss, a sudden increase in boundary layer thickness 
5 at the transition point ( Figure 4) . Both are characteristics of an apparently 
increased viscosity. The thickness 6 of the laminar boundary layer increases 
with the V 2 power of friction length i, and the thickness of the turbulent boundary 
layer increases with the % power. This considerable thickness increase 
because of decelerated boundary layer material finally leads, with declining 
drag effect of the outer flow, to "turbulent separation" at the curved upper side 
of the airfoil, again with a reverse flow wedge detaching the flow from the wall 
before the trailing edge ( Figure 5) . The point of separation migrates only a 



slight amount at first with increasing angle of attack and then suddenly forward 
at a critical angle of attack; "breakaway of flow. " 



Figure 6. Drag coefficient c^ of surface friction on flat plates 

iti the laminar and turbulent friction layer as a function of 
Reynolds number. 

4. Dependence of the Critical Reynolds Number 

A transition in flow from the subcritical to supercritical state, which 
the critical Reynolds number Re^ characterizes, has not been studied very much 

to date for the airfoil because Re fe , lying in the region of model flight, is of 

little importance to the large airplane. The mathematical treatment of this 
region is confronted by much greater difficulties because farther the consider- 
able separation phenomena remove the flow from the easily investigated 
quantitatively theoretical potential flow, the poorer is the drag -lift ratio of the 
airfoil. In contrast to this, the subcritical and supercritical phenomena in flow 
in tubes along flat plates and spheres has been largely explained and calculated; 
for this reason, a brief study will be made of the conclusions resulting from this 
for the model-airplane airfoil. 


a) The Critical Re Number for the Flat Plate 

If a plane, infinitely thin, smooth plate "friction plate" is in i. laminar 
flow parallel to the wall, so that no forces perpendicular to the wall but only 
frictional drag occurs, then a boundary layer with a laminar approach flow 
forms; its transition takes place after a friction length ( Figure 6) a. 
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Re = (v • i/v « 500,000*)) (see [3J, III, p. 5). This value of Re, on the 
k k 

friction plate cannot be related to the model airplane, however, as can be seen 
immediately: the model airplane used in the following measurements has a 
chord length of t » 90 mm. To achieve a transition in a 90-mm friction length, 
the velocity v would have to be 500,000/(90 • 70) « 80 m/sec. The smooth 
flat steel plate used in the measurement here, 2. 6 mm thick and 0. 5 mm nose 
radius at the leading edge, exhibits supercritical flow at the lower limit of the 
measured range, that is, even at v -3.2 m/sec (Re = 20,000); this is there- 
fore produced by influences other than surface friction. If the plate is thick and 
the leading edge rounded ( Figure 7) , then through local acceleration in combina- 
tion with a local separation (transition vortex) , this produces an earlier transi- 
tion of the flow, as can be seen from the suddenly increasing boundary layer 
thickness in Figure 8. Wieselsberger (see [3] , I, p. 123) found in his first 
measurement of friction drag on plates 
with rounded leading edge in the low- 
turbulence Goettingen wind tunnel only 
supercritical friction drag (corre- 
sponding to the upper line in Figure 6) . 

To permit measurement of laminar 
friction, according to Hansen [ 11] the 
plate must be tapered to a sharp edge 
(Figure 8). This situation immedi- 
ately reverses itself, however, if the 
sharpened plate is set at an angle as 
an airfoil so that flow takes place at a 
slight angle to the sharp leading edge. 

The stagnation point then moves away 
downstream from the leading edge by 
a small amount along the lower surface, 
with the consequence that the lower sur- 
face flow becomes turbulent because the 
sharp leading edge now has the effect of 
a ’’turbulent edge." This makes it 
understandable why the aerodynamic 
force measurements on plates and thin 
profiles with sharp leading edge give 
supercritical values even at the small- 
est angle of attack. The turbulence 
edge on the flat plate of course has 

*) According to more recent research on thin plates, this number dis 
places toward about 1 million, depending on the shape of the leading edge in 
laminar approach flow, as shown by the dashed line in Figure 6. 



Figure 7. Pressure distribution 
at the rounded leading edge of a 
flat plate with flow parallel to the 
wall. 
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Figure 8. Development of boundary layer thickness 6 along 
plates of various thickness and with various shapes of the 
leading edge ( according to Hansen) . 


this favorable effect only to about a «* 7 to 10° , because then the sharp edge 
becomes a "knife edge, " as described elsewhere. With the plate standing per- 
pendicular to the flow, and thus at a = 90° , at any Reynolds number and at any 
state of turbulence of the airflow, the same drag coefficient is measured 
( Figure 9) . Sharp-edged bodies have no critical Reynolds number and they are 
"insensitive to Reynolds number" because the point of breakaway is fixed, so 
that here the model law of Reynolds has either no or limited validity. 



b) The Critical Re Number for Rounded Bodies 

In contrast to this, all rounded bodies - sphere, ellipsoid, streamlined 
body and round-nosed airfoil - all show a c jump at a Re. and have "Reynolds 

number and turbulence sensitivity, " as Figure 9 shows. According to Figure 9, 
Re decreases with increasing fineness. Applied to airfoils the following con- 

K 

elusion results: the finer the profile, the lower the Reynolds number at which 
the transition from subcritical to supercritical flow occurs. Because basic 
conclusions, both for the explanation of problems with the model airplane as 
well as for the utility of the wind tunnel employed for the existing measurement, 
result from known phenomena on the sphere, it will be valuable to summarize 
these phenomena briefly here. The flowing particles approaching the stagnation 
point Si of the sphere’s leading side ( Figure 21) come to rest, so that their 
kinetic energy is converted into pressure. At a hole through the sphere center 
and from the stagnation point this pressure can be measured by a manometer 
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Figure 9. Reynolds -number sensitivity and turbulence sensitivity 

( according to Prandtl-Tietjans) . The circular plate ( upper part 
of the figure) is not sensitive to Reynolds number: = constant. 

The rounded bodies are 1) sensitive to Reynolds number: indicated 

by c jump; 2) turbulence sensitive: with a turbulence grille 
w 

( dashed lines) , the c jump occurs at a lower Reynolds number. 

w 

For the sphere, that Reynolds number is critical at which 

c = 0.3: here it is at 2. 58 • 10®. The turbulence factor of the 
w 

wind tunnel is: 


Re, 


TF = 


k air 


4.05 


Re 


k tunnel 


2.58 • 10 


10 6 

5 » 1.6 




connected to the rear side of the sphere; this represents the "total pressure" 

P, similar to a Pitot tube measurement. After subtraction of the "static 
pressure" p^ of the undisturbed flow from it, there results the "dynamic 
pressure" q^, in accordance with the Bernoulli equation: 

P=p+q or q=P-p . 

In the free -flow wind tunnel, the static pressure is negligibly low. In the 
Cologne wind tunnel which was used, it was only 0.4 percent of the dynamic 
pressure ( p =* 0. 004 q^) after the conversion of the installation into undis- 
turbed flow, as has already been described In Subsection II 2, so that there can 
be written: q^ » P. In flow past bodies, on the other hand, the static pres- 
sure, as the wall pressure, plays a considerable role. In Figure 21 the space 
around the sphere is divided into flow tubes whose cross section is reduced 
most strongly at the sphere equator by constriction of the flow between the 
sphere and the outer flow. It is assumed that no fluid passes through the 
imaginary walls of the flow tubes, and therefore in a unit of time the same 
quantity of flow must pass through each cross section ( continuity requirement) : 
for a cross-section reduction to half, therefore, the velocity must be twice as 
large. Since the dynamic pressure increases with the square of velocity, how- 
ever, the dynamic pressure at this place in the flow must be four times as high: 

dynamic pressure q = — v 2 and q = •— , 

because for air at zero-meter height the air density p * Vg. According to the 
law of Bernoulli, however, the total pressure P in the flow tube remains 
constant: 


P = p + q - constant. 

Therefore, with increasing dynamic pressure the wall pressure must corre- 
spondingly decrease, and therefore at the sphere equator there is a negative 
pressure. These phenomena at the leading side of the sphere always take the 
same course no matter whether an ideal flow of a frictionless liquid is involved, 
corresponding to Figure 21 , or whether flow takes place around the sphere sub- 
critically or supercritically ( Figure 10) . The streamline pictures of the three 
flow forms named are first differentiated by the phenomenon of flow on the rear 
side of the sphere. 
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SUPERCRITICAL 


Figure 10. a) Subcritical flow around the sphere: the boundary 
layer flows in laminar form; separation occurs somewhat ahead 
of the sphere equator; negative pressure prevails at the rear 
side of the sphere and the large turbulent area produces high 
resistance: c^ = 0.48. b) Supercritical flow around the sphere: 

in the turbulence-free airstream the boundary layer on the smooth 
sphere is turbulent at Re = 405,000; the flow remains attached, 

K 

and the turbulent region is considerably smaller. Positive pres- 
sure prevails at the rear side; = 0. 08. 


At the rear side of the sphere in a frictionless fluid, that is, a fluid 
without viscosity, the streamline picture would be a mirror image of that at the 
front side of the sphere, in accordance with Figure 21. The conversion of 
pressure into velocity energy, effective up to the equator from the front stagna- 
tion point Si, would be reversed without loss at the rear side, so that a second 
stagnation point S 2 would occur there. This lossless, practically impossible, 
flow ( called potential flow) leads to the suprising finding that the resistance = 0 
because there are no friction or vortex losses: the front-side pressure and the 
rear -side pressure would balance one another. *) 


The streamline picture of subcritical flow around a sphere shows the 
greatest deviation from this ideal flow ( Figure 10a) in the large turbulent 
region at the rear side of the sphere, in which a reverse flow wedge separates 
the boundary layer of the sphere front side even upstream from the equator. 


* ) This streamline picture of potential flow can be exactly calculated 
and plotted by means of the stream function. In addition, this picture can be 
produced by "creeping" laminar flow around a thin circular sheet between glass 
plates, for example, in the Pohl streamline apparatus, even though the physical 
causes are different. 
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This subcritical state is caused by the oontinued laminar nature of the boundary' 
layer on the sphere front side; it transports no energy to the layer decelerated 
by wall friction and therefore cannot overcome the pressure increase beginning 
at the equator but instead separates from the sphere surface. At the rear side 
of the sphere, a negative pressure exists in the turbulent region (as indicated 
by the minus sign in Figure 10a) and can be measured by a second manometer 
whose tube connection terminates at the rear support tube ( Figure 13) , in the 
corner between the support tube and the sphere. 

The flow picture of the supercritical state ( Figure 10b) approaches the 
ideal picture more closely, because the rear side turbulent region is considerably 
smaller in comparison with that for the subcritical state. In the same manner as 
in model -airplane airfoils, the transition from subcritical to supercritical flow 
occurs when the boundary layer becomes turbulent and energy of the external 
flow is transported to the layer near the wall and decelerated by friction and thus 
in part overcomes the pressure increase at tl?e rear side, so that the flow 
remains largely attached, and the resistance tecomes considerably smaller, 
which is the decisive factor. Supercritically, the drag coefficient is only %: 

subcritical c = 0.48 , 
w 

supercritical c^ =0.08 , 

based on the measurements of the DVL by Hoerner [9] in turbulence-free air 
( Figure 11) . 

The turbulent transition does not occur suddenly but, as with the flat 

plate, in a transition region in which that Reynolds number at which c =0.3 

w 

is considered the critical Reynolds number of the sphere. Twenty years ago, 
considerable trouble was caused for research laboratories when each one found 
a different critical number for the sphere, until Professor Prandtl explained 
the reason: the transition must occur earlier the more turbulent is the airflow 
of the wind tunnel in question, and the critical Reynolds number of the sphere 
can thus serve as an indication of turbulence. According to Hoerner, the transi- 
tion occurs for the sphere towed in still air at Re^ » 405,000, and in gusty air 

in areas surrounded by trees slightly earlier, at about 380,000. Thus in regard 
to its influence on the boundary layer, free atmosphere can be considered free of 
turbulence; the fine structure of vorticity of a boundary layer only millimeters 
thick is not effectively related to the large turbulence nodes of the free atmo- 
sphere. Strictly speaking therefore, only measurements in turbulence-free air 
correctly correspond to free flight. This is especially true for the measure- 
ments involved here at low Reynolds numbers, because in turbulent flow the 
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Figure 11, Measurement of the critical number for a sphere 

at c =0.3 serves as an indication of turbulence of the airflow 
w 

in wind tunnel. In turbulence-free air there was measured: 

Re « 4.05 ♦ 10 5 . The ratio of this number to the number for 
k 

a sphere measured in the wind tunnel gives the turbulence fac- 
tor of this wind tunnel. For the NACA variable-density tunnel 
TF = 4.05:1.55 « 2.6; this corresponds to the sphere mea- 
surement with a turbulence grill in Figure 9. For the Cologne 
wind tunnel after rebuilding, TF was 4. 05:3. 8 = 1. 06. 


subcritlcal state generally does not appear and the maximum value of lift is 
completely in error. Therefore all the previously measured polar curves in 
turbulent flow are worthless for model flying work, especially the American 
NACA airfoil measurements which were made in a strongly turbulent variable- 
density tunnel down to Re = 42,000. 


In the wind tunnel used for the following measurements the critical 
transition initially occurred at Re fc = 207,000. Only after a year of effort did 


improvements in the utility of the wind tunnel occur, allowing the subcritical 

state to become measurable on the model-airplane airfoil and rewarding the 

laborious preliminary efforts; after rebuilding, a critical sphere characteristic 

of Re = 380,000 was measured. This favorable result was achieved by: 
k 






X. incorporation of screens between the straighteners an- the nozzle, 
to convert the coarse turbulence into fine turbulence, 

2. by a long stilling run between the last screen and the nozzle, to 
cause attenuation of the fine-structure turbulence, 

3. incorporation of a narrow nozzle with a pronounced contraction, 

4. change of diffuser and blower. 

It was now easy to make virtue of necessity by using a turbulence grille 
( Figure 12) to make the flow completely turbulent, that comparison of measure- 
ments with and without the grille ( Figures 9 and 82) would show the influence of 
tunnel turbulence and in addition allow test of the effectiveness of measures for 



Figure 12. Turbulence grille before the nozzle of a wind tunnel, 
to make the airstream before the airfoil model completely 
turbulent. 

artificial turbulence, as described in Subsection VI 4. At this point, the suc- 
cessful attempt of Prandtl and Wieselsberger should be referred to; in this, a 
wire hoop placed around the upstream side of the sphere ( Figure 13) will induce 
the supercritical state at smaller Reynolds numbers, proving that it is merely 
necessary to make the boundary layer turbulent. A rough surface achieves the 



Figure 13. Fundamental experiment of Prandtl by which he 
proved that making the boundary layer turbulent by means of 
a wire hoop can compel supercritical ( attached) flow around 
a sphere at low Reynolds number, [from Prandtl-Tietjens , 

A eromechanik (Aeromechanics) , vol. 2, Springer-Verlag, 
Berlin] 


same effect ( Figure 11) . In addition, it may be concluded from this that 1. the 
sphere should be as smooth as possible, 2. the air stream before and after the 
sphere should not be disturbed by any kind of support wires or the like, and 
therefore the measurement is correct only when the sphere is held by a support 
tube at the back. 

In conclusion, reference should be made to determination of Re fc from the 

measurement of the pressure transition at the rear side of the sphere. This is 
considerably simpler than drag measurement, for which the proper determination 
of the amount of resistance to be subtracted for the support is laborious. With the 
second manometer mentioned, the pressure p prevailing at the rear side of the 
sphere can be measured. In the subcritical state, it is negative but becomes 
positive in the supercritical state, as shown by the plus sign in Figure 10b. 
Pressure measurement of this type is shown by plotting the dimensionless 
coefficient p/q^ versus Re, as Figure 14 shows schematically. That value at 
which p/q^ = 0, that is, at which the manometer on the rear side of the sphere 
shows p = 0 upon the change from minus to plus, serves as the critical Reynolds 
number. The results of these measurements were briefly described in 
Subsection II 2. 
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Figure 14. Pressure at the rear side of the 
sphere, plotted as pressure coefficient p/q^ 
versus Reynolds number. That Reynolds num- 
ber at which p/q^ passes through the zero line 
going from minus to plus values serves as the 
critical sphere characteristic; this agrees with 
the critical sphere characteristic determined 
from drag measurement at =• 0.3, 


c) Comparati ve Conclusions for Model -Airplane Wings 

1. Round-nosed, thick wings are sensitive to Reynolds number and turbu- 
lence; thin sharp-nosed airfoils are insensitive to Reynolds number and 
turbulence. 

2. The lower the Reynolds number, that is, the smaller the model airplane 
or its speed, the thinner must the profile be to achieve the supercritical 
flight state. 

3. To achieve the supercritical flight state, it is sufficient if the upper sur- 
face flow is turbulent. 

4. The critical Reynolds number of an airfoil sensitive to Reynolds number 
and turbulence can be reduced by artificial creation of turbulence in the 
upper-surface boundary layer, through 

a) pointing of the wing nose ( knife edge) , 

b) by a rough-surfaced wing nose, 

c) most effectively, by stretching a turbulence wire or thread parallel 
to the wing's leading edge, and in the wind tunnel by the use of a 
turbulence grille ahead of the model, 

d) in addition, the critical Reynolds number can be reduced by sonic 
vibrations, for example, by an intense whistle note. 


d) Convergent and Divergent Flow 

The conditions under which the boundary layer on the flying model air- 
plane is laminar or turbulent are supplemented by the following study. Repro- 
duction of the flow around the model airplane in the wind tunnel assumes that the 
influences stemming from turbulence of the wind stream are eliminated to as 
great a degree as possible, because wind-tunnel turbulence falsifies the entire 
picture, as Figures 9 and 11 show. An effective means for reducing tunnel 
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turbulence is, as is known, the constriction of the wind stream by a nozzle 
(see [5] , vol. IV, part 2, p. 74) , This can be explained approximately as 
follows: In addition to its motion in the direction of the principal flow, the small 
vortices of the turbulent flow cause disturbance motions directed transverse to 
the principal flow; these are damped by the constriction of flow. Conversely, 
widening of the flow causes increase in the disturbance motions of the vortices 
directed transverse to the principal motion, and thus produces an increase in 
turbulence. It has been observed in flow in tubes that a very small convergence 
considerably increases the critioal Reynolds number, and a small convergence 
decreases it (see [4] , vol, 2, p. 66) . The same is true for flow around the 
airfoil, as will be described in the next section. 


5. Flow Around the Airfoil 

The airfoil positioned for lift ( Figure 17) causes a downward deflection 
of the air mass involved, and this in reaction produces lift. As long as the 
flow is attached to the profile on all sides, the detailed shape of the airfoil's 
upper and lower surfaces is less Important than the course of the profile mean 
line as the governing, central streamline. Since, however, in the range of low 
Reynolds numbers various separation phenomena depend upon the details of 
profile form - as a cause of the drag -lift ratio decrease - several of these 
phenomena will be described here. 


a) Streamline Form and Bernoulli's Equation 

The streamline pictures ( Figures 15 and 16) made visible in flow 
apparatus clearly show the ohanges in convergent and divergent flow, and the 
separation phenomenon characteristic of the latter, as angle of attack increases. 
If ttie profile shown in Figure 17 is considered to be a uniform cross section of 
an airfoil of infinite span, then spanwise there is no flow change. At two cross 
sections separated by distance b in the span direction the streamline pictures 
are congruent ( two-dimensional flow around an airfoil) . If "a" designates the 
distance between two streamlines ( Figure 17) , then the flow tube delimited by 
it has the cross section f = a w • b in the undisturbed flow far from the wing. 

If the velocity changes from v m to v along this flow tube, then there results for 
an incompressible fluid, from the continuity requirement, 

a • b • v = a • b • v 
00 00 
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Figure 15. Attached flow around an airfoil at a low angle of 
attack. On the upper surface of the wing, however, at the half- 
chord point, the thickening of the boundary layer can be seen, 
and at the bright points the deceleration in the boundary layer 
is the beginning of the inverse flow and subsequent separation, 
[from Prandtl, Abriss der Stroemungslehre (Outline of Flow 
Mechanics) , Verlag Vieweg & Sohn, Brunswick]. 



Figure 16. Breakaway of flow on the upper surface at a high 
angle of attack ( stalled flight condition) . [from Prandtl, 
Abriss der Stroemungslehre ( Outline of Flow Mechanics) , 
Verlag Vieweg & Sohn, Brunswick ] . 
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Figure 17. Streamline picture of the wing during lift. 


In other words, the velocity is inversely proportional to the spacing of 
the streamlines. On the upper surface of the airfoil the streamlines are closer 
together, and on the under surface farther apart, than upstream from the airfoil: 
a Q < a B < a^: the approach flow runs convergent to the airfoil upper surface 

and divergent to the airfoil under surface. Convergence of streamlines denotes 

acceleration, and divergence deceleration, because v > v > v . In addition, 

o 00 u 

along a flow tube the velocity v and static pressure ( wall pressure) p are inter- 
related by Bernoulli's law: 

P +-f v 2 = Pgo + | yj = P = const. 

P + q = P* + q,* 

In other words, the sum of the static pressure p and dynamic pressure q 
is constant. The dynamic pressure of the external flow is easily measurable 
by means of the dynamic indicator discovered by Prandtl. There then results 
from the dynamic pressure, v m = *J(2/p) q^, or with p » % for air at 0 m 
height, v w - W q . °° 

In the study of flow pictures, the following holds: 

Convergent flow: decreasing streamline spacing means acceleration and 
pressure drrp; it acts in the sense of stabilization of laminar flow or reduction 
of turbulence. 
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Divergent flow: Increasing streamline spacing means deceleration and 
pressure increase; it acts in the sense of transition to turbulence, or, if turbu- 
lence already exists, intensifies it. 

In Figure 17 the streamline running to the stagnation point divides the 
flow into the flow on the upper and lower surfaces. The upper surface flow at 
first is convergent, especially over Ihe nose of the airfoil, and it therefore 
produces an intense negative pressure which forms a principal part of the lift. 
From about the highest point of the upper surface camber, that is, from the 
pressure minimum, the flow is divergent to the trailing edge. This is the criti- 
cal region of the profile, because there deceleration and excessive pressure 
increase can cause separation of the flow. On the lower surface, the approaching 
flow shows divergence and therefore pressure increase up to the stagnation point; 
from the stagnation point to the trailing edge there is a slight convergence and 
therefore decreasing pressure. At the trailing edge, there is a pressure 
equalization of the two flows. In reference to the effect on the boundary layer 
state, it must now be concluded from Subsection II 4d that along profile runs 
with convergent flow, that is, on the lower surface from the stagnation point to 
the trailing edge, and on the upper surface from the stagnation point to the pres- 
sure minimum, there is a tendency to maintain the laminar boundary flow. The 
divergent flow, on the other hand, tends toward formation of turbulence. The 
latter is favorable for the separation region on the airfoi 1 upper surface (that is, 
from the pressure minimum to the trailing edge) for the model airplane, 
because the turbulent boundary layer can better overcome the pressure increase 
in this region than can the laminar, and thereby the separation does not occur 

until a higher angle of attack, and therefore higher c values are achieved. 

ft max 


b) Pressure Distribution Along the Chord 

If the upper surface of the airfoil model is provided with small holes 
(0.4 mm diameter) then the effective wall pressure p can be measured at a 
hole if the hole is connected with a manometer by a fine tube ( through the 
interior of the airfoil or from the other side of the airfoil) . The plot of the 
measured wall pressures as a function of the chord then gives the curve of 
"pressure distribution" along the chord. In Figure 18, the pressure distribution 
is plotted for a thin profile and in Figure 19 for a thick profile, from 
Goettingen measurements ( see [3] , II, pp. 44 and 46) . Not only is the wall pres- 
sure p plotted, but also, as customary, the dimensionless pressure coefficient 
p/q . Because the pressures perpendicular to the chord are plotted, file 
resultant pressure area from the upper and lower sides corresponds to normal 
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force N or to the coefficient of this aerodynamic force c^. * ) The velocity 

distribution obtaining at the free edge of the boundary layer can be determined 
from the pressure distribution: 

v = v „o - (y-) = n/-j - <±p) 1 w 4 - (±p) » 

in which negative values of p give an increase of v compared with v M . At the 
stagnation point, where p/q^ =* 1, v becomes 0; at pressure equalization where 
p/q^ a 0, v is v B . Strictly 0 speaking, the Bernoulli equation is valid only for 
turbulence-free flow, so that it cannot be applied to profile elements with 
separated flow. 

At low angles of attack, the pressure distribution is fairly similar for 
both profiles. As the angle of attack increases, there forms on the thin sharp- 
nosed profile ( Figure 18) , as a result of a local transition vortex, a "suction 
peak" which at a high angle of attack lies exactly at the foremost profile point. 

At the suction peak, p/q^ = -3.7; the pressure minimum therefore is 3.7 
times the dynamic pressure q^. For the thick, round-nosed profile ( Figure 19) , 
on the other hand, the pressure distribution is "more complete"; the pressure 
minimum extended only to x = 0. 2 t of the leading edge (= smaller movement 
of the center of pr essure) , an d it is p/q^ = -2. 4. At the suction peak in the 
thin profile, v = *>/l - (-3.7), which is 2.16 times the speed of the outer flow; 

at p . of the thick profile, v = 1.85 v . 

*min 00 

The pressure distribution surfaces give an important idea in regard to 
the relation of the airfoil upper surface and under surface in the production of 
lift. The negative pressure areas in Figures 18 and 19 are five times as great 
as the positive pressure areas ( - plotted above , and + under the zero line) , so 
that an airfoil can be rightly called a "suction wing. " 


* ) The coefficient c determined on the same airfoil in a wind-tunnel 
n 

measurement is somewhat larger; the difference corresponds to the friction 
drag, which does not appear in the pressure measurement. 




420,000 










c) The Reynolds Number of the Boundary Layer, Re^ 

The pressure distribution according to Figures 18 and 19 shows the wall 
pressure and thus the pressure state in the boundary layer. The local velocity 
v calculated from it is valid only to the edge of the boundary layer; however, 
while in the boundary layer itself the velocity gradient is determined by the 
laminar or turbulent state ( Figure 4) . If the local velocity v effective at the 
edge of the boundary layer is now used, along with the friction length t as the 
distance from the stagnation point, to form a Reynolds number of the boundary 
layer Re fi » v • t/v, this number differs from the Reynolds number of the 

plane friction plate of Figure 6 only in that now v changes with the pressure 

state. At the stagnation point, this Reynolds number is zero, and at the trailing 

edge it is approximately equal to the profile Reynolds number. Because on thj 

lower surface the velocity v, proceeding from a zero value at the stagnation 

point, first reaches the value v^ of the external flow at the trailing edge, the 

increase of Re. to Re on the lower surface occurs more slowly and on the upper 
o 

surface more rapidly, than the increase of the Reynolds number along the plane 
friction plate. Accordingly, the laminar approach stretch on the lower surface 
of airfoils is longer, even with a flat under surface, and shorter on the upper 
surface, than on thin plane friction plates. 

To obtain a picture of this, the Reynolds numbers of the boundary layer 
were each determined for two angles of attack from the pressure distributions 
given in Figures 18 and 19, and the ratio Re^/Re wa - plotted in Figure 20 as a 

function of the chord t. The diagonal represents the local Reynolds number for 
the flat friction plate. As a consequence of circulation (acceleration on the 
upper surface, retardation on the lower surface) the Re number line for the 
suction side flow goes above the diagonal and that for the pressure side below 
the diagonal. It is immediately seen that this representation clearly shows the 
thickening of the boundary layer, or the separation point, on the upper surface 
as a sudden increase or bend. Although on the thick profile 382 the separation 
point changes only slightly with angle of attack change from 5. 6 to 14. 5® , with 
the thin profile and for the same a -range there is a pronounced migration from 
0. 75 t to 0. 3 t. From the higher velocity v in flow around the nose of the thin 
profile, a higher local Reynolds number would be expected here along the upper 
surface than for the thick profile. Instead of this, however, Figure 20 shows 
higher local values of Re jRe for the thicker profile. 


The conclusions from pressure distribution (Re = 420,000) are: 

1. The sharp leading edge gives a high suction peak but the large pres- 
sure gradient is countered by an equally great pressure increase; in 
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Figure 20. Evolution of local Reynolds number Re^ on the free 

edge of the boundary layer for profiles 389 and 382, shown in 
ratio to the Reynolds number of a friction plate of length t = t. 
The diagonal corresponds to the Reynolds number of the friction 
plate. The Reynolds number of the pressure measurement used 
for calculation of v is 420,000. 


consequence, the pressure conversion and the flow around the nose 
involve an energy loss which is greater as the leading edge is 
sharper; as a high angle of attack and high Reynolds numbers fills is 
unfavorable, because reduced energy content of the boundary layer 
and higher pressure gradient effect an earlier separation of the flow. 
In the large airplane, accordingly, the nose should not be pointed but 
should be kept as round as possible to achieve a high c ft max value. 

2. If the surface friction were to be regarded as file sole turbulence 
cause for the upper-surface flow, then the transition would have to 
happen sooner for the thick profile than for the thin, because in spite 
of higher acceleration on the flow around the nose the local Reynolds 
number is higher than for the thin profile. The turbulent transition 
leading edge is based on another cause, however, which justifies the 
sharp nose on the model airplane in spite of the above-mentioned 
losses. 
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d) Movement of the Stagnation Point 


With the increase in circulation, the stagnation point S t moves counter 
to circulation on the rotating cylinder in parallel flow; on the other hand, the 
rear stagnation point S 2 , which in the transition from ideal to friction flow 
becomes the separation point A, moves with the circulation direction 
( Figures 21 and 22) . At u/v “4,8! and A have approached each other to such a 
degree that they coincide (see (4], vol. 2, p. 97) ;the body has been drawn com- 
pletely into a negative pressure region. The streamlines move on the leading 
side toward the point of maximum negative pressure and on the rear side 
divergently with increasing pressure. As soon as the circulation becomes 
smaller, the convergence decreases: the stagnation point moves back. Phenom- 
ena on the airfoil are similar to this. The theory bases the origin of lift on the 
airfoil in the effect of the superimposition of a parallel low on a circulatory 
flow; their velocities are additive at the airfoil upper surface because there they 
are in the same direction and thus produce an acceleration and accordingly a 
negative pressure, and on the other hand they are subtractive on the airfoil under 
surface because they are counter to one another; that is, they produce retarda- 
tion and positive pressure. Lift is proportional to the circulation. Analogous 
to the process on the cylinder, the stagnation point on the airfoil moves counter 
to tiie circulation, that is, from the forward profile point along the lower surface. 
At a high angle of attack and high lift, the suction effect of the upper surface is 
several times greater than the stagnation pressure on the lower surface, so that 
with the wing the flow ahead of the wing is convergently drawn up, so that the 
streamline striking on the stagnation point is bent sharply upward ( Figures 15 
and 17) , If the flow detaches, then with decreasing circulation, the convergence 
of flow and the stagnation point backward displacement become smaller. *) 

Even though the movement of the stagnation point is only a few percent of the 
chord, nevertheless this has a considerable influence on production of turbulence, 
as will be shown. With increasing angle of attack, the stagnation point moves 
back, as far as t/2 at a = 90®. 

To obtain an idea of the movement of the stagnation point from the pres- 
sure distributions given in Figures 18 and 19, the pressures measured at the 
five holes farthest upstream are plotted on the development of the profile nose 
in Figures 23 and 24. The thick profile has the higher stagnation point move- 
ment. For an a-range of about 15® on the thick profile, this is about 4 percent 
of the chord, and on the thin profile 2. 5 percent. If the wing nose is replaced 
approximately by circles, then the center angle traversed by the stagnation point 
for this a -range of 15® is about equal to 90® for the two profiles. 


*) The phenomenon of hysteresis can be explained by this. 


40 


-^^AiaagMl»iiij l | ' l i. l . l l '.'Hi—WM 




Figure 21 


Figure 22 



Figure 21 . Theoretical flow around a cylinder ( or a sphere) in 
an ideal friotionless fluid (potential flow) . 

Figure 22. Theoretical flow around the rotating cylinder in a 
parallel flow. The front stagnation point moves counter to the 
circulation and the rear one with the circulation. 


e) Change in Pressure and Velocity in Flow About the Nose 

In addition to center -of-pressure movement, Figures 23 and 24 show the 
change in pressure in the vicinity of the stagnation point. The maximum pres- 
sure gradient can be seen on the thin profile at a = 14.6°, where the stagnation 
point coincides with measurement hole 9, but the pressure mi nimum p = 3. 7 
coincides with hole 0 at the leading edge 4 mm away ( Figure 23) . Accordingly, 
from the stagnation point to the leading edge the total pressure drop is 4.7 q^. 
Because the pressure measurement was made at v = 30 m/sec, at the leading 
edge ( according to Subsection n 5b) the velocity v = v w \/l - ( -3.7) = 30 • 2.16 
w 64 m/sec prevails. This velocity develops at the edge of die boundary layer 
along the short stretch of i = 4 mm, so that the acceleration from v = 0 at the 
stagnation point to v = 64 m/ sec occurs in the time of 

t « 2 l/v = 2 . 0.004/64 = V 8000 sec. 

The acceleration, assumed uniform, is accordingly 

b = v/t = 64 • 8000 « 512,000 m/sec 2 . 

This means that if the acceleration were to last only one second the 
velocity would be 512,000 m/sec; however, since it lasts only V 8000 sec, v reaches 
only 64 m/sec. Close behind the stagnation point, however, the acceleration is 
even greater, as the slope of the osculating tangent shows. This acceleration is 
several times greater than that of a rifle bullet, which, for example, at the 
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muzzle of an 0. 8~m long rifle barrel has reached a velocity of v = 800 m sec in 
Vm sec an ^ th us achieved an acceleration of 200,000 m/sec 2 . 

On the thick profile 382 and at a = 14.5° the velocity increases from the 

stagnation point out up to the leading edge from 0 to 34 m/sec along a distance 

of 7 mm in Vi>4oo see * Acceleration here is 82,000 m/sec 2 , and thus less than 

V 6 that of the thin profile. The pressure minimum is farther above the back of 

the nose, however, and v =55 m/sec is not reached until after 38 mm path 

max 

length from the stagnation point, while on the thin profile the velocity decreases 
from directly above the leading edge. Accordingly, the flow on the thin profile 
is more likely to be turbulent sooner and separate sooner, as will be explained 
later. 


f) The Transition Vortox 

The suction peak and the velocity connected with it in the flow around the 
wing nose are higher as the profile is thinner . In the thin, plane plate with 
sharp leading edge, therefore, these phenomena must be especially prominent, 
so that the flow around the flat plate will be studied somewhat more closely. 

At a = 0° angle of attack on the thin plane plate only friction drag is produced, 
and at a = 90° only pressure (form) drag, while at an angle of attack as an 
airfoil, lift and both kinds of drag are produced. The increase in lift per degree 
of change in angle of attack even considerably exceeds that of the customary 
profiles, although of course only to about 10 s , because then the flow has 
separated. The separating tendency is strongest at a = 90* . On the plate to 
which flow approaches perpendicularly ( Figure 25) , in the first instant of 
beginning flow the flow follows the rear side of the plate ( potential ow) . At 
increasing speed, the particles of the boundary layer flowing on the front side 
toward the edge attempt to go out transversely to the principal direetion over the 
edge; they are diverted by the neighboring layer of fluid and forced to change 
direction by 180 4 . A point of mass p, which moves at velocity v and radius r 
on a curved path, has a centrifugal effect 



If the edge is knife-sharp (radius r — 0) then the centrifugal force in a 
frictionless fluid would have to be infinitely large and thus produce an absolute 
vacuum, and the flow around the edge therefore would have to occur at infinitely 
high velocity. Instead of this, however, the boundary layer detaches itself as 
a "separation surface" from the edge, forms a small vortex, which quickly 
becomes larger ( Figures 26 and 27) and then is carried away by the flow to 
make space for a newly forming vortex. 


DIRECTION OF PRINCIPAL FLOW 



Figure 2D, Starting Figure 26. Development Figure 27. Detaehed se pa- 
sta to: potential flow. of a curling vortex. ration surface: progross- 

ing vortex. 


Figures 2D-27. Flow around the edge of a plate perpendicular 
to the principal flow. 


If the plate is set, as an airfoil, at an angle of attack from 0 to 3® , for 
example, then a substantial lift is effective, with the flow ahead of the plate 
acting convergently and the stagnation point moving backward, so that the 
boun da ry layer must flow from the stagnation point around a sharp leading edge 
with a directional change of almost 180® ( Figure 28) . Similar to what has been 
described above, there then forms over the leading edge a small eddy as a 
small local separation, which at first does not lead to separation of the flow, 
because on the rear side of the eddy the boundary layer again attaches itself 
turbulently ( Figure 29) . Because the boundary layer detaching itself in 
laminar form above the leading edge becomes turbulent through this eddy, it is 
called the "transition eddy." As the angle of attack increases, the vortex 
broadens and flattens itself, and the external flow forms a bridge over it, while 
a profile with a fluid upper side forms to a certain extent ( Figures 29 and 30) . 
With increasing angle of attack, the bridge span approaches the trailing edge 
more closely, and the character of the vortex changes; it is no longer local but 
increases. The change is from the stationary transition vortex to the curling 
vortex, or from the "fluid profile" to the "separation surface. " Individual 
vortices under the bridge move backward ( Figure 31) , and the flow separates 
for an instant but again attaches after passage up to the next vortex. This 
phenomenon manifests itself unpleasantly as heavy vibration on a plate hanging 
in the wind tunnel. After a small increase in the angle of attack, the flow then 
remains separated, the vibration of the model decreases, and it then hangs 
again almost quietly in the wind tunnel. 
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Figure 28 


Figure 29 
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Figures 28 and 29. Flow around the leading edge of a flat plate 
at low ~ngle of attack; origin of the turbulent boundary layer on 
the upper surface by stationary transition eddy. 


GREATLY THICKENED BOUNDARY 
LAYER 


CURLING VORTEX 




Figure 30. Extension of the transi- 
tion eddy with increasing angle of 
attack. 


Figure 31. Transition from stationary 
transition eddy to periodically separat- 
ing curling vortex. 


The negative pressure at the leading edge acts as a "suction force" 
counter to the drag component; theoretically both cancel out in potential flow of 
the ideal fluid with circulation. The drag would then equal zero. 

The phenomena were first described in principle by Prandtl ( see [2] , 
p. 145) in connection with a discussion of the behavior of the flat plate in a 
potential flow with circulation: "In the actual flow around a plate with a sharp 
leading edge, infinitely high negative pressures do not occur, of course, but 
there is a separation of the flow; this at small angles of attack through the 
action of a turbulent process again attaches itself to the upper surface so that a 
flow form not entirely dissimilar in general terms to the theoretical occurs, 
with a similarly high lift. In the absence of the suction force, a resistance is 
naturally obtained here which exceeds the friction drag by A • tan a ( A = lift) . 
The equivalent of this additional drag is to be sought in the loss in speed which 
occurs through the turbulent process on the upper surface. " In the same place, 
it is said concerning the significance of the suction force: "With well-rounded 
profiles the suction force takes on practical meaning. If the resultant of an 






airfoil ino lined at a 6° angle of attack is only 2® to the perpendicular, as is 
frequently observed, then the effect of the suction force can be clearly seen. " 

If die velocity of sound ( c » 340 m/ sec) is achieved in the flow around 
the sharp leading edge at the free border of the boundary layer, then the 
mechanism of nose turbulence is influenced by Mach waves ( see [2] , p. 184) . 
Because the flow at this point occurs in a region of lower pressure, these must 
be rarefaction waves. It can be assumed that this phenomenon intensifies 
turbulence. Its range of influence must be very small, however, because the 
velocity decreases rapidly again to that of the outer flow 


g) The Upper-Surface Flow 

This begins at the wing nose at the splitting point of the flow; the stagna- 
tion point, and therefore, because of stagnation point migration on the airfoil 
undersurface of an airfoil at the angle of attack, so that the distance to the 
trailing edge for the upper surface flow in die thick profile 382 at 14. 5° , for 
example, is about 17 percent greater than the path along the lower surface. 

Four characteristic points characterize the flow state along the upper surface: 
stagnation point, pressure minimum, transition point, and separation point. 
After the pressure minimum, a pressure increase begins which, depending on 
the magnitude of the suction peak, can reach four times that of the stagnation 
pressure ( Figures 18 and 19) against which the boundary layer must flow 
and in the most favorable case overcome it as far as the trailing edge. Over- 
coming of this pressure increase is the most important phenomenon of flow 
around the airfoil. The closer the separation point can approach the trailing 
edge the greater will be die lift achieved, the less the resistance, and thus the 
more favorable the drag-lift ratio. The theoretical maximum would be 
achieved by potential flow, that is, by vortex free and frictionless flow around 
the airfoil in an "ideal" noncompressible fluid, as can be calculated scientifi- 
cally by extensive support of mathematics. Through superposition of two 
potential flows - a parallel flow and a circulation flow - a lift, but no drag, can 
be demonstrated by calculation for an airfoil with an infinite span. This would 
be the most favorable limiting case of the drag-lift ratio c ^/c a = 0/eo) - 0 » 

that is, practically a gliding flight without loss of altitude in still air. In 
reality, however, even with a sound flow, that is, one attached all the way to 
the trailing edge, because of friction and vortex losses, a drag-lift ratio of 
1:30 today must be looked upon as the maximum value. From a theoretical 
potential standpoint, the airfoil lift increases with the angle of attack according 
to the equation 

c = 2 7 t • sin or , 
a 
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so that at a « 90° , a c = 2 ir » 6. 28 would be achieved. Since, however, 

’a max 

the actual flow depending on profile form and Reynolds number detaches pre- 
maturely supercritically in the range of a - 10° to 30° on the wing upper 

surface only c » 1 to 2 is reached; with the aid of flaps 2 to 3 is reached, 

J a max 

with suction airfoil » 4, and with an airfoil out of which air is blown, > 4. A 
brief description will be given here of the three types of force which operate on 
the straight airfoil in straight flight: 

1. Pressure forces act on all particles in the pressure field of the 
airfoil. 

2. As a result of the viscosity of air, friction forces act from the wall 
on the particles near the wall and thereby produce a friction layer. 

3. Drag forces from the outside flow act on the friction layer when it 
flows turbulently. 

Considerable clarification is still awaited on the details of the manifold 
interplay of flow phenomena on the airfoil. A good idea of individual processes 
is obtained from flow photographs and measurements made in a closed water 
tunnel of the Aerodynamic Research Laboratory at Goettingen by P. Jordan 
( see [12] , p. 191) . These tests were done on the symmetric Goettingen profile 
409 at Re » 150, 000, and thus still in the region of model flight. They are 
therefore especially worthwhile for the study at hand. Profile 409 ( see [3 ] , I, 
p. 94, and HI, p. 29) has a thickness d = 0.127 t; the radius of the nose circle 
is r « 0. 017 t. The test model with an aspect ratio of A = 2. 5 was bounded by 
flat walls at the airfoil ends to realize two-dimensional airfoil flow as much as 
possible. To make visible the phenomena in the boundary layer, colored liquid 
was supplied out of six holes, 0.44 mm diameter, in the middle section of the 
wing. * ) At the negative angle of attack a » -3°, Figure 32a shows die laminar 
flow along the lower surface; although at -4, 5® it probably is completely attached 
to the airfoil, at -3® the already existent pressure increase cannot be overcome 
and consequently causes separation at about 3 / 4 t. The outflow point of the 
colored liquid at hole 6 in Figure 32b is shown by the arrow; the reverse flow 


*)On the side walls a secondary flow (see [3], I, p. 89 and [6], p. 103) 
influences the flow; this is less noticeable in the center of the wing the higher the 
aspect ratio. The boundary layer flow can be seen most completely in the closed 
water tunnel with dyestuffs emitted from the airfoil; in addition, higher speeds 
can be reached than in the open water tunnel where surface tension and, at high 
speeds, wave formation, influence boundary pehnomena. 


/wedge according to Figure 6 is clearly visible. From this phenomenon, it can be 
concluded that it would be inadvisable to camber the airfoil of an airplane model 
substantially on the underside as is done on modern profile sections for large 
airplanes. 

The laminar separated boundary layer is mobile, as is any separation 
surface. Through the sudden increase in speed at the separation surface, it 
tends to form waves whose swells then cause vortices to roll up ( Figure 33) 
and finally to break down in disordered succession of vortices: the laminar 
separated boundary flow becomes turbulent. At a = + 1 0 , the upper side has 
become the suction side and the stagnation point has begun its movement to the 
lower surface so that the turbulence effect of flow around the airfoil nose 
becomes effective. At 0.7 t a laminar separation can be observed ( Figure 32c) ; 
the laminar separated boundary layer becomes turbulent, however, after a 
distance of 0.2 t, and at 0.9 t again attaches itself, because the impulse exchange 
occurring from the external flow just suffices to overcome the pressure increase. 
The separation, which is only local, can be seen through its shadow in the 
enlarged Figure 32d. 

Jordan states further, in a note, that separated flow can again attach 
itself for other reasons. "If, for example, the wall forms a reentrant corner 
( Figure 34) then in this corner there occurs a local dead space containing 
vortices. At the site where the flow again attaches itself after this dead water 
a stagnation point occurs. The boundary layer proceeding from this stagnation 
point then begins in laminar form. " This is an interesting example of the 
reduction of turbulence as a consequence of convergence of flow. 

The transition location clearly indicated for a = +1° contracts with 
increasing angle of attack and moves toward the wing nose. At a = 4.5° in 
Figure 32e, laminar separation can no longer be seen. It can be concluded that 
this phenomenon is possible down to molecular size; probably every turbulent 
transition point on a wall lies behind a laminar separation point. The further 
observation that at a = 10. 5° the transition point extends almost to the leading 
edge ( Figure 32f) , and here the apparently laminar separation and the turbulent 
transition are again clearly separated, as can be seen from the enlargement 
( Figure 32g) , speaks further for this assumption. The colored liquid emerges 
at hole 1, indicated by the arrow, and thus flows, counter to the main flow, 
toward the leading 1 edge. Between the wing and the outer flow above the airfoil 
leading edge a transition eddy acts; at its rear side the boundary layer again 
attaches itself turbulently. It is a phenomenon similar to that at +1°. Behind 
the transition vortex the boundary layer is of course much thicker, which per- 
mits a conclusion of an early breakaway. This transition point is not visible in 
the picture, because the nearest colored liquid hole 3 lies farther back, at the 
righthand arrow. The scheme is shown in Figure 35. 
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For Figures 32a through 32g: laminar and turbulent 
phenomena in the boundary layer, made visible in 
the water tunnel by colored liquid, which emerges 
from holes in the airfoil. 



3° • Dyestuff emission at the arrow to ind 
reverse flow in the reverse 






f) a = 10. 5° • The transition point has reached the vortex, in front the laminar separation, transition 
leading edge and expands as a transition vortex * region not visible* Second dyestuff hole shows a turbu- 
Turbulent separation at 0* 9 t* lent thickened boundary layer* 






Conclusions: 


LAMINAR SEPARATION 


TRANSITION 


The transition point of laminar 
into turbulent boundary layer 
flow moves along the upper 
surface from the trailing edge 
to the leading edge: 


a) with increasing Reynolds 
numbers in small degree 
through friction turbulence, 

b) essentially, however, with 
increasing angle of attack 
as the result of "nose 
turbulence, " through the cooperative effect of stagnation point move- 
ment and acceleration at the flow around the nose, and then further 
by flow divergence as a consequence of the upper-side curvature 
and angle of attack. 

The transition point reaches the leading edge earlier, that is, at smaller 
Reynolds numbers and smaller angle of attack, 

a) the smaller is the camber of the upper side, and therefore earliest 
for the flat plate. 

b) the sharper is the leading edge. 

The suction peak manifesting itself on thin or slightly curved profiles 
with increase in angle* of attack (cf. Figure 18 at 11.6° and 14. 6°) shows 
that the transition point has reached approximately the foremost position 
and that its broadening into a transition vortex is occurring. 

In general, the transition point lies behind the pressure minimum. For 
a high-speed airplane which flies in normal flight with a very small c a , 

and therefore uses a thin, approximately symmetrical profile, it is 
therefore advantageous to move the maximum camber back so that the 
transition point lies far downstream, to achieve as long as possible a 
laminar run of low surface friction: "laminar profile. " In contrast to 
this, a "turbulent profile" with a high mean camber but with a smaller 
rearward displacement of the maximum camber and a sharp leading 
edge is favorable for the airplane model, as bird wing profiles show. 



4, \l 


Figure 35. Scheme of the transition 
vortex. 




5. 


For the thin flat plate at the minimum angle of attack, the nose turbulence 
effect is very great and the pressure increase is small toward the trail- 
ing edge; a small energy oontent of flow at Re < 20,000 suffices for the 
boundary layer to attach turbulently and to displace the transition vortex 
to the leading edge, in contrast to the round- nosed, thick, high -camber 
profile where at low angle of attack the nose turbulence is small and the 
pressure increase is great as a result of high wall curvature, so that it 
is not until Re « 100,000 that the transition is forced; then, hov/ever, the 
attachment occurs suddenly along a considerable profile stretch. 

6. Because the effect of nose turbulence increases fester than pressure 

increase for the thin plate, the approach flow in the wind tunnel can 

begin at any angle of attack; in addition, the flow condition is reproducible 

with increasing or decreasing speed and at increasing or decreasing 

angle of attack, since the measurement then gives the same course of 

aerodynamic force plots, as a proof of "lack of sensitivity to Reynolds 

number. " The critical angles of attack a. nndc always are the basis 

k a max 

for the influence of Reynolds number, since both increase further with 
increasing drag effect of the external flow, and therefore with increasing 
Reynolds number, as is valid for all profiles in the Reynolds number 
range studied. 


h) Separation Phenomena of the Upper -Surface Flow 

The separation phenomenon in the presence of a transition vortex, which 
has been briefly discussed for the flat plate, may be described in further analysis 
of the Goettingen measurements of Jordan on profile 409. At first on profile 
409, intermittent separation and attachment of flow occur at a - 11. 25®. The 
^ oscillation of flow, which at first occurs slowly and then increasingly rapidly 
with increasing angle of attack, leads above 11.75° to a final separation. As 
shown by pressure recordings on a strip chart, the separations during oscilla- 
tion of flow begin at the leading edge; the transition vortex broadens into a 
bridge, then separates from the leading edge with periodically detaching vortices, 
and then again attaches itself. Although separation of the boundary layer pro- 
gressing from the rear forward was not seen here, it can be assumed that with 
a thicker profile having a sharper leading edge both separation phenomena can 
probably occur simultaneously. 

An instructive supplement to these flow observations is supplied by the 

pertinent measurements of pressure distribution on profile 409, which shows in 

Figure 36 the transition from the attached state at c through the oscillation 
b a max 

phenomenon to a fully separated state at measuring holes 1 to 6. The suction 
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peak is connected with the "attached 
state, and the negative pressure, 
increased with time, at hole 4 Is con- 
nected with the transtion between 
attached and separated flow during 
oscillation. Both have the same pres- 
sure area; but the "center of pressure" 
of the resultant aerodynamic force 
changes between */4 and t/3. After 
final break way, the same pressure 
obtains everywhere along the upper 
surface ( dot-dashed pressure area) , 
The reverse flow wedge has become 
a "dead water area" filled with vor- 
tices, extending to the leading edge 
and sucking reverse flow fluid around 
the trailing edge from the lower 
surface, so that the pressure dis- 
tribution on the lower surface is 
essentially changed, with negative 
pressure prevailing from about t/3 
on. As the angle of attack increases, 
the negative pressure again dis- 
appears on the lower surface. 

It may be said in summary 
about separation phenomena: 



Figure 36. Pressure distribution on the 

symmetric profile G 409 at the transition 

from c through the oscillation 
a max 

phenomenon up to completely separated 
upper-surface flow at 14.3®. 


? . In the subcritical flight state, the boundary layer remains detached in 
laminar form on the profile, and the separation point lies just behind 
the pressure minimum. 

2. Separation of the supercritical flow is preceded by a thickening of the 
boundary layer. 

3. The higher the Reynolds number the later does turbulent separation begin 
at the trailing edge; with increasing angle of attack and increasing profile 
thickness it then progresses to- the leading edge. 

4. At a low Reynolds number, a transition vortex is characteristic; it 
boradens with increasing angle of attack upon reaching the foremost 
position, and brief turbulent shocks give enhanced lift peaks, and then 
separation starts at the leading edge and proceeds along the entire chord. 


54 


If the leading edge is sharp and the camber of the upper surface small, 
then after separation the cause of turbulence remains, and accordingly 
the characteristic bridge building of the separated upper-surface flow 
exists: lift does not decrease but instead increases somewhat further 
with increasing « -value. 

On the thick round-nosed profile, the separation point passes the transi- 
tion point as angle of attack increases, so that a separation of flow 
occurs; the separation point then Jumps to tin point of pressure minimum 
corresponding to laminar separation. Lift decreases suddenly to its 
suboritlcal value and the flow then first increases suddenly at a consid- 
erably smaller a -value or at a higher Reynolds number (hysteresis loop 
in the polar plot) . 

5. The sharper the leading edge, the earlier does the transition vortex 
reach the leading edge, and the earlier doos its broadening begin; the 
flow separates at a smaller critical angle of attack the smaller is 

the achievable c . It can be concluded from this that the incorpora- 
a max 

don of a sharp edge on an airfoil leading edge is favorable only on a 
model airplane, to compel the supercritical flight stale, but on the large 
airplane, as angle of attack increases the drag-lift ratio worsens and 
premature separation occurs. 

To keep the picture simple, up to now only the two-dimenu ional flow 
around an airfoil has been considered. In the Influence of phenomena by span- 
wise flows, especially in separation pehnomena, the effect of the transition 
vortex on the twisted wing at a high angle of attack is worth mentioning in this 
respect. Here the flow which first separates in the center of the wing flows to 
the location of low pressure of the attached flow, while parts of the flow, using 
the transition vortex as a tunnel for transverse flow, again become boundary 
layer material, with the result that although the suction peak sinks there, a 
more complete pressure distribution (similar to that in Figure 36) occurs, 
permitting, for example, 2 to 4* higher angle of attack before the flow 

separates up to the wing tip. In addition, as a result of this, the nose-down 
moment is greater before separation, so that the model airplane tips more 
favorably forward instead of stalling. 

The flow picture of the transition vortex is similar to cavitation 
(cavity formation) on ships' propellers, where local water bubble vortices or 
vapor bubble vortices - on thin profiles going from the leading edge and on 
thick profiles from the pressure minimum - flow downstream and have a harm- 
ful action on the material of the propeller through powerful impact blows ( see 
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[5], vol. IV, parti, pp. 473-476, contribution by Ackeret) . Just as the cavity 
formation decreases with increasing external pressure of the water ( for 
example, on a diving submarine) , so with increasing Reynolds number a decrease 
in the transition number is probable. Broadening and separation then begin as a 
result of the higher drag, effect of the external flow at a higher angle of attack. 
Accordingly, the transition vortex has a greater importance in the Reynolds num- 
ber range of model flight than in the Reynolds number range of the large airplane. 
For ships’ screws it may further be concluded that measures to displace the 
beginning of broadening of the transition vortex to higher angles of attack are 
suitable for reducing cavitation, for example, by the use of profiles with 
rounded leading edge instead of the sickle profile with its sharp leading edge 
( Figure 56b) ; this profile's transition vortex and high suction peak are unavoid- 
able at a low angle of attack. 


i) The Flow Along the Lov/er Surface 


Although it suffices for achieving the supercritical flight state if merely 
the boundary layer on the upper surface is turbulent, it is still worthwhile to 
see why the boundary layer of the lower surface remains laminar in general, 
under what circumstances it becomes turbulent, and how the transition manifests 
itself. 

The profile underside can essentially be 

1. cambered concavely, as in the bird's wing; 

2. plane, as was first done by Junkers to achieve simpler manufacture; 

3. moderately negatively cambered as in the case in modern high-speed 
airplanes; 

4. made with S curvature, like the profiles of sailplanes. 

The velocity is reduced on the under surface of the airfoil (v^ < v w in 

Figure 17) , and the streamlines are convergent. Both phenomena are more 
effective in respect to maintenance of laminar state the higher is the angle of 
attack. Because nose turbulence is effective only on the upper surface, only 
friction turbulence will consequently be expected on a flat or moderately convex 
lower surface. 

The use of profiles with convex lower surface, especially symmetrical 
profiles, is not to be recommended for model airplanes, however, because at 
angles of attack in the range of normal flight the laminar boundary layer 
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separates from the lower surface, as is made clear in Figure 32a for the 
symmetrical profile G 409 at a = 3° and Re = 150,000. The separation dis- 
appears as the angle of attack increases. If this attachment suddenly takes place 
over a lengthy profile stretch, then the change in the values of c and c must 
show a sudden increase in the polar plot. w a 

For a flat lower surface, at transition from negative to positive lift, the 
profile is at an angle of attack where the pressure equalization p/q^ => 0 is 
achieved shortly after the stagnation point and, as for example in Figure 18 on 
profile 389 at a ~ 1* , it remains up to the trailing edge, so that on this stretch 
Re^ of plane friction must be exactly valid. Accordingly, the length i of the 

laminar initial stretch can be determined here. According to Figure 6, the 
beginning of the transition lies between Re = 500,000 and 1,000,000. For the 
airplane model, because of flight in the turbulent region near the ground and 
because of the greater relative roughness of the wing upper side, the smaller 
number is more likely to hold; the larger number is more probable for the 
large airplane. With Re^ = 500,000, the length of laminar initial run: 



is calculated for five examples: 



v, m/sec 

i, mm 

Model airplane 

6 

1200 

Model airplane 

12 

600 

Sailplane 

15 ' 

480 

Small sport plane 

40 

125 

High-speed engine- 
powered plane 

100 

70 


For Re^ = 1,000,000, the laminar initial run is twice as large in each 

case. In general, it is true that with increasing Reynolds number the laminar 
initial run becomes smaller; that is, the transition point achieved because of 
plane friction turbulence moves from the trailing edge toward the leading edge. 
In any case, the flow along the lower surface for model airplanes with profiles 
similar to that in Figure 17 remains laminar in the range of normal flight. 

The following concept is important for understanding of flow around an 
airfoil: between laminar and "fully turbulent" is a transition region which 
according to Figure 6 extends from Re = 5 • 10 5 to 20 times this amount or 
Re = 107, so that actually three flow states must be differentiated. Behind 



the above-mentioned transition points, the fully turbulent flow state would not 
be reached until after 20 times the friction length. It can therefore be concluded 
that behind the transition point a flow form corresponding more or less to the 
transition state prevails at first. 

Because in the usual aerodynamic force measurement in the wind tunnel, 
a is moved by steps to a positive angle of attack ( raised up) from ihe range of 
negative lift, and thus a turbulent state is initially present on the wing under 
surface, at a setting in the neighborhood of c & =* 0 the nose turbulence, which 

up to then influenced the airfoil lower surface, now becomes effective for the 
upper surface. The boundary layer of the lower surface therefore becomes lam- 
inar and at times separates. This change becomes noticeable in every polar 
plot as a transition, jump, or point of inflection. 

This phenomenon is especially apparent with bird wing profiles which are 
highly concave on the lower surface of the wing. Here, at a negative angle of 
attack, there is a transition eddy on the wing lower surface extending from the 
leading to the trailing edge and becoming smaller when the model is turned higher 
in + he wind tunnel; then at tangential approach flow to the front part of the profile, 
and thus at relatively great positive angle of attack, it disappears, as shown in 
the polar plot by a c jump at about c = 0. 5. 


III. MODELS AND MEASURING METHODS 

1. The Models Used 

The profiles were chosen to show extreme contrasts in behavior in regard 
to Reynolds number and thus to show the limits of influence of the problem. The 
model wings are rectangular in plan. 

Model dimensions: 

Chord length t = 90 mm 

Span b = 450 mm 

Aspect ratio A = b : t = 5 

The models are not rounded at the tips but are cut off sharply. The 
axis of rotation is at t/4 on the mean line chord. The sting for the rear mount 
is located here on the lower surface to avoid disturbance on the upper surface. 
Trolon plastic was used as a material for the profile model. The rough-milled 
profile outline was cut out by spoke shave according to a glued-on template 
photograph, to an accuracy of 0. 1 to 0.2 mm, and was worked to a good surface 
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smoothness. The flat plate and the curved plate were made of steel sheet 
2. 6 mm thick. Figure 37 shows the dimensions and ratios given the following 
lists 3 and 4, in percentages of chord length, witr the exception of a few angles. 

The camber f of the profile mean line is measured upward from the 
theoretical chord, which serves as the x-axis in profile calculations. At the 
forward end, this theoretical chord goes through the intersection of the profile 
mean line with the point r/2 of the nose circle and to a corresponding point at 
the rear, that is, through the apex if the trailing edge is sharp. 


2. The Characteristic Measuring Method 

The lower limit of the measured range of Reynolds numbers - 
Re « 20,000 - corresponds, with the chosen chord of t = 90 mm, to a stagna- 
tion pressure reading of q = 0. 75 mm water column. The aerodynamic force 
changes to be measured at this low stagnation pressure are of the order of 
magnitude of a tenth of a gram at low angles of attack. Every angle of attack 
setting produces a small change in the tare balancing in experiments with wire- 
hung models; this was reduced, to be sure, by journalling of the model axis in 
two small rollers, but it caused serious dispersion in the subcritical range of 
Reynolds numbers. Accordingly, the conventional polar measurement, in which 
speed is held constant and the angle of attack a in the wind is changed, was used 
only for a check. Instead of this, another measuring method was used success- 
fully. In it, at a fixed angle of attack a the velocity v was changed by stages in 
two test series, first with increasing speed and then with decreasing speed. 

This method will be called "characteristic measurement" to differentiate it from 
the polar measurement. The result of the two characteristic tests series for 
a = 10° on profile N 60 is shown in Figure 39 as an example. The coefficients 

c , c , and c _ „ are plotted as a function of Re. The c line shows with 
a w mO. 25 & 

increasing v, and thus increasing Reynolds number, subcritically from A' to B* 

an increase in lift coefficient and then at B’, at Re = 147,000, a jump of 

c & = 0. 76 to the supercritical line at c & = 1. 06. With increasing v, the c & line 

increases only slightly, but then with decreasing v at point E’, where 

Re = 82,400, it suddenly drops down again to the subcritical line trace A'B\ 

The cause is that from D to E the boundary layer remains turbulent and therefore 
attached at the upper surface of the profile. At E*, the separation point of the 
boundary layer reaches the transition point on the profile; the flow separates 
and remains separated subcritically from A' to B' because the boundary layer 
remains laminar. The points B* and E* correspond on the profile to the boundary 
layer points where the transition and separation points coincide. If the flow has 
attained the super or subcritical state, then the flow remains separated below B 
and attached above E, no matter whether increasing or decreasing v is used in 
measurement. 





Ratios of the Wing Sections Used For Measurement 
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The two lines measured for drag coefficient c^ show contrary behavior 
in comparison with c ; at separation (point E") , the resistance suddenly 

8L 

becomes greater, and therefore the moment, In this manner each profile was 
measured in steps of 2° in setting. With increasing angle of attack, the separa- 
tion and sudden increase displaced toward higher Reynolds numbers, so that the 
available velocity of the test installation no longer sufficed to cause the flow to 
reattach. To permit supercritical measurement of the air forces at these high 
anglos of attack, the flow had to be made to attach artificially by making the 
boundary layer turbulent. To cuase it to attach, a stick was placed horizontally 
In the flow at the height of the model axis for a few seconds. 

The compilation of these characteristic test series foi some 10 to 25 
different angles of attack in which the a span ranged from -25* to +25 
(e.g. , for profile N 60 in Figures 57 through 59) gave a continuous overall 
picture for each profile. The polar plot can be taken from this as a cross 
section for any given Reynolds number of the measured region (compiled for 
N 60 in Table 4) . 



3. The Hysteresis Loop 

The three-dimensional coordinate system c versus c and Re in 

J a w 

Figure 40 makes clear the relationship between polar measurement and 
characteristic measurement. The rectangle CEFB represents a longitudinal 
section through the hysteresis loop, which appears in cross section in file polar 
plot. The measurement of the supercritical polar curve occurs at increasing 
a ( raising up into the wind) , beginning with attached flow at low angles of 
attack, until the flow separatee, o, g. , at 19. 5* for N 60 in Figure 40 from 
point M to G in the polar plot for Re = 147, 000. As measurement continues with 
a decreasing, the flow remains separated from H to B, and then first at 10° 
springs to point C of the polar curve. The following interrelationship exists 
between characteristic measurement and polar measurement: 

1. Measurement with attached flow up to the point of separation: this 
corresponds to 

a) Characteristic measurement with decreasing v ( line DCE) , 

b) polar measurement at high values of a ( line MCG) . 

2. Measurement with nonattached flow up to the point of sudden 
increase: this corresponds to 




Figure 40. Three-dimensional coordinate system c versus 

ft 

c and Re to show the relationship between the characteristic 

W 

measuring method ( longitudinal section) and the polar measur- 
ing method ( cross section) . 


c) characteristic measurement with increasing v ( line AFB) , 

d) polar measurement with decreasing a for separated flow 
(line HB). 

Conclusion: if a wind tunnel measurement begins in a low turbulence 
flow for any reason with an a -setting in the region between f he limit of the 
sudden increase CL and the separation line EG, then - in both the characteristic 
and the polar measurement - the flow after attainment of the desired velocity 
must be ’’attached" by making the boundary layer turbulent, for example, by 
briefly placing a stick upstream. 
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In the negative range of angles of attack, only very thick profiles with 
intensively rounded nose show a hysteresis loop. The test rule must be noted 
as above, and thus the flow must be attached or, for the polar eurve, ot must 
be reduced. 

The lines of sudden increase and separation have their origin at the 
critical Reynolds number Re^ (e.g. , for profile N 60, Re^ is 63,000) . In the 

subcritical region, the transition necessary for attachment of the laminar 
boundary layer flow remains in the turbulent region, so that no hysteresis 
occurs. 


Most wind tunnels are afflicted with an undesirably high inherent turbu- 
lence. With increasing degree of turbulence, Ro , as well as the sudden 

K 

increase and separation belonging to a given angle of attack are displaced to a 
lowor Reynolds number, and therefore show excessively favorable test results. 

In the fully turbulent stream, e.g. , with an upstream turbulence grille 
{ Figure 12 ) , there is no subcritical flow and no hysterosis in the range 
measured here; in addition, instead of the sudden separation at high a -values 
there is a more rounded path of the polar curve. It should be mentioned in 
addition that in the subcritical region and after separation corresponding to F 
to H the turbulence bar is effective only as long as it remains held in front. 
After it is removed the flow separates again both the two cases. 

The flow behavior important for the study of model flight in the sub- 
critical region is found in the characteristic measuring method to be simply 
the prolongation of the lower hysteresis line BF, or the area HBKF, to A. 

For model flight, hysteresis has practical importance for the "stalled" stage 
of flight, for spin conditions, and for the subsequent sudden increase after a 
dive {acceleration at low a) for the purpose of restoration of the normal flying 
position. 


IV. PLOTTING OF THE TEST RESULTS 

The aerodynamic forces were measured on the models at various angles 
of attack lift A ( kg) perpendicular to the direction of airflow, drag W ( kg) 

parallel to the direction of airflow, and moment Mq, 25 ( m • kg) about the axis 

of rotation of the model flying on the theoretical chord 0. 25 t. After conversion 

these are plotted as dimensionless coefficients c , c , and c . as a 

a wm mu, 25 

function of Re. In the Reynolds number pertaining to each test point, the 
atmospheric fluctuations of temperature and air pressure, or the changes or 


air density p and kinematic viscosity v, were taken i nto account in calculation 
of the Reynolds number Re - (t • v/v) - (W2 q/pj/v. The influence of air 
humidity in the test area is negligibly small. 

Equations for calculation of the coefficients in the notation introduced 
by Prandtl are: 


Lift coefficient c « -- ; 

a q • F 


Drag coefficient c = ; 

W q * F 


Moment coefficient c _ __ <= ■ — . ; 

mO, 25 q • F • t ’ 


Stagnation pressure q 


y • v e 


(kg/m 2 ) . 


1. G625 Wing Profile 

The influence of Reynolds number is strongest in round-nosed, thick, 
pronouncedly cambered profiles. As a first example, therefore, the thick 
profile 625 of the Goettingen Aerodynamics Research Laboratory ( AVA) was 
chosen, because this wing section was used in Goettingen as well as by the 
German Aviation Research Institute in Adlershof [19] for characteristic 
measurements, whose results the present measurements supplement. The 
profile has a camber of the profile mean line of f = 0. 06 t at x = t/3 and a 
maximum thickness of d = 0. 20 t at x « 0.30 t ( Figure 38) . 


Plot of c^ Versus Reynolds Number 


The rapid increase in c values at the critical Reynolds number, which 

here is Re - 105,000 ( Figure 41) , is especially striking. In the supercritical 

region, the two most important phenomena "separation" and "sudden increase" 

appear as perpendicular lines: separation with the arrow downward, and 

sudden increase with the arrow upward. The beginnings and ends of the lines 

are delimited by four dot-dashed line segments which have their origin at Re . 

k 

The two upper lines correspond to the supercritical state of "attached flow, " 
and the two lower to the state of "unattached flow. " 




from the subcritical to the supercritical state for profile 625. Arrows pointing 
downward indicate "separation”; arrows pointing upward indicate "sudden 
increase. " A turbulence wire stretched parallel to the leading edge of the airfoil 
displaces the critical Reynolds number from 105,000 to about 20,000. 







The behavior of several a -lines will be explained. The measurement of 

the 4-degree line, for example, starting at Re » 21,000, shows only the low 

c value = 0.2; at Re, there is a stronger increase, and then at Re = 154,000 
a K 

the increase is from c = 0.5 to c = 0.8. Thus it increases to four times the 

a a 

initial value. With attached flow, the further increase of the 4-degree line is 
small. Going in the direction of decreasing Reynolds number, at Re » 108,000 
the flow separates, the lift coefficient drops on the 4 -degree line back to 
•'unattached flow. " For angles of attack higher than 4® , the upper Reynolds 
number limit available in the test installation was not sufficient to cause the 
flow to attach. For example, to measure in the "attached" condition at 12®, a 
stick was held upstream for a few seconds, and making the boundary layer turbu- 
lent in this way caused the flow to attach. At Re = 141,000, the 12® line 
separates and drops back to the 12® line of unattached flow with = 0.5. 

Below 155,000, as long as a turublence grid is held in place or a turbulence 
wire stretched across upstream ( Td) , a premature displacement of the critical 
Reynolds number to about 20,000 results. 

The separation line is identical with the c line up to 14® . At higher 

& lilctX 

angles of attack, the separation line is below the c line. Between the 

ft max 

separated flow - supercritical compared with subcritical - at high angles of 

attack, there was no difference recognizable from the test point series, which 

showed pronounced dispersion here; the c line for 20® was plotted as a horizon- 

a 

tal, although from comparison with other measurements at higher Reynolds 
numbers an increase would be expected here, too. 

In the region of negative lift, the picture shows by the remarkable 

S-shaped curve of the c line an apparent reversal of the phenomena, 

a max 

because here the negative lift coefficient increases with decreasing Reynolds 

number. The value c = 0 results for the -8® line at Re = 100,000 and for 

a 

Re = 21, 000 at a = -2® . The value of c = -0. 38 at Re = 21, 000, calculated 

a max 

from c = 0, is reached with an a setting of -6°. On the other hand, to produce 
a 

an equally large positive lift at Re = 21,000 about twice the a -setting, +11.5®, 
is necessary, that is, the less cambered lower surface is more effective at low 
Reynolds number because the flow here remains "attached" up to -8® , while on 
the strongly cambered upper profile half, the flow has largely separated. The 
setting -2® corresponds roughly to airflow on the direction of the mean line chord 
so that the profile half lying below the mean line chord at Re = 21,000 balances 
the lift produced by the four times more strongly cambered line of the upper side 
for c = 0. In the degree to which at an increasing Reynolds number the point of 

ft 


separation moves toward the trailing edge on the strongly cambered upper sur- 
face, so that the flow produces a higher + component in accordance with the 

stronger camber, the smaller becomes the residual -c value or the stronger 

a 

must the negative angle of attack of the lower surface be for equalization to 
c = 0, e.g. , -8" at Re = 100,000, The influence of Reynolds number on the 

9l 

change in angle of attack at c = 0 is therefore the greater tile greater is the 

a 

difference in the pattern of camber of the two halves of the profile, and thus it 
is least for symmetrical profiles and plates. 

The more that separation phenomena cause the flow picture to depart 
from the potential flow with circulation, and the less does the flow centerline 
coincide with the profile mean line, the more must the study of the separation 
phenomena on the profile halves be drawn upon for explanation of the relation- 
ships. It may be stated here that the flat plate and the symmetrical profile at 
Re = 21,000 and a = 6° reach exactly the same value mentioned above, 
e = 0.38, and this then remains almost at the same level for the flat plate as 

A 

Reynolds number increases ( Figure 47) . For the symmetrical profile it 
increases slowly, while in the picture at hand for profile G 625 on the path of 
the -8° line, -c a slowly decreases. Above -8° the flow separates on the lower 

surface, so that now the slowly increasing stagnation effect of the upper surface 
( previsouly separated in laminar form) causes the phenomena to intersect at 
<x = -18.4°. 


b) Plot of c as a Function of Reynolds Number 
wm 


In the supercritical region, the lift increase achieved by a sudden increase 

or attachment of flow is connected with a sudden reduction of drag, because the 

vortex zone behind the separation point, displaced farther back as it is, has 

become smaller. The phenomena of "sudden increase" ( Figure 42) and 

"separation" ( Figure 43) are delimited by dot-dashed lines. Subcritically the 

-4° line, and supercritically the -6’ line, form the c . value. The 12° line 

a min 

for separated flow runs horizontally, while in the c measurement this did not 
occur until 20°. a 

In the region of negative angle of attack ( Figure 44) the peculiar reversal 
of the phenomena is again evident, with drag first decreasing at small Reynolds 
numbers and large negative angles of attack. With increasing Reynolds number 
and decreasing minus angle of attack, the picture changes from the character of 
the symmetrical profile at Re^ = 105,000 and ^ = -4° to the region of 

influence of the strongly cambered upper surface. 
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Arrow — , measurement with increasing v, gives * sudden change, 
decreasing v gives t separation of flow ( -6° to 12° , profile 625) . 





















Figure 44. Test result: Drag coefficient c wm as a function 

of Reynolds number ( -6° to -20° , profile 625) . Subscript m 
means measured value without a flow correction calculation. 










t 


c) Plot of c . n . as a Function of Reynolds Number 
ml). Jo 


The moment coefficient c . is referred to an axis of rotation which 

m0.25 

(in agreement with the model axis in the measurements) lies 0. 25 t behind the 
leading edge and on the theoretical chord of the profile mean line ( Figure 38) . 
The measured moment M 0 , 25 (m ♦ kg) , the moment coefficient c^^ and the 

formerly customary value c^, whose axis of rotation passes through the lower- 

surface tangent and through the foremost projection point 0 of the profile, are 
connected by the equation: 


M 0. 28 = c _ c 

q • F • t m0.23 m 


0. 25 c + 
n 


yp^26 


(see Figures 37, 45, and 67) . 

If the resultant aerodynamic force or its coefficient 

c = /c 2 + c 2 
r N a w 


is resolved ( Figure 38) 1. at the intersection with the lower-surface tangent - 
normal and tangential to it - then the coefficients of the normal and tangential 
force components are: 


c = c • cos a + c • sin a 
n a w 


c. = c • cos a - c • sin a 
t w a 

2. at the intersection with the theoretical chord - normal and tangential to this - 

then instead of a , the angle a = (a + a) must be substituted: 

6 


c = c ♦ cos (01 + a) + c • sin (a + a) 
ns & w 


c = c • cos (a + cr) - c • sin ( a + a) 

tS W £L 


For conversion of c . to another moment reference axis, such as 
m0.25 

the point of rotation h, there is valid, in view of the sign data of Figure 45: 


x y 

1 S C — — • 0 . JL q 

mh m0.25 t ns t ts 


j 

S; 




•i* 
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Figure 45. 
axis. 


Conversion of c A „ to another moment reference 
mo, 25 


Profile G 625 has the coefficient c _ __ - 0.1 in Figure 46, supercriti- 

mo, 25 

cally with attached flow, from -8° to about +10° . In the measured range of 
Reynolds numbers, this thick profile - in contrast to the other profiles 
( Figure 39) - shows an increase in moment coefficient at the sudden upward 
increase and a reduction at separation from Re = 100,000 to 130, 000. The 
latter is for the range of angles of attack from -8° to +10° . This range of 
Reynolds numbers lies in the area of medium-sized model airplanes. In the 
flying model, the flow in fids range separates in stall form even at low angles 
of attack. The separation, however, does not have the effect of an increase 
of c . measured on other profiles, causing the downward rotation of the 

airfoil nose, desirably for the stabilization process; instead, the tendency to 
turn upward is even increased, by reduction of the moment coefficient. This 
viewpoint is evidence against the use of thick, strongly cambered airfoil 
sections for model flight. The moment coefficient of the -8° line shows the 
highest influence of Reynolds number: At Re = 21,000, c^ 2g is about equal 

to 0 - corresponding to the behavior of the symmetrical profile - and super- 

eritically it is 0. 11. With increasingly negative angle of attack, the change 

in moment coefficient again becomes smaller up to the point -18. 4° , for which 

the position of the axis of rotation at gives a minimum for c „ constant 

mu.zo 

as for the profile insensitive to center of pressure, and insensitive to Reynolds 
number, so that the model when turning free like a weather vane always auto- 
matically takes this position and maintains itself stably in the entire range of 
Reynolds numbers measured here. 
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Mgure 46. Test Results: moment coefficient c 
function of Re for profile 625. 


d) Determination of Profile Drag c 


The polar plots shown in Tables 1 through 5 for the five wing profiles 

were obtained as cross sections from the plots of the coefficients c , c , 

a wm 

an d c mQ 25 versus Reynolds number (for example, Figures 41 through 46) , by 

taking the coefficients at intersections of ordinates Re = 21, 42, 63, 84, 105, 

126, 147, and 168 • 10 s * ) with the a lines. These values are summarized in 

m 

the numerical tables ( Section VIU) nd were converted to two-dimensional flow, 

i.e. , to infinite span (b = «) or to the aspect ratio A = b 2 /F = » according 

to the airfoil theory published by Prandtl in 1918 (see [3], edition I, p. 35, 

Prandtl, ‘'Outline of Airfoil Theory") . For measured c , the conversion 

a 

furnishes the corresponding profile drag c and the pertinent "true" angle of 

W oo 

attack . After about 20 years of publication of airfoil measurements based 
on polar plots for A = 5, now the presentation for A = « is more and more 


*) A series of Reynolds numbers divisible by 70 was chosen to facilitate 
calculation, because for rough calculations the equation Re = 70 • v • t 
suffices ( t in mm, v in m/ sec) . 













becoming the basis of comparison for profile characteristics*) . In addition to 
conversion to A = *> , the influence of finite dimensions of the wind tunnel air- 
stream and the influence of the rectangular plan of the airfoil model must be 
taken into account, as shown in the following sections. 

1 . Influence of finite dimensions of the wind tunnel flow. The pressure 
field of a wing has a definite range of action which is influenced the more by the 
finite diameter of the airstream in the wind tunnel the greater is the wing area 
F in relation to the flow cross section F q and the greater the wingspan b in 

relation to the flow diameter D q . The free jet nozzle of the wind tunnel used in 

the State Engineering School at Cologne is a regular octagon with an inner 

diameter of 700 mm and F =40.6 dm 2 Jet cross section. With a dimension of 

o 

the model airplane of t = 0. 9 dm and b - 4. 5 dm, the ratio becomes: 




airfoil area _ F _ 4,05 _ _1_ . 

jet cross section F 40. 6 10 

0 

The pressure difference producing lift is, as is known, increased in an 
enclosed tunnel by the tunnel wall, and reduced in the free jet. The lift measured 
in the free jet is therefore smaller than the lift effective in open air. In Prandtl'e 
correction, however, the lift is not used but instead the angle of attack a m 

belonging to the measured lift and the drag or c wm * For the free jet with circu- 
lar cross section (see [5], vol. IV, part 2, p. 169, Seyferth, "Studies of Air- 
plane Models in the Wind Tunnel") there results: 


-c 2 • F 
a 

Drag correction Ac = — = — * 8 

& w 8 • F 

o 


57.3° -c F 
a 

Angle of attack correction Aa = — 8 - F * 9 

’ o 


Correction factor 9 « 



+ . • » 


*) used for the first time in 1918 by the author in the research laboratory 
of Professor Junkers for airfoil measurements on the study of the "Influence of 
Thickness Ratio on Profile Properties. " 


a 
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To the octagonal jet cross section F = 40.6 dm 2 there corresponds ft 
diameter of a circle of equal area: 


D = 
o 


7. 19 dm . 


With b/D = 0.626. 8 - 1.032, and the correction then is: 
o 

Ac => c 2 • 0. 01287 
w a 

A a ~ c • 0.738 . 
a 

The drag coefficient corrected for the influence of finite jet cross- 

sectional area, and the pertinent angle of attack which corresponds to the lift 

coefficient c measured for the rectangular airfoil model ( A = 5) are then: 
a 

o=c - c . _ = c -0. 01287 c 2 
w5 wm wk5 wm a 

a = a - of. _ — o' - 0.738 c . 

5 m k5 m a 

2. Influence of aspect ratio A = b 2 F. Pressure equalization at the ^ 
wing tips results in lift reduction, too, so that to achieve the same lift as on a 
wing of infinite span, the true angle of attack must be increased by the 
induced angle of attack g^. In addition, the edge vortices arising at the wing 

tip produce an energy loss which is designated induced drag ( Wi = c^ • q ♦ F) 

For a given lift and given aspect ratio the induced drag is smallest 
when the lift is distributed along the span in the form of a semiellipse, as is 
the case for wings with an elliptical pU n. According to Prandtl, then: 


a, - a - a 
i 00 


57.3° • c 
i 

IT • A 


e = c - c 
wi w w«o 


IT • A 


For example, on the elliptical wing, for the transition from A = 5 to A = •*. 

c « c _ - 0.0637 c 2 ; <* = a g - 3.65 c . 

woo w5 a 49 a 


\ i 


\ 


On the elliptical wing, the conversion of given test values for a wing of aspect 
ratio A to any other aspeot ratio follows the equations: 


9 



3. Influence of roctangular model plan. Lift distribution is more com- 
plete on the rectangular airfoil, and the induced resistance at A - 0 is about 
4 percent greater than on the elliptical airfoil. According to Glauert ( sec [8] , 
p. 132) the conversion equations: 



hold for the rectangular airfoil for the transition to another aspect ratio A . 

x 

The values 5 and t can be taken from the following table: 


3 

4 

5 

6 

7 

8 

9 

10 

11 

0.016 

0.026 

0.037 

0.046 

0.055 

0.064 

0.072 

0. 080 

0.088 

0.097 

0.122 

0.145 

0.163 

0.183 

0.201 

0.216 

0.228 

0.240 


These values agree exactly only for those cases where the slope of the 
lift line dc & / da w - 2ir, corresponding to theory ( cf. Figure 63) . 


Conversion of existing measurements on airfoil models of rectangular 
plan to infinite aspect ratio then gives as a very important result of the study 
of profile drag: 
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and die pertinent "true" angle of attack: 


67.3|° 


a = O'-.*" 

oo 5 


J a 1 + 0.145 


IT 


- of - 4. 18 c 
5 a 


Compared with the elliptical wing, the profile drag is 3. 7 percent nigher, 
and the pertinent angle of attaok is 14.5 percent greater. 

To summarize the three conversion calculations, there result for the 
Cologne measurements: 

c ~ c - 0.07897 c 2 a. ~ a -4. 9? 8c . 
w<w wm a 00 m a 

The numerical tables given in the Appendix contain the values: 

a , c , c i c , « ; c , a and c 

m a r n tiwiA *• • K r k *«• _ r Art « 


nr a wm w5 5 woe 00 


’'mO. 25 * 


In the subcritical as well as in the supercritical separated stats of flow, 
where the elliptical lift distribution is seriously disturbed, the conversion 
equations give uncertain results. Although the two states mentioned here have 
no practical use, the calculation was done to show a uniform comparison picture 
and the relationship of the subcritical with the separated supercritical flow. 

Table 1*) gives a survey of wing profile G 325. Three Goettingen 
measurements at Re = 121, 226, and 430 * 10 s are added as conclusion and 
prolongation to the eight polar plots of the authors' own measurements. 


e) c as a Function of a 
a 00 


The plot of c ft as a function of shows clearly the significant lift 

difference between the subcritical and supercritical attached flow and also the 
relationship between the subcritical and supercritical detached flow, and for 
the latter the "transition" to the sudden increase. Because the conversion of 
induced drag coefficient c ^ and the a change connected with this are based on 

the constant lift coefficient - * c aoe ) » the perpendicular separation and 

sudden increase lines for A = 5 show a slope which is greatest at A — 00 . 


*) Tables 1 to 5 appear at the conclusion of this report. 



The differences now appearing at the beginning and end of the sudden-increase 

lines are proportional to the c differences. These a differences would not 

& 

appear if the equation of airfoil theory could be transformed into a uniform 
angle of attack. 

The picture further shows especially clearly the influence of Reynolds 

number, 1. on the zero-lift /9, which varies from -2° to -15°, and 2. on the 

maximum value of lift, which increases from c » 0. 6 to 1. 5. Up to 

a 

Re = 147,000 or a =7°, the separation occurs at the c peak; then above 

’ oo • ^ a max 

these values, it occurs on the other side of c , after a transition charac- 
’ a max 

teristic of turbulent separation and rounded off. The slope of the lift line 

dc /da agrees closely with the Goettingen lines. The lines from 

ft 

Re = 150, 000 to 430,000 coincide on the straight sections of the lift lines, which 

extend about to the point of sudden increase. In this region, the efficiency rj of 

the profile ( corresponding to the expression dc /a = 2mj) is constant. Subcriti- 

& 

cally, the efficiency is very poor; only in the region of negative lift does it 
become more favorable with decreasing Reynolds number as a result of increas- 
ing influence of the attached flow on the profile underside. At -18°, all lines 
intersect in the neutral point. 


f) The Polar Plot c as a Function of c 
a w« 


Because lift and drag are both changed favorably when the boundary 

layer becomes turbulent, the differences between subcritical and supercritical 

state appear more pronounced in the polar plot. The minimum of the profile 

drag coefficient, which is about 0. 017 at Re = 430, 000 and remains constant 

in a 10° range of angle of attack from c = Oto 0.9, is about double at 

ft 

Re = 105,000 and about five times as large at Re = 21, 000. With decreasing 
Reynolds number, the range of constant drag becomes smaller, so that on the 
elongated semiellipses of similar polar plots the upper and lower curvatures, 
which characterize the beginning of separation of flow, become more 
pronounced. Because c is plotted to five times the scale of c ft , the polar 

plots of the subcritical flow state appear as segments of an ellipse with 
inverted axis ratio. The minimum profile drag for 21,000 occurs in the nega- 
tive c region valid for inverted flight, where - again opposite to the supercriti- 
ft 

cal trend - at small Reynolds number a more favorable (even though still 
nominally poor) best drag-lift ratio c /c , as the tangent to this polar plot 

Woo ft 

shows in the positive c & region, occurs. 
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For most profiles, the line of the moment coefficient referred to the 

t/4 axis shows a gentle S-shape for its basic figure; it starts at negative angles 

of attack, passes through the value c = 0, increases rapidly, and then 

mu. jSu 

becomes almost constant in the entire supercritical range up to separation, 

which generally gives an increase in moment coefficient. Thus, profile 625 in 

the mean range of angle of attack has a value differing only slightly from 

c A = 0.1 and a somewhat decreasing value with increasing angle of attack, 
mu.io 

While separation at Re = 105, 000 reduces the moment coefficient, at 
Re = 168,000 the generally observed increase can be seen. The two separately 
shown subcritical lines for 21,000 and 84,000 show, increasing from zero to 
0. 1 , the S-shape in rough form. 


2. The Flat Plate 

In contrast to the case of the thick strongly curved profile, the influence 
of Reynolds number is least for the thin uncambered plate. The plots of 
Figures 47 and 48 show lines c & and c^ ^ as a function of Reynolds number; the 

lines are almost horizontal over the entire measuring range without a critical 

sudden increase, and show low Reynolds number sensitivity. At Re = 20, 000, 

no subcritical state could yet be seen in the test, because the flow was already 

turbulent over the sharp leading edge. The sharp leading edge acts as a "knife 

edge. " For the result summarized in Table 2, the low sensitivity to Reynolds 

number is indicated by the closeness of the lines. While the measuring range is 

represented by eight polar curves for the thick profile, here two polar curves 

suffice, and for comparison the Goettingen measurement is added again. As the 

lines for c versus a show, the flow separates at a =5.5 and 6. 5° . On the 
a 00 00 

flat plate, the stagnation point movement has the strongest effect, so that the 

leading edge at = 5.5° (« m = 8° in Figure 47) has the effect of a 

"separation edge, " and thus in spite of the steep slope of the lift line at 

Re = 42,000, c is only 0.5, and at ten times the Reynolds number it is 
a max 

0. 7. Separation results in a pronounced increase in drag, as can be seen from 
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Figures 47 ( left) and 48 ( right) . Test results: c and c as 
a function of Re for the flat plate. a wm 

the polar curve; here lift increases further up to a « 25°. In the polar curve, 
the Goettingen measurement shows a higher drag at low angle of attack, but at 
higher angle of attack there is a smaller drag increase because of the greater 
pulling effect at Re = 420, 000. The Goettingen plate has a somewhat smaller 
thickness 6/t and semicircular roundings at the leading and trailing edge, while 
the plate of the present measurements - resembling a thin symmetric profile - 
has an elliptical nose and a sharply pointed trailing edge from % t. 

The lines of the moment coefficients which show a minimum at 

Re = 20,000 in file plot against Re in Figure 49 up to a = 12° , in Table 2 

m 

increase perpendicularly in the plot of c versus c - at c . __ = 0 to 

a m0.25 mQ.25 

<x ao = 4 to 5° . This indicates: up to = 4° the flat plate is fixed in regard to 

center of pressure, and the point of attack of the resultant aerodynamic force is 
at t/4. The result is that the flat plate is longitudinally stable in flight, when 
the center of gravity lies at t/4. Lanchester in about 1890 experimentally 
demonstrated the inherent stability of the flat plate with his "flying board. " 

The process of stabilization about the transverse axis is further promoted by the 
following. The rapid increase of c n _ to higher plus values above 5° indicates 

an increase in the nose-down moment. The increase in moment begins about 
1 to 2° before achievement of the separation angle of attack, as the graph shows. 
If the flying model, is, for example, threatened with stalling in a gust, then 2° 
before the separation, the restoring moment takes effect. By trimming of the 
center of gravity in the range of 0. 25 to 0. 30 t, the angle of attack at which the 
flat plate should fly can be varied to achieve the best drag-lift ratio or, what is 
more important to achieve the minimum sinking velocity for the model airplane, 
the best ceiling factor reciprocal. 
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Figure 49. Test result: c as a function of Re for the 

flat plate. m0 ‘ 25 


For the understanding of longitudinal stability, that is, the stability about 
the transverse axis, the following is important. Departing from the normal 
flight position with decreasing a (= nose down) , velocity v becomes greater, 
and at increasing a (= nose-up) it becomes smaller. Thus a definite speed 
belongs to an angle of attack desired for the neutral flight position ( the resultant 
aerodynamic force passes through the center of gravity) . How should the 
stabilization about the transverse axis of the airplane model be affected: 

Exactly as for any airplane ( Figure 50) : 

1. The plane shoudl fly neutrally in the principal flying position at 
normal speed, and thus without turning moment. 

2. At increasing speed or decreasing angle of attack, the airplane 
should be tail-heavy. 

3. With decreasing speed or increasing angle of attack, it should be 
nose-heavy; in brief, it should always automatically return to the 
neutral flight position. 
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Figure 50. The properly stabilized model airplane is nose-heavy 
at excessively low speed, and tail-heavy at excessively high speed, 
so that it always automatically arrives at the proper neutral posi- 
tion at the right speed. 


The moment about the center of gravity S (expressed by the moment 

coefficient c ) of a flat plate (A = 5) is plotted in Figure 51 for various center 
m 

of gravity rearward positions x/t as a function of a . Along the line c m “ 0 

there is neutrality. Above it, there is nose heaviness, and below it, tail heavi- 
ness. If the center of gravity is too far forward, that is, ahead of the point of 
application of the aerodynamic force, say at 0. 2 t, then the airplane noses over 
and accelerates and the angle of attack is reduced until, at zero position, that 
is, in the vertical diving position, the neutral position is reached. If the center 
of gravity is at 0. 25 t, then neutrality occurs at 5° . If a disturbance causes 
higher angle of attack, or if the model is launched with a higher angle of attack 
than 5 s , then it becomes nose-heavy as desired. The range useful for the model 
is shown hatched in the graph; at a low angle of attack or at high speed it should 
preferably become tail-heavy, and then have a neutral position at that angle of 
attack where the optimum ceiling factor reciprocal prevails, as shown by the 
polar plot. Finally, it should become nose-heavy. The center of gravity loca- 
tion between 0. 25 and 0. 3 t meets this requirement. Practical flight tests gave 
the best flying characteristics at 0. 28 t. Stabilization occurs as an oscillation 
about the neutral position at 7® . If the center of gravity is at 0.35 t, however, 
then the change from pronounced tail heaviness to pronounced nose heaviness 
occurs over a small increase in the angle of attack. The oscillation leads to 
"festoon flight" (the phugoid flight curve) . The restoring moment has been made 
excessively large, so that it makes the model dynamically unstable. At even 
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«f=0.25 

*-=0.30 

-■=0.35 


--=0.40 


Figure 51. Explanation of the longitudinal stability of the "flying 

board. " Moment coefficient c referred to the center of gravity 

m 

as a function of angle of attack a. 


farther rearward positions of the center of gravity - as exemplified by a falling 
postcard - the phugoid flight becomes zig-zag flight, and the postcard oscillates 
downward or spins about the transverse axis. *) 

The following may be said in summary about flow over die flat plate: 

The phenomena are controlled by the transition vortex occurring over 
the sharp leading edge and producing a high suction peak at that point. The fol- 
lowing phenomena are explained by the rapid increase of the suction peak above 
the leading edge: 


*) For a detailed treatment of these phenomena see: Schmitz, F. W. , 
"Stability Experiments on Model Airplanes" (Luftfiahrt und Schule, vol. n. 






1. The high increase in lift per degree of angle of attack, dc /da. 

a 

2. The advance of the pressure point to t/4 ( on the cambered profile 
at low angles of attack it moves toward the leading edge to about 
t/ 2 and at large angles of attack to t/3) . 

3. The transition vortex is effective even at the smallest angle of 
attack, such as from a - 0 to 1°, as can be concluded from the 
lift increase and from the center of pressure advancement to t/4. 

4. The bridge building extending to the trailing edge as the transition 
vortex broadens gives, prior to separation at a further increase in 
lift, a periodic sinking of the suction peaks and arching of the pres- 
sure distribution (which, dampened by the scales, is probably 
indicated as an average value) as an explanation of the early origin 
of the center of pressure movement to the trailing edge and the 
simultaneous large increase in drag. 

When bridge-building is illustrated by a profile in liquid, it shows 
clearly how, by increasing the camber of the plate by filling out the transition 
vortex space into a sickle profile or, better, by cambering the plate especially 
at the leading edge, the separation following the bridge-building can be 
displaced. 


3. Cambered Plate 417 a 

Otto Lilienthal showed the superiority of the cambered plate compared 
with the flat plate by model tests, as indicated in his book Der Vogelflug als 
Grundlage der Fliegekunst ( Bird Plight as a Basis of Aviation) (see [14], 
p. 74) . Nevertheless, the result of the present measurement in the range of 
small Reynolds numbers is surprising and interesting. 

First of all, the lines in Figures 52 through 55, by their straight course 
without a sudden critical increase, again show the character of low Reynolds- 
number sensitivity of the model, expressed in Table 3 by the small difference 
between the lines for Re = 42, 000 and the lines of ten times that Reynolds 
number of the Goettingen measurement. Accordingly, two polar curves at 
Re = 42, 000 and 168 , 000 suffice to represent the test series here. At 
Re = 42,000 and 0 = 8°, the cambered plate at a steep angle of attack reaches 
a high c & max = 1. 06 and thereby is superior to the flat plate by a factor of 

two, and superior to the thick profile 625 by a factor of 2. 6 at the same angle 
of attack calculated from c =0, because subcritical flow still prevails on the 
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latter. This difference is compared in Figure 1 (converted to A - 5) . The 
slope of the lift line for the cambered plate surpasses the lift Increase for the 
flat plate, up to now designated the optimum. An equally favorable picture 
results at Re = 42,000 from the behavior of the profile drag in the polar plot. 
The minimum is c^ = 0.027 at c & = 0.7, which corresponds to an optimum 

"profile drag-lift ratio of c Wof / c a - V^e* Th Q position of this optimum drag- 

lift ratio at as high as possible c gives good soaring ability; high drag at 

ft 

c 0 means a low diving speed. These two properties give an extraordinarily 
ft 

favorable situation for the model airplane. On the other hand, the pronounced 
mean-line camber f/t = 0. 06 results in an undesirable center of pressure 
movement, as shown by the moment curve of c & versus c^ For longitudi- 
nal stabilization, large rear control surfaces are required, while the flat plate 
flies with inherent stability ( with the proper center of gravity location) . 

For construction reasons, however, a certain profile thickness is 
required for a wing section on medium-size models. The measurement given 
in Table 3 for Re = 420, 000 gives an idea of how a thickening of the thin, cam- 
bered plate 417 a works compared with the bird-wind profile 417 at the same 
mean line and a thickness d = 0.063 t. The polar curve shows, as do most 
bird-wing profiles, a sudden increase in c^, because the resistance here 

decreases above c & = 0.5. The slope of the lift line is less up to this turning 

point than for 417 a, but above it the two are equal. The influence of the thick- 
ness ratio on the moment curve can also clearly be seen: the moment coeffi- 
cient c „ is smaller throughout and is nearly constant over a wide c range. 
mu,65 ft 

The superiority of the thin, cambered airfoil in the Reynolds number 
range under 100,000 has three causes: 

1 . Favorable interaction of tangential approach flow to the leading 
edge at large angles of attack and nose turbulence, 

2. The large camber of the lower surface, 

3. The relatively small curvature of the rearward upper surface. 

In supplementation to the basic laws described in Section II, several 
essential viewpoints for profile study result from the explanation of the causes. 
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Figure 52. Test result: c as a function of Re for the cambered 
plate 417 a. 
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Figure 53 ( left) and 54 ( right) . Test results: c as a 
funotion of Re for the cambered plate 417 a. Arri 


Concerning 1. If a tangent is drawn with the theoretical chord at the 
leading edge intersection of the profile mean line, the front tangent angle </> = 22° 
is found for cambered plate 417 a ( see Table 3 and List 4) . This means that at 
22® angle of attack the leading edge is parallel to the main flow. Tangential 
approach flow on the leading edge occurs even at low angle of attack, however, 
because of the convergence of flow before the leading edge ( Figure 56a) , since 
even at a w = 0° a large lift is present (c « 0.4 at Re = 42,000; Table 3) , 

fit 

which along with the wing also draws the flow up high ahead of the leading edge. 
With plate 417 a, the flow separates at ~ 8° , and the Rat plate at 5.5°. If 
this 5.5° range of the flat plate (from tangential approach flow to separation) is 
used as a comparison, then the leading edge of the cambered plate would 
„ experience tangential approach Row at about a M = 8 - 5. 5° = 2. 5° . At this 
angle of attack, the lift coefficient has a relatively high value c - 0.7, and the 

fit 

drag coefficient c reaches its minimum of 0.027. This gives a profile drag- 

lift ratio c /c of 1:26 at Re - 42,000; converted to the aspect ratio A = b: t 
w a 

= 10 it becomes. c ^/c & = 0.042/0.7 = 1:16.5. Corresponding to Figure 1, at 

A = 5, it is c V / C SL ~ 1*12. This drag-lift ratio, unusually favorable for a model 

wing, is little changed even at Re = 20,000, and it surpasses any other profile. 
This is the explanation for the astonishing flying performances of the indoor air- 
plane models with a curved thin film skin stretched over them. 
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Figure 55. Test result: c as a function of Re for the 
cambered plate 417 a. m ' * 


The stagnation point movement here, too, along with the nose turbulence, 

controls the phenomena, which can be seen more clearly on the lines of profile 

417 in Table 3 than on the lines of the cambered plate 417 a. Beginning at high 

negative angle of attack, the flow separates on the lower surface. If the negative 

angle of attack is smaller, then at -8.2° (or -5°, for 417 a) the lower-surface 

flow attaches in part behind a transition vortex extending to the leading edge. 

At -8. 2° , the center of pressure is at t/4, because then o n =0 (suction 

m0,25 

peak on the lower surface, pressure peak on the upper surface) . The transition 
vortex on the lower surface now collapses, so that the resistance becomes 
smaller and the pressure distribution more complete, and thus the moment 
becomes larger until at -2. 7° the c mQ g5 line has reached approximately the 

constant value of 0.11, because the flow now adheres along the entire lower 
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Figure 66. Circularly cambered (are] plate (a) and sickle 
profile ( b) . 


surface. Between and -0. 8°, the transition vortex remainder disappears 
under the leading edge, so that now the stagnation point can move around the 
leading edge at this point. Now the upper -surface flow can first become turbu- 
lent, and at « 0.24 it attaches at the trailing edge, as the turning point of 
the line c^ versus at 0,24° shows. From 1" to 5°, the profile drag io nearly 

constant. At about 5° , the transition point has reached the leading edge. Then 
the broadening of the transition vortex begins which leads at 9° to separation. 

The three characteristic intermediate points of profile 417 are not 
recognizable on the cambered plate 417 a. The difference is probably caused 
by the profile thickness d = 0. 063 t and by the sharp leading edge of the profile, 
with r = 0. 007 t against d = 0. 029 t and r = 0. 0145 t for the plate. The three 
measurements of the plate all In all agree well with one another. 

Concerning 2 and 3. An important conclusion for profile design of the 
model airplane is given by study of the profile lower surface. If a circularly 
cambered (arc) plate with radius r ( Figure 56a) has flow approach at a = 0% 
then in contrast to the flat plate a large lift is produced. The explanation of the 
lift as a force directed transversely to the principal direction of flow can be 
visualized from a mechanical pressure standpoint as a centrifugal effect. In a 
homogeneous parallel flow no force directed transverse to the flow is effective, 
because every curvature of a flow produces a force component directed radially 
outward. The particle of specific mass p moving on the curve path experiences 
a centrifugal effect Z = ( p • vVr) , which manifests itself on the circularly 
cambered plate as a positive pressure on the lower surface and a 3 a negative 
pressure on the upper surface, and in both cases is greater the more pronounced 
is the curvature or the smaller the r/t value, and thus from a theoretical 
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potential standpoint at r/t — 1/2. If tfte concavity of the lower-surface is filled, 
then the sickle profile used for ship's propellers ( Figure 56b) occurs, where 
the straight lower surface obviously reduces the centrifugal effect because now 
the flatter profile mean line serves as mean line of curvature of the flow. 
Although at increasing angle of attack of the circularly cambered plate the flow on 
the lower surface now improves, separation begins earlier on the downstream 
half of the upper surface the more intensively it is cambered and the smaller 
the Reynolds number. The effort to camber the lower surface more and the 
upper surface less then leads to the thin, cambered plate as a favorable form to 
give high lift with low drag at small Reynolds numbers. 

Forward camber, that is, the displacement of the highest camber, which 

is, for example, at x -= 0.5 t on the circularly cambered plate, to x = 0.38 t 

for plate 417 a, by enlarging the leading side tangent angle and reducing the 

trailing side curvature gives higher c values and displacement of the 

q. max 

region of the best drag-lift ratio to high c values, c = 0. 8 to 1 here. 

a a 


4. Wing Profile N 60 and N60R 

For these profiles there are comparison measurements made in 
America in the NACA variable-density wind tunnel at Re » 3 • 10 6 ( NACA 
Technical Notes, No. 388) . Profile N 60 has at 0.4 t a mean-line camber 
f = 0. 04 t and a thickness d = 0. 124 t. Profile N 60R is derived from N 60; 
from 0. 3 t on, *he mean line is given an S curvature according to the equation 

y = 0.3 x(l - x) • - x) 

and the raised trailing edge part of the profile, acting as a tail plane, makes the 
profile proof against center -of-pressure change. *) In addition to making clear 
the action of the S curvature, the measurements will show the influence of 
Reynolds number and e specially to what extent resistance to center-of-pressure 
change still applies at low Reynolds numbers, making profile N 60R suitable for 
flying-wing models. 


* ) These profiles were chosen because A. Lippisch, in Profile Collection 
No. 4 and 5 of the periodical Flugsport ( Sport Plying) ( Frankfurt a. M. ) , using 
these airfoil sections as an example, has given a readily understandable calcula- 
tion of the lift and moment coefficients c and c on the basis of the 

ao mo 

Birnbaum-Glauert equation; this is widely available because of the broad circula- 
tion of Flugsport. 
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a) Profile N 60 


The critical Reynolds number is at Re = 63,000, from Figure 57. Up 

to G! = 16. 5° , the envelope curve of the separation points is identical with the 
m 

c line, whose steep slope ends at 16.5° or c = 1.3. This range from 

a max a 

Re, = 63,000 up to the end of the steep c rise at Re = 110,000 is similar 
k a max 

to the "transition" in flow on the friction plate ( Figure 6) from the laminar 

state at Re, to the fully turbulent boundary layer flow. From Re = 110,000 

K 

there follows then c of the 17° line, while separation occurs at larger a , 
a max 

but smaller c . The origin of a test series for profile N 60 for a - 10° has 
a m 

already been given in Figure 39 as an example. In the plot of c versus Re 

wnx 

( Figure 58) , the c lines for supercritically detached flow, which result as the 
w 

horizontal prolongation of the subcritical lines, are omitted between 10° and 19°. 

In the graph of the moment coefficient ( Figure 59) is a gathering point for 

Re, = 63,000 at c „ „„ = 0.093, at which the two dot-dashed envelope curves 
k m0.25 

of separation points have their origin. In contrast to profile G 625, separation 
here produces the desired increase in c ^ that is, increase in nose-over 
moment. 

Eight polar curve -j of the measurement from Re = 21,000 to 168,000, 
and for comparison the NACA measurement Re^ » 8 million are plotted in 

Table 4. The difference between the polar curves for 21,000 and 168,000 is 

considerably smaller than in the polar curves of the same Reynolds numbers for 

the thick profile 625 in Table 1. The drag coefficients are about V s to V 2 

smaller, but the "c range of approximately constant profile drag" is consider- 
St 

ably higher. The slope of the lift line is steeper, and even in the lower part of 
the line for Re = 21,000 it is almost parallel to the 8 million line. The c 

st max 

values are higher than for the 625, and the influence of Reynolds number on the 
change of the zero-lift angle is only half as great. The lines for "unattached 
flow" intersect at a point at at^ = 16.5° or in the polar curves at c = 0.275. 

These intersection points of the influence of Reynolds number are to be observed 
at larfe positive and large negative angles of attack in all measurements of 
Reynolds number influence. 

The measurement gave a more uniform curve for moment coefficient 
than for the 625, so that the transition from the subcritical to the supercritical 
flow and the supercritical mean value c mQ ^ = 0. 08 can clearly be seen. 


93 













0.14 1 
C n»0.25 



140 160,000 


Figure 59. Test result: c A as a function of Re for profile 
N 60. m0 * 25 

b) Profile N 60R 

The critical Reynolds number is Re = 73, 000 ( Figure 60) , and thus 

somewhat higher than for N 60, because the upper surface from 0.3 t on is 

somewhat more sharply curved. From Re. on, the slope of the c line is 

ic & niEX 

less steep, and it is still not ended at Re = 168, 000. Thee point of N 60R at 

a max 

168,000 is about 15 percent lower in lift coefficient, in spite of the 50 percent 
higher Reynolds number, because through the S curvature the camber rise of 
the mean line is reduced from f = 0. 04 t to 0. 03 t; the stabilizing property of 
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the profile is thus purchased at the cost of lift. Conversely, the values for 
a given are smaller ( Figure 61) ; in particular, however, the c^ values 
have become considerably smaller ( Figure 62) . 


In the measurement of Re .. = 8,000,000, the point c =1.4 for 

eff a max 

N 60R ( Table 5) is about 12. 5 percent lower than the 1. 6 for N 60 ( Table 4) 

while the slope of the lift line dc /da remains the same. The plot of c 

a 00 a 

versus a ^ shows a remarkable path, with bending and crossing of the lines. 

The measurement was repeated several times to verify the correctness. Super - 
critically, the Reynolds number increase of 84,000 does not cause a parallel dis- 
placement of the c lines, as in the cambered profile, but instead, as for the flat 
a 

plate, there is an angular rotation about the intersection point at - 1°, so 
that the line Re = 168,000 has the same slope up to c & = 0. 8 as the line of the 

NACA measurement. At negative angles of attack the picture changes suddenly, 
because now the influence of the profile lower surface predominates; seen from 
the mean line chord out, it appears as the profile half of a symmetrical profile 
with a sharp point on the trailing edge. As in profile 625, there are the effects 
which in profile 625 disturbed the course of the lower polar curve half because 
of the alternating play of separations; here, however, the polar curves in the 
-c region show a smooth curve very similar to the polar curves of the flat plate 

3 

( Table 2) . The span of uniform profile drag extends from c = 0 to about 0. 7. 

3 

The c . value is of course not at c = 0.4, as it is for profile N 60, but, for 
w min a 

example, at Re = 84, 000 it is at c - 0. 1, so that the model has a large and 
undesirable diving speed. 


The lines for the moment coefficient, which for clarity are shown 
separately in Table 5 for subcritlcal and supercritical flow, give the most 
significant results. The line of the NACA measurement follows the ordinate 

0. 8. The profile shows little change in center of pres- 


c „ n = 0 up to c = 
m0.25 ^ a 


sure up to oi 


S 00 


= V 


at Re _ = 8,000,000. Above c =0.8, an increase in 
eff a 


'm0.25 


occurs and quickly increases in the neighborhood of the separation point. 


The effect of the S-shape here is the same as that for the flat plate: at 
R e e ff = 8,000,000 a rectangular airfoil with this profile would, like the flying 

board, be longitudinally stable at the correct center of gravity location of 
x = 0. 25 to 0.3 t. As can be seen from the graph, however, with decreasing 
Reynolds number there is an increasing nose -heaviness moment. In addition, 
it is unfavorable for the process of stabilization that the moment coefficients 
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Figure 62. Test result: c as a function of Re for profile 
N60R. m0,J8 


become smaller with increasing angle of attack, becoming larger only after 
separation. For this profile to be used on a flying-wing model, accordingly, a 
sharp sweepback and twist would be additionally necessary to achieve longitudinal 
stability. 


V. COMPARISON OF THE TEST RESULTS 

1. The Lift Slope dc a /da* 

For the flat plate, the theoretical slope of the lift line is the familiar: 


dc /da = 2ir 
a 00 


( a in degrees of arc) 






or 


dc^/do! 0 ^ = 27r/57.3 = 0,11 . ( or in radians) 

This value is used as a comparison value for lift efficiency r/ of a profile 
form through the equation: 


dc /d a = 27T»7 . 
a 00 

According to Figure 03, at Re n 420,000 the line of the Goettingen 
measurement for the flat plate coincides exactly with the slope 2 tt; it corresponds 
to the profile efficiency tj ~ 1. Frequently the lift slope 27r is called the theoreti- 
cally achievable optimum; actually, however, the cambered plate always exceeds 
the value 2ir, with a surprisingly steep slope of 2 . 777T at Re - 42,000, and in the 

o region of -0. 1 to +0. 6 (cf. Table 3) , which corresponds to a profile efficiency 
a 

i) - 1.38. 


The three profiles 625, N 60, and N 60R, which are very poor producers 

of lift subcritically, have the same lift slope of 1. 67ir at Re = 168,000 in the 

range of c = 0.3 to 0. 8; at Re .. = 8,000,000, the two profiles N 60 and N 60R 
a ett 

achieve the maximum value of 1. 8ir up to c = 0. 9, and thus 77 = 0. 9. Only thin 

& 

profiles similar to birds' wings, such as 417, approach the cambered plate. A 
disadvantage is found along with the advantage here: the steeper the lift slope, 
the greater is the sensitivity to gusts, because small changes in the direction of 
approach of air cause large fluctuations in lift. The thin highly cambered pro- 
file therefore requires a larger tail plane for balance than does a thick profile 
possessing the same camber. Lilienthal was the victim of insufficient gust 
stabilization of his cambered airfoil. 

It should further be pointed out that for all profiles and plates in the 
supercritical range the lift slope can be only slightly improved with increasing 
Reynolds number, as shown by study of the plots of versus in Tables 1 to 5. 

On the other hand, the increase in the c range or a range, in which dc /da is 

3 3 - 

constant, is also greatly dependent on the Reynolds number. 
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2. The Ze.o Lift Angle ft 

In Tables 1 to 5 every lift line of the plot of c^ versus intersects the 

li» e c = 0 at a different angle of attack. Those zero lift angles ( designated 
0 - 

by ft in Figure 37) differ from one another to a Uogree depending on how sensitive 
to Reynolds number the profilo is. In Figure 64 the zero lift angles of the five 
profiles are plotted as a function of Reynolds number; the angles are not 
expressed in terms of the profile chord but of the theoretical chord and therefore 
are designated as fl a . Tho flat plate has the lowest Reynolds number influence, 

D 

and the thick profile 63ii has tho greatest, Suberitioally /9 g & 0" for all. At 

Re , a sudden increase occurs and tho camber difforonoo grows larger between 
K 

the upper and lowar surface. Aftor an ineronso of ft in tho transition region, 

b 

there is generally a decrease in zero-lift angle above tho transition. 
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Figure 64. Change in the zero-lift angle of attack /3 g of the mean 
line chord as a function of Re. 










The zero-lift angle plays a role in theoretical lift calculations. The true 
or aerodynamic angle of attack in terms of the direction of approach of air- 
flow at c & = 0 ( in Figure 37) was used in Figure 63 for a comparison of lift 
slopes. 


3. c a max Slope As a Function of Re 

The c lines of the measurements are shown in Figure 66. The 
3. max 

steep c & max slope in the transition from the subcritical to the fully turbulent 

condition occurs for N 60 from Re^ = 63,0 ft 0 to Re = 110,000, and for 625 

from Re^ = 105,000 to Re = 157,000 as can be seen from the deviation of the 

dot-dashed separation line from the c line. According to Figure 60, there 

a max 

is no separation yet at the upper range of measurements for N 60R, so that the 

transition region extends further here. On the flat plate, no subcritical condition 

could be found, but its c value is too small in comparison with the cambered 

a max 

plate, which is superior to all others up to Re = 100,000. At this Reynolds 
number, the superiority of the profiles begins. At Re = 420, 000 and 8, 000, 000, 
according to Figure 63, the cambered plate is far inferior to the profiles. The 
compulsion of the premature supercritical state is at the expense of the c 

a max 

value at high Reynolds numbers; this is a general rule for other measures of 

artificial turbulence. The gradual c increase found here at all Reynolds 

a max 

numbers measured is based on the pulling effect of the outside flow. 

4. The Plot of the Minimum Profile Drag c 
As a Function of Reynolds Number w«oimn 

The line of minimum drag for the flat plate in Figure 66 corresponds 
approximately to the laminar friction coefficient c^ on the friction plate of 

Figure 6, which reaches the transition at Re, = 500,000. Thee , line of 

k w*> min 

the cambered plate is almost horizontal. The thick airfoil section 625 is again 
subject to the strongest influence of Reynolds number. The "transition" from 
the subcritical to the supercritical state occurs here below Re, , in contrast to 

K 

the c slope, where the transition lies above Re, . With the thick airfoil 
a max k 

section 625, the transition ranges from Re = 95,000 to Re. = 105,000, with 

K 

c becoming % smaller and the angle of attack for c , chf jged from -4 to 
w Woo min 
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Figure 65. Comparison of the five profiles with respect to change 
of their supercritical maximum lift value c in the Reynolds 
number range of 20, 000 to 168, 000. a maX 


-6° . For the two intermediate airfoil sections, the influence of Reynolds num- 
ber on their profile form is correspondingly less. The c /c ratio 
determines the "velocity lange" of the airplane. a ma * w mIn 


5. Movement of the Center of Pressure e/t 

The change of moment consists of the simultaneous change of aerodynamic 
force and its point of application or the change of the lever arm, as a result of 
movement of the center of pressure. Although the moment coefficient suffices 
for calculation of longitudinal stability, still the representation of the movement 


105 







•if— 





20 40 60 80 100 120 140 


160,000 


Figure 66. Comparison of the five profiles with respect to change 
of profile minimum drag c , in the Reynolds number range of 
20,000 to 168,000. Weo mln 


of center of pressure or the center-of-pressure distance e/t from the leading 
edge, gives clearer comparative pictures. According to Figure 67 there can be 
set up: 


C nS ‘ (e _ 0,25 


C m0.25 


t . 



Figure 67. Sketch for determining the center-of-pressure 
distance e. 
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Accordingly, the center-of-pressure distance Is: 


+ o.25 . 

‘ c nS 

The center-of-pressure movements resulting for five profiles as a 
function of angle of attack and Reynolds number are shown in Figures 68 
through 72, for airfoils of aspect ratio A « 5 . Because the reference axis for 
°m0 25 * e ^ * 0,1 theoretical chord, here c^ and a should be referred to 

this chord, and thus, since A - 5, c ng g and oj^ should be used in calculation. 

For wing section N 60, the center of pressure at Re = 168,000 and -4° 
is at e/t = 1, that is, at the trailing edge; at 0° it is at 0.43 t, at 12° at 0.3 t, 
and then at 19° it suddenly increases to 0.41. Below -4°, the line moves toward 
plus infinity and at -8° returns from minus infinity and then at -12° reaches 
0. 25 t. With increasing Reynolds number, the lines generally show an equidis- 
tant displacement to the right, and the movement becomes greater; subcritically, 
the movement is smaller, but die center of pressure is farther back. The 
NACA measurement shows the greatest movement, with a center of pressure 
point moving from 0.5 to 0.31 to from 0 to 20°. 

On profile N 60R, on the other hand, at.Re Qff = 8,000,000, the center 

of pressure is fixed at t/4 from 0 to 10°; then at 22° it moves back favorably to 
0. 3 t. ( Before the flow threatens to separate in a stalled condition, a stalled 
condition, a restoring moment is effective which moves the wing back into the 
normal flying position. ) As Reynolds number decreases, the range of file fixed 
center-of-pressure position rapidly becomes smaller; at Re = 84,000, e reaches 
the t/4 position at only one point, 12° , and then suddenly increases to t/3 as 
flow separates, hi comparison with N 60, the movement for low Reynolds num- 
bers is even greater; there can no longer be center-of-pressure point fixity 
spoken of here. 

For the thick profile 625 ( Figure 70) , the course of the line for 
Re = 168,000 is almost the same as with N 60, but the difference, as far as the 
line for Re = 21,000, 1 s considerably greater, however, and the jurve in the 
region of large negative angle of attack changes suddenly. As expected, the 
flat plate ( Figure 71) , which has a fixed center of pressure from -6° to +6® , 
shows the most uniform behavior. After a further increase in angle of attack of 
3° , the center of pressure moves from t/4 to t/3, and after another 7° , to 
0. 4 t. In the measurement with the plate, it is seen that the line for 
Re = 42,000 even reaches a center-of-pressure position e/t = 0.23 at a = 3°. 
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Figures 68 ( above) and 69 ( below) . Center-of-pressure movement e/t on air- 
foils of A = 5 with profiles N 60 and N 60R, as a function of angle of attack at 
various Reynolds numbers. o»g is the angle of attack of the theoretical chord of 

the profile mean line. NACA measurement, Re = 3 million or Re „ = 8 million. 
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Figure 70. Center -of-pres sure movement on the thick Goettingen 

profile 625. 
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Figure 71. Center- of-pres sure movement on the flat plate 
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For the cambered plate 417 a ( Figure 72) , the lines show little deviation among 
themselves; nevertheless, the movement of the center of pressure is large 
corresponding to the high camber of f = 0. 06 t: the line coincides with the lino 
of profile 625, which has the same camber. 



Figure 72. Center-of-pressure movement for the cambered 
plate 417 a. 

6. Profile Drag*Uft Ratio c n a and Profile Ceiling Factor 
Reciprocal c we ^c a 3 as a Function of Re 

The plot of profile drag-lift ratio c /c and profile ceiling factor 

Woo a 

reciprocal c 2 /c 3 as functions of Reynolds number furnishes the most 
woo a 

important comparison. These profile evaluation coefficients are shown in 

Figure 73 for file five Reynolds numbers 42,000, 84,000, 168,000, and 420,000 

or up to 8 • 10 6 shown in Figure 73. To avoid fractions, the reciprocals are 

plotted, namely, the lift-drag ratio c / c as abscissa and the ceiling factor 

c 3 /c 2 as ordinate. a Wo ° 

a w oo 




The result: at Re = 42, 000, the cambered plate 417 a predominates 
over all others (camber f = 0.06 t; thickness d = 0.029 t) ; at 84,000 it is still 
superior to profile N 60. On the other hand, at Re = 168,000, profile N 60 
( f = 0. 04 t; d = 0. 12 t) has surpassed every other one and from then on main- 
tains the superiority. The thick profile 625 (d = 0.20 t) is comparatively poor. 

In spite of this high c value as a result of the large camber f - 0. 06 t 

a max 

(cf. Figure 65) , it is inferior because of its high profile drag. At Re = 168,000, 
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profile 625 is equal to the cambered plate 417 a, but at higher Reynolds number, 
for example, 420,000, it is again inferior. Because of its relatively great 
thickness and its low camber, the S-shaped profile N 60R is comparatively 
unfavorable in the Reynolds number range of model flying, although it is 
favorable at higher Reynolds numbers. *) 


7. Comparison of Measurements at Very Low and Very High Velocities 

If flow at an air velocity v approaching the speed of sound (c w 430 m/sec) 
passes a profile, then the profile properties become much poorer, as is well- 
known, and this occurs first for thick profiles {see (8), vol. IV, part 1, p. 402; 
A. Busemann, Gasdynamlk (Gas Dynamics)}. It therefore occurs to compare 
the effect of this deterioration with the present deterioration measured for pro- 
files in the transition to the subcritical condition. Although the causes are 
basically different, nevertheless some similarities may be pointed out. In 
Figures 74 and 76, the measurement for profile G 625, converted to the same 
aspect ratio, is shown compared with NACA measurement for a propeller pro- 
file of the same thickness (d/t 0.2) . Just as a sudden increase to the sub- 
critical flow occurs from Re / 105,000 to 84,000 in Figure 74, in Figure 75 a 
similar increase to a ''hypercritical" state can be seen frcwtp v/c = 0.6 to 0.6. #*) 
Its polar curve at v/c = 0. 85 coincides with the polar curve for Re = 21,000 of 
subcritical flow, t) Although in Figure 74, from Re = 105,000 to about 150, 000, 


*) For the practical case with the airplane or the model airplane, where 

the total resistance coefficient c = c + c + c .the differences in the 

Wg W°o Wl ws 

ceiling factor and drag-lift coefficient are considerably smaller, and thus the 
influence of profile differences is not so clearly evident. Here c represents 

WS 

the proportion of the sum of harmful drags of fuselage, tail planes, etc. , per 

unit of wing area F, at dynamic pressure q=l;c =Ze • f/F. 

ws w 


**) This flight speed v = 0,6 c * 204 m/sec corresponds approximately 
to the absolute speed record of 1939 of the Messerschmidt airplane which 
reached 210 m/sec. 

t ) More recent measurements have shown that the NACA values found in 
the closed tunnel do not apply until velocities about 50 percent higher; this 
quantitatively shifts the comparative picture. 
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Figures 74 ( left) and 75 ( right) . Comparison of measurements 
at very low and very high velocity. 


separation does not occur for incompletely turbulent boundary layer flow without 
a transition suddenly over the entire wing, the upper polar curve branch at 
Re = 430,000 shows the polar curve's characteristic rounded-off behavior for 
a fully turbulent flow before separation, the same as the polar curves in 
Figure 75. v 

The Reynolds number influence upon approach to very high velocity as 
well as to the lowest velocity is less the thinner the profile. Thus the same 
requirement holds for both: the higher the velocity above v/c = 0. 5, and the 
smaller the Reynolds number below Re = 100, 000, the thinner must be the 
profile and the sharper must be the leading edge (see [2], p. 215) , because 
thick profiles are unsuitable as a result of the very high profile drag. 
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VI. PRACTICAL CONCLUSIONS 


1. Choice of Profile for the Model Airplane 

A good range and a long duration of flight are required for the model air- 
plane; that is, it should glide as far as possible and remain up as long as 
possible and thus fly with as low a sinking velocity as possible. This assumes a 
choice of a wing section which gives the maximum attainable values for profile 
drag-lift ratio and ceiling factor reciprocal for the Reynolds number at which the 
model airplane flies. All in all, for profile choice, only the representation 
shown in Figure 73 is valid. These profile evaluation numbers depend upon tho 
profile parameters f/t, d/t, their location x/t, and in addition, on the aose 
radius r/t. Determination of these profile parameters is the practical goal of 
this study. The effort to choose profile thickness as large as permitted by 
aerodynamic evaluation, to achieve maximum spar depth, leads according to 
Figure 73 to a limiting valuo of profile thickness for the model airplane of 
d s 0.12 t, with a camber f = 0.05 t, at Re = 200,000. If, in addition, the 
parameters of the cambered plate 417 a are chosen as optimum for Re = 50, 000, 
then with these guiding values in Figure 76 there results a system for profile 
selection or for the incorporation of existing profiles in the suitable range of 
Reynolds numbers. In this system, some Goettingen profiles have been plotted 
for choice (see [3] editions I, III, IV; profile dimensions in m, p. 27, IV, 
p. 63) . Camber and thickness for these are generally at :: = 0. 3 t. More 
recent profile studies [13] showed as the most favorable location of maximum 
camber 0.15 t, for achievement of a high c value. For the profiles chosen 

here, the maximum camber actually lies farther back, but the camber of the 
mean line at first increases steeply, for example for profiles 375, 381, 417, 
and others, with a high tangent angle at the wing leading edge, so that about the 
same effect is achieved. 

In the comparison of these profiles, it becomes clear that they are more 
like the sections of birds' wings than the wing sections of modern airplanes. The 
measurements corroborate on one hand the correctness of the model tests of 
Lilienthal who, starting from thin, curved bird wings, determined that these and 
the cambered plate are the most favorable profile forms. Since he did not know 
the influence of Reynolds number, however, he assumed erroneously that this 
result was the most favorable form for large-scale designs, too. In 1910, 

Junkers showed for the first time the advantages of the thick profile for the air- 
plane of that time and thereby showed aviation the way to the cantilever monoplane 
Up until about 1918, however, the biplane and triplane, with thin wings, ruled 
the field as the most favorable from the construction standpoint, although they 
were aerodynamically unfavorable. On the other hand, modern model airplane 
building often makes the error of copying the wing section of the modern air- 
plane as accurately as possible for the model airplane, using wing sections with 
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Arith convex lower surface and sections with too small a mean- 
example, 2 percent) or using sections from sailplane examples, 
but too thick, American characteristic measurements ( for 
tribution by Seyferth, p. 205) have given impetus to the use of 
is, too; in a plot of c versus Reynolds numbers, at the lower 

a 

,000 the thick heavily cambered profiles USA 35 A ( d = 20 

ttingen 387 (d = 15 percent) reached the highest c values 

a max 

it the same time the lines for decreasing Reynolds number 

ther increase in c increase to be imagined. At 
ci max 

owever, the plot of the minimum profile drag shows the most 
s for these thick profiles. It is therefore to be recommended 
value be used to obtain a lower profile drag, because 

IllltA 

shows that the profile drag-lift ratio and ceiling factor 
ire decisive for effectiveness. 

in Figure 76 can serve only for a preliminary estimate 
d on only two profiles, N 60 and 417 a. Accordingly, further 
cambered plates and profiles in this range are necessary, 
jrs f/ 1 and maximum-camber locations x/t, to make more 







The profiles shown in Figure 76 all have 6 percent camber, and thus 
high moment, which must be compensated by a suitable tail plane or by sweep- 
back. For calculation of moment compensation by a tail plane, the existing 
measurements of the flat plate can be used. 

The "rigging angle of incidence" of the airfoil, that is, the angle at 
which the profile chord is attached to the longitudinal axis of the fuselage, plays 
an important role. The best drag-lift ratio of a wing occurs at a very definite 

angle of attack, which, after conversion of the c and a values to a chosen 

w» 00 

aspect ratio (and after addition of the resistance of the fuselage, tail planes, 
etc. , in terms of the wing area) , is given by the pola* line tangent to the polar 
curve in the polar plot. The optimum ceiling factor reciprocal lies at a some- 
what larger angle of attack. On the Thermik model airplane, designed for 
optimum climb, the rigging angle of incidence should be chosen higher than for 
the model airplane intended for the best distance performance. The optimum 
value of the ceiling factor reciprocal is found in the polar plot (see Figure 1) 

by plot of a c line as a function of c 3 /c 2 . 

a a' w 

An interesting possibility of profile shaping to obtain high spar depth is 
given by the wing section of the hawk ( Figure 77b) which is known to be an 
excellent soarer (see [16], p. 212) . The hawk profile has on the lower under 
surface and behind a steeply cambered thin nose, a step which is also found in 
the stork, crane, and flamingo. Behind the step lies the upper-arm humerus 
or the forearm bone and the thick extensor. The step is recessed farthest at 
file elbow. At the airfoil nose are the flexors, especially thin for the land 
soarers. Upper arm, forearm and flexor form a triangle which is covered by 
the leading-edge wing skin. In birds with beating wings (ornithopters) , the 
flexor is more strongly developed and more thickly feathered. According to 
the observation of Hankin, birds with a section called the eagle profile in 
Figure 77a are superior in stroking flight; on the other hand, birds with the 
hawk profile ( Figure 77b) are faster and can climb at a higher rate in soaring 
flight. Hankin found that the step is most strongly developed in the soaring 
birds having the best efficiency and the highest wing loading, such as the kite. 

He states that the kites, in spite of higher wing loading, always attack the 
eagles in play and are faster in turning and climb better than the eagles in 
horizontal paths. He further observed on birds with the hawk profile a tendency 
during soaring to buck (gust sensitivity) and to create noise, which could no 
longer be heard in steep gliding flight. Idrac (see [15], p. 41) was able to 
measure on the kite the lowest sinking velocity v =0.42 m/sec, for a horizon- 

y 

tal velocity v^ = 7 m/sec, corresponding to a drag-lift ratio of: 


116 

t 










mmemt* 


Figure 77. a) Eagle profile, b) hawk profile. 
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These values, measured in an upwind and converted to gliding flight in still air, 
probably cannot be achieved by model airplanes in still air. Only flight 
measurement in calm - for the flying model airplane, too - gives a satisfactory 
basis for calculation, because only here are all wind fluctuations with respect 
to time and place eliminated. If the model airplane has demonstrated good 
performance in a calm, it will give an optimum performance in a head wind, too, 
and then the wing loading can be adapted to the wind strength, by additional 
weight, say. Instead of calm, it is better to say '’nonmoving air," because below 
a certain height thermals can be effective locally without a noticeable horizontal 
wind movement on the ground. 


2. The Wing Plan for the Model Airplane 

Is it advisable, following the example of the large airplane, to use the 
elliptical or trapezoidal plan for the model airplane? The elliptical plan, as is 
well known, gives an elliptical lift distribution and thus the lowest induced drag 
for airfoils of a given aspect ratio b 2 / F, according to the airfoil theory of 
Prandtl. The trapezoidal wing with arcs on the end represents an approach to 
the elliptical wing and is simpler to make from a construction standpoint. The 
following disadvantage counters these advantages. If the wing chord at the tips 
is reduced for example to t/2, for a trapezoidal airfoil, then the airfoil tips fly 
at half the Reynolds number. On foe large airplane this influence is small, but 
on the model airplane it is of greater importance the smaller is the supercritical 
Reynolds number span in which the airplane model flies, because the airfoil tips 
may find themselves in the subcritical separated state even in the normal flying 
condition ( Figure 78) . 
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Figure 78. At the wing tip of the trapezoidal wing, the Reynolds 
number corresponding to the narrowing to t/2 is only half as 
large as at the root of the wing. In the rectangular wing, the 
Reynolds number is constant along the span, making this plan 
more favorable for the small model airplane. 


The result is that the airplane with the tapered plan works correspond- 
ingly unfavorably at the wing tip because of the low Reynolds number. For the 
model airplane, the rectangular plan is therefore the most favorable form. 
Hawks and eagles have an almost rectangular wing plan when soaring. 

Extremely tapered designs have been attempted in various models of 
high-performance sailplanes; these are to be considered unfavorable because 
the wing tips at low Reynolds numbers worsen the drag-lift ratio more than it 
is improved because of the intended reduction in induced drag. To maintain 
supercritical flow at these wing tips and to reduce the danger of separation in 
the stalled flight state, sometimes the use of wires stretched ahead of the air- 
foil nose at the wing tips is of advantage ( cf. Subsection VI 4c) . 


3. The Model Airplane Propeller 

Practical conclusions can be drawn from the test results for the propeller 
of the model airplane, too. What is true for the airfoil holds for the propeller. 

Of course, here centrifugal forces act on the flow, since particles of the boundary 
layer are directed outwardly so that the flow on the model airplane propeller is 
improved and thus the influence of Reynolds number is reduced here (see (17], 
p. 33) ; this influence will be ignored here. The Reynolds number of the pro- 
peller is least at the hub and increases up to the blade tip. As a consequence of 
the test results, the blade cross section would have to have an especially thin 
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and pointed profile at the hub, but for strength reasons it is generally made 
thicker; nor is the blade. plan of the propeller of the large airplane a suitable 
prototype for the model airplane propeller. As shown in Figure 79, instead of 
the blade width decreasing outwardly, a uniform blade width - or perhaps one 
with a semicircular shaped blade tip - is to be recommended. The table fan 
with its thin, curved sheet blades and blade width increasing outwardly thus is 
a better example for the model airplane propeller than is the large airplane 
propeller. The part of the propeller lying near the hub works least efficiently. 
On the model airplane with a thick Juselage, it is therefore advisable to cover 
the hub with a hood as in the large airplane since this, as a turning fuselage 
point, improves approach flow. With respect to blade profile and hub covering, 
this also holds for the propeller of the small wind tunnel fan in instruction in 
aeronautic physics. 




Figure 79. In contrast to the airplane propeller of the large 
airplane ( upper picture) , on the model airplane as in the 
table fan, constant or outwardly increasing blade depth is 
more favorable for achieving high Reynolds number on the 
most effective part of the propeller lying at about 0. 7 r. 


It must further be required of the model airplane propeller as an airfoil 
that it work supercritically. Near the hub, of course, a subcritical region 
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must be accepted. Because the thin cambered plate is supercritical even at 
Reynolds numbers of 20,000, supercritical work is possible with the proper 
choice of pitch and speed for this blade profile even at relatively small diameters, 
as will be shown in an example. 


A propeller has the following dimensions: diameter d = 200 mm, chord 
length t = d/8 = 25 mm, speed n = 1300 rpm. The model airplane has a flight 
velocity of v = 5.4 m/sec. The radius at which Reynolds number 20,000 pre- 
vails is to be determined, under the assumption that 20,000 corresponds to the 
critical Reynolds number of the thin blade profile used. 

Velocity v in the propeller circle is somewhat greater than the air 
velocity v w , because the propeller must accelerate the air to produce a thrust, 
as the test on the stand, where v^ = 0, makes very clear. This velooity dif- 
ference is called the slip of the propeller. With a slip = (v - v )/(v) = 0.1 
= 10 percent, v becomes 6 m/sec. If u designates the peripheral speed of the 

blade element at radius x, then the speed of this blade element on the propeller 
path is: 



The Reynolds number of this blade element would be 

w ♦ t /u 2 + v 2 • t 

Re — = 20,000 . 

Because u = 2 x 7r n/60 ( m/ sec) then: 



It is customary in aeronautics to refer the initial design of a propeller 
to a mean radius of 0.7 r, because here is found the highest energy transfer on 
the propeller blade. It appears analogously useful here to give the Reynolds l 
number of the blade element at 0. 7 r as the mean Reynolds number of the entire ? 
propeller. Using | 

I 

u 0 .7 = 2 • 0. 7 r • it • n/60 = 0.073 • r • n | 

i 
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then 


Re 0 . T a n/( 0.073 rn) 2 +. v 2 • ^?*1 . 

v 

The 200-mm diameter propeller calculated in the example here would 
therefore have a Reynolds number Re 0 7 = 20,000, so that the connection with 
profile measurement can easily be seen. F. Gutsche (see [7], p. 54) gives 
Re fc = 150,000 as the lower critical limit for the usual profiles. With the 

values in the example above, an Re of 150,000 would be achieved only at a 
radius of x = 0. 63 m, and thus would be unattainable in the model airplane pro- 
peller. The transition from laminar to the turbulent boundary layer flow for 
the usual profiles of the measurements above would occur more favorably, in 
the range of 50,000 to 100, 000. If, for example, profile N 60 were used uni- 
formly along r, then according to Figure 57 the transition would begin at 
63 , 000 and would be finished at 113 , 000. This would have the following practi- 
cal results: 


1. The part of this propeller operating at Re 63,000 would be suberiti- 
cal, that is, would have high resistance and low thrust. 

2. Flow adheres to the blade element operating at Re = 63,000, but 
it separates at an approach flow angle of 4°. 

3. The blade element operating at Re = 113,000 works with a fully 
turbulent boundary layer and does not experience separation until 
c a = 1*3, or at 16.5°. For N 60 the Reynolds number of 113,000 

would be die lower limit of the supercritical condition, which would 
be required at least 0. 7 r for a propeller. 

4. The c & max increase according to Figure 57 is still small on the 

propeller section up to the blade tip; the separation angle of attack 
increases further to 19° at Re = 150,000, however. 

For small model airplanes, the gain achieved by a large Reynolds 
number through the rectangular plan of the propeller blade ( Figure 79) and the 
wing ( Figure 78) is larger than the reduction in induced drag resulting from 
tapering of the plan. 
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4. Measures for Artificial Turbulence 

With turbulence-free approach flow assumed, the boundary layer is 
made turbulent ; 

a) by wall friction, 

b) by knife edges, 

c) by turbulence grilles and turbulence-inducing wire, 

d) by sound waves. 

Convergence reduces turbulence of the flow and divergence increases it. 

Concerning a) , on the plane, completely smooth, friction plate, the 
boundary layer first becomes fully turbulent at Re « 10 7 , according to 
Figure 6; the rough wall becomes turbulent at lower Reynolds numbers, so that 
roughness along the wing nose can sometimes be favorable for the model air- 
plane because the artificially increased turbulence attaches the flow sooner than 
on the smooth model wing. 

As is well known, the rough wing surface on the large airplane has 
exactly the opposite effect ( Figure 80) . On the rough wing, the maximum lift 
value drops to half at high Reynolds numbers. At Re = 118,000, the influence 
of roughness is zero according to the chart; it has an improving effect at 
Re < 100,000. Figure 81 shows especially how a rough wing nose behaves at 
higher Reynolds numbers. Measurements on wings with rough surface are 
currently being made; they will show whether the gain achieved at Re < 100, 000 
is not countered by a worsening in the drag-lift ratio as a result of the increased 
supercritical resistance. 

It should be mentioned at this point that birds' feathers are far removed 
from the smoothness of a painted surface. Their relative roughness is not 
much less than the waves in the sheet-metal skin of a Ju 52, with the difference 
that the ribbing of the bird's feather runs in every direction and transverse to 
the direction of flight. Dragonfly wings even have transverse folds whose depth 
is about 0. 1 1; in addition, the leading edge is sharply notched to make the 
boundary layer turbulent. 

Figure 11 shows how the value of Re^ for a sphere can be reduced by 
roughness. 

In respect to b) , the effect of pointing the wing leading edge has already 
been discussed in detail in connection with stagnation point movement and the 
transition vortex formation. 
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Figure 80. Influence of surface roughness on c at various 

a max 

F.eynolds numbers. ( From Technical Note No. 457; grain size 

« 0. 12 mm; relative roughness k/t = 8 . 10“ 4 ; NACA variable- 

density wind tunnel. ) 3 • 6 • 10® corresponds to the Reynolds 

number of a small sport plane, 1. 78 * 10® to that of a sailplane. 

118,000 corresponds to the Reynolds number of a model airplane 

of medium size; the influence of roughness on c here eauals 

a max H 

zero. At even lower Reynolds numbers, the rough wing has a 

relatively small but higher c value than the smooth wing. 

£i max 


In respect to c) , all effective measures on the model airplane for 
artificial turbulence are analogies of the classic sphere experiment of Prandtl 
and Wieselsberger ( Figures 10 and 13) , in which the boundary layer is 
turbulent by a wire loop; in turbulence-free approach flow, the critical transi- 
tion is reduced thereby from Re = 405,000 to a fraction of this Reynolds number. 
Through the mounting of a turbulence grille ahead of the wing ( Figure 12) , it is 
very simple to make the airstream of a wind tunnel completely turbulent to 
allow study of the effect of turbulence by comparison of measurements with and 
without the grille. As an example, Figure 82 shows a Goettingen measurement 
(see [5], vol. IV, part 2; contribution of Seyferth, p. 189) with and without 
grille for profile 387 (d/t = 0.14) at Re ■ 84, 000. Without a turbulence grille, 

the maximum lift value is c = 0. 8, and with a grille 1.4, nearly twice as 

max 


!3 








Figure 81. Influence of local roughness on the polar curves. 

{ From Wiese lsberger, Ergebnisse der AVA (Goettingen 
Aeronautical Research Institute Reports [3]) ; Re = 630,000.} 
For the large airplane, a rough wing nose, for example, a 
small amount of ice, can halve lift through premature separa- 
tion as a result of large turbulent friction loss in the thickened 
boundary layer. 


large, while the profile drag is % to V g smaller. This attractive result makes 
it seem desirable to achieve a similar resul* by measures for artificial turbu- 
lence in the flying model airplane, too. The tests of course have not been 
finished, but a few essential results have clearly been shown ( Figure 83) . As 
already mentioned, it is sufficient if the boundary layer on the wing upper sur- 
face is turbulent, so that instead of the grille a single wire or twine stretched 
along the wing nose suffices. The result achieved is just as favorable as with 
a turbulence grille. According to a report by Lippisch, model airplanes with 
wing profile 387 showed a peculiar festoon flight path whose explanation is now 
simple. The model airplane flew with a Reynolds number only slightly above 
the critical Reynolds number Re fc = 80,000, so that in horizontal flight, for 

example, according to the chart here, the flow separated at -0.3°. The flying 
model airplane '’pancaked" and became faster, so that the Reynolds number and 
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Figure 82. Polar ourves of an airfoil taken with and without 
turbulence grille. ( From R. Seyferth [5], Goettingen profile 
387, 6/t = 0.14, t = 60 mm, b * 600 mm, v = 20 m/sec, 

Re = 84,000; material, unpolished wood. ) 


thus the lift became greater; the model airplane again pointed upward and the 
process began all over again. When a new model airplane of this type was built, 
it was possible to reproduce die process. When a turbulence wire was stretched 
ahead of the wing, however, the same model airplane showed very good flight 
performance even in still air. For example, on a 1-km long gliding flight it had 
the excellent drag-lift ratio of 1:12. A single twine piece was used as a pre- 
stretched wire. 

There is an Italian patent according to which prestretched wires or tubes 
delay the separation of flow in the stall condition. The patent claim does not 
cover the model airplane in whose Reynolds number range the wire is effective 
but applies only to large airplanes where it increases ( ?) the separation 
tendency. 
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Figure 8a. Polar curves taken with and without turbulence wire. 
Author’s measurements? t = 90 mm, b = 450 mm, v = 13 m/sec, 
He = 82 , 000; material , Trolon, shaped with profiled steel razor 
blades. 


To permit better understanding of the effect of the turbulence wire, 
profile 625, which has been pronounced too thick for model flying, was tested 
with a prestretched wire ( Figure 41) , The result showed that the turbulence 
wire displaces the critical Reynolds number from 105,000 to Re » 20,000, so 
that thick profiles become advantageous, too, because they give new construction 
possibilities. 
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•ast to wires arranged perpendicular to the flow direction, wires 
w direction reduce turbulence, as stated in a patent by M. Kramer. 

t and most interesting finding is the measurement on profile 
>) . This profile is derived from wing section N 60, with an 
j for the mean line from 0.3 t on. The effect of this is that the 
the if 4 axis equals zero, and the profile has a fixed center of 
a high angle of attack, as shown by NACA measurements at 
000. It is therefore suitable for the construction of tailless 

planes. This airplane form has its particular advantages for 
> building, too, but the results with this form have been Insignifi- 
aasurements described here make clear. The lift values are coa- 
ler than on profile N 60, but It is the moment curves which show 
•tant results, because unfortunately the moment becomes greater 
? Reynolds numbers. The smaller the model airplane, the more 
le measures additionally necessary to insure flight stability. 

i as the prestretched wire makes the boundary layer turbulent 
owever, It makes the rear part of the wing profile effective to a 
cizontal stabilizer. The moment to be equalized does not com- 
ar, but it is only half as large. The following four advantages 
usly obtained from the prestretched wire: 

3 flight state is supercritical at Re = 20,000, and thus even for 
smallest model airplane. 

38 sweepbaok and wing twist suffice to stabilize die flying-wing 
del airplane. 

3 relatively thick profile makes possible great construction 
ength and the use of a rubber motor. 

3 prestretched wire in addition insures effective protection in 
lision against a tree or a wall edge. 

tice shows, with a strong up current any symmetrically con- 
. with any profile flies. The turbulence wire, however, furnishes 
simple means, easily applied to the model airplane, to make a 
3 able to fly even in a weak up current or in calm. 


owered model airplane, the propeller fulfills the same purpose 
ce wire. It is thus understandable why American tests show 
s, with propellers of large diameter, about % of the wingspan, 
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Figure 84. By the turbulence wire (or twine) stretched parallel 

to the airfoil leading edge, the moment in this Reynolds range is 

cut about in half for profile N 6 OR, ( The c values are some- 

a max 

what low here with the turbulence wire: the wire is mounted some- 
what too high. ) 


to have especially favorable flight characteristics. In this sense, the use of a 
one-blade propeller makes possible, in addition to a large propeller circle 
diameter, the advantage of a high (propeller) Reynolds number because of high 
blade chord length. In the canard, the forward stabilizing surfaces create 
turbulence for the center airfoil part. 

In bird flight, too, various manifestations of artificial turbulence creation 
can be seen. For example, the soaring flight in column of the hawk is to be 
explained by the turublence created by the hawk soaring ahead for the one behind. 
In the V-formation flight of birds of passage, in addition to the utilization of the 
ascending induced vortex and the airwave of the bird flying ahead in ornithopter 
flight, the turbulence in the wing boundary layer of the succeeding bird plays a 
role, too. Utilization of the two latter effects is also true for the file flight of 
ducka. In addition, the swarm-like flight of birds and insects can be explained 
as turbulence effect - in addition to the use of a local thermal. 


5. Notes for Instruction in Aeronautical Physics 

To show the effect of measures of artificial turbulence by means of 
turbulence grilles, turbulence wire, and so forth, as well as the subcritical and 
supercritical flow around spheres and thick, round-nosed airfoil models, the 
airflow of the wind tunnel fan must be made turbulence-free to some extent. A 
simple means for making turbulence visible is the flame probe ( Figure 85) . In 
the turbulence-free stream, the gas flame is completely smooth, Hlta the ihm<nm> 
flow of a slightly opened water faucet ( Figure 86) , but in the turbulent flow it is 
disordered and undulant, especially in the peripheral zone, where outside air is 
entrained. The author's own experiments made a fan of conventional design 
(by a firm selling aids to education (casing * 0.5 m diameter, nozzle 0.3 m 
diameter) so turbulence-free by the incorporation of a wire screen g ( Figure 87) 
and a stilling stretch h that with a round-nosed thick airfoil model (profile 625, 
t = 50 mm, b = 100 mm) the placing of a turbulence grille or wire in front, 
gave two to three times the lift value, depending on the angle of attack, as shown 
by the two-component balance ( Figure 88) . 

The following is worth noting for the fan. Inlet edge a and inlet hub b 
improve the approach flow to the propeller c, and the rear part of the stream- 
lined body e improves the outflow. The spin of the airstream caused by the 
propeller rotation is reduced 

1. by a curved guide plate d, 

2. by the honeycomb f. 
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LAMINAR TURBULENT 
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Figure 85. Flame probe for mak- 
ing visible the turbulence state of 
an airstream by a gas flame. 



Figure 86. Water flow fe.m a 
water pipe, a) laminar at low v or 
Re, b) turbulent at high v or Re, 
c) laminar flow made turbulent by 
a nail. 


The propeller causes large globules of turbulence; these are; 

1. converted into smaller globules by the honeycomb, and then 

2. converted into fine -structure vorticity by the wire screen g, and 

3. attenuate in the stilling stretch h, and finally 

4. are eliminated as far as possible by constriction in the nozzle. 



Figure 87 . Wind tunnel for school experiments to produce as 

turbulence-free an airflow as possible. 


When several screens are used, it is advantageous to place them a certain dis- 
tance apart and arrange them so that the air flows through the finest screen last. 
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As the experiments showed, the screen here was effective only in connection 
with the stilling stretch. As the screen, insect screen of as fine wire as possible 
was used. The stilling run was made of 1.5-mm thick plywood sheets. For the 
turbulence grille, a simple wooden frame with twine stretched over it at about 
10-mm spacing sufficed. 


6. Order of School Tests 

a) Laminar and Turbulent Flow 

First Experiment. Water flow in a water pipe, laminar at low Reynolds 
number and turbulent at higher Reynolds number ( Figure 86) . 

Second Experiment . Turbulence of the airstream of a fan, like the one 
in Figure 87, made visible by a flame probe ( Figure 85) . 

a) Laminar gas flame in a turbulence-free airstream. 

b) After the placing of a turbulence grille upstream ( Figure 12) , the 
airstream and thereby the flame become turbulent. 

Third Experiment. Stepwise elimination of a fan’s turbulence ( Figure 87) 
shown by means of a flame probe. 

a) Fan without honeycomb and nozzle. 

b) With honeycomb, 

c) With honeycomb and nozzle, 

d) With honeycomb, screen, and nozzle, 

e) Without screen, but with honeycomb, stilling run, and nozzle, 

f) With honeycomb, screen, stilling run, and nozzle. 

5 

t 

b) The Sphere in Low-Turbulence and Turbulent Airs tream 

( The sphere must be fastened to a tube at the rear. No support wires or 
the like can be permitted to make the air turbulent on the upstream side of the 
sphere.) 


Fourth Experiment. Measurement of the drag of the sphere with the 
drag balance: 

a) In subcritical flow in the low-turbulence air stream, 

b) At supercritical flow by the placing of a turbulence grille upstream. 

This experiment shows only the "sensitivity to turbulence" of the sphere. 
To aid in understanding of "sensitivity to Reynolds numbers" of rounded bodies, 
therefore, it must be emphasized that at sufficiently high speed or high Reynolds 
number, as a result of the greater "entraining effect of the external flow" 

(" stronger impulse exchange) the boundary layer becomes turbulent or super- 
critical by itself, at Re^ = 406, 000 in laminar flow and correspondingly earlier 

in turbulent airflow, and the critical sphere Reynolds number of a wind tunnel 
therefore serves as an indication of turbulence for the tunnel. The turbulence 

factor is TF = 405, 000/Re . The following calculations are to be recommended 
here: * 


1. The Reynolds number of the fourth experiment, after the velocity is 
calculated from the dynamic pressure. 

2. The velocity which at the given diameter of the sphere used in the 

fourth experiment would be necessary to attain Re = 405 000 

k ’ 

Fifth Experiment. Measurement of sphere drag without and with a rubber 
ring around the sphere equator in the low-turbulence airstream to demonstrate 
that only the turbulence of the boundary layer and not that of the total airstream 
is important ( Figure 13) . 

Sixth Experiment. A pressure measurement to demonstrate the effect 
of the transition to supercritical flow on the pressure state at the downstream 
side of the sphere ( Figure 10) . The measuring tube should be led along the 
rear sting, with the opening within a few millimeters of the sphere. 

a) Subcritical: negative pressure ( - sign) . 

b) Supercritical: positive pressure (+ sign) with a turbulence grille 
located upstream. 

Reference may be made here to determination of the critical sphere 

characteristic Re, : 
k 
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1. At ohange of the rear pressure from - to +, which means the value 
of p = 0 ( Figure 14) . 

2. In drag measurement at the value •- 0.3 ( Figure 11) . 

Seventh Experiment. Compulsion of the supercritical transition on a 
sphere with rough surface: 

a) By coarse sand glued to the leading side of the sphere, 

b) By a long-haired fur ring on the roar side of the sphere. 

Eighth Experiment. Visualization of the subcritlcal, separated flow and 
the supercritical attached flow by threads glued ahead of the sphere equator 
( or by smoke; Figure 13) . 

Ninth Experiment. Sharp-edged bodies ( flat circular plate and so forth) 
have the same resistance in laminar as in turbulent airflow, because the 
separation point is fixed. Sharp-edged bodies are insensitive to turbulence and 
thus are not sensitive to Reynolds number. This relationship can be demon- 
strated if the drag coefficient is determined for circular plates of various 
diameters and at various velocities, that is, at various Reynolds numbers; in 
all cases the coefficient will be c = 1. 1 ( Figure 9) . 


c) Airfoil Model in the Low- Turbulence and Turbulent 
Airstream 

The discussion of Figure 9 permits an instructive transition from the 
sphere to the airfoil. Sensitivity to turbulence and to Reynolds number are 
shown there in connection with the critical Reynolds nui.ijer and its dependence 
on the fineness ratio 6/ 1. This figure permits the following basic conclusions 
to be drawn: 

1. Progressing from the sphere to the slender airfoil profile, Re fc 

becomes smaller and enters the range of flight of the bird and the 
model airplane. There it has a considerable influence on ability to 
fly, while the flow state for a large airplane which, for example, 
has 10 to 500 times as large a Reynolds number accordingly is 
always supercritical. 
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2. 


The smaller the Reynolds number the thinner must the profile of 
the airplane model or the flying animal ( bird, butterfly, or insect) 
be to permit it to fly through a turbulent boundary layer with super- 
critical attached flow. 

Tenth Experiment. Lift and drag are measured with a two-component 
balance on an airfoil model with an aspect ratio A ~ b : t = >2 : 1 and a round- 
nosed, thick wing profile, for example, 625: 

a) In a low-turbulence airstream, 

b) In a turbulent airstream with a turbulence grille mounted upstream, 

c) With only a turbulent boundary layer on the upper surface, with a 
turbulence wire mounted ahead of it. 

The three measurements are done at a constant angle of attack of about 
6° and at constant velocity and then at another angle of attack and velocity. For 
an aspect ratio smaller than 2, the two tip vortices keep the stream attached, so 
that the subcritical region does not occur. 

Eleventh Experiment. Hardly any lift change can be seen on a thin airfoil 
model of the same aspect ratio no matter whether it is with or without a turbu- 
lence grille. As a result of nose turbulence of the sharp leading edge, the 
condition is supercritical even without a turbulence grille. 

Twelfth Experiment. To verify the effect of nose turbulence, we use in 
this experiment an airfoil model of the same thickness as profile 625 but with a 
sharp leading edge. 

Thirteenth Experiment. We measure two polar curves for the airfoil 
model with profile 625, from 0 to 30* a and then from 0 to -30° a at a constant 
dynamic pressure: 

a) In the low-turbulence airstream, 

b) With a turbulence grille. 

We plot the values a, A, and W' in a table and measure the resistance of the 

sting; after the subtraction of the latter, W results. We know the dynamic 

pressure q and we calculate F and plot for each measurement c and c on 

ft w 

millimeter graph paper, c at five times the c scale ( Figure 88) . In addition, 

W ft 

the Reynolds number should be calculated after determining v. 
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Fourteenth Experiment. We attempt to determine the most favorable 
position of the turbulence wire at various angles of attack for profile 625. 


Fifteenth Experiment. We construct two similar model sailplanes with 
profile 625, t » 150 mm, b = 1200 mm: 


a) With a prestretched wire (cf. Figures 83 and 91) 

b) Without prestretched wire, 

and start both at the same time. Differences in construction should be elimi- 
nated by alternate application of the wire. 

Sixteenth Experiment. To obtain data for calculation of the model flight, 
an attempt must always be made to measure the flight in still air ( Figure 89) : 
a stop watch, a measuring tape, and a protractor for gliding angle, somewhat 
like that shown in Figure 90, are the most essential aids. It is necessary to 
determine the gliding angle, the gliding velocity, the sinking velocity, the 
horizontal velocity, the wing loading, c , c , the ceiling factor reciprocal, the 

Reynolds number, and all these at various rigging angles of incidence, to find 
the best values . The author has given instructions , and simple calculation 
examples, in the magazine Luftfahrt und Schule. * ) vol. VII, No. 4. 


STARTING AND 
STABILIZING RUN 


UNACCELERATED GLIDING 
PLIGHT IN NEUTRAL 
.POSITION 



PULLOUT STRETCH 
NEAR THE EARTH. 


■^—8=20 m— ' L 







Figure 89. Flight measurement of unaccelerated gliding flight in 
still air by means of stop watch, protractor, and measuring tape. 


2 . 


* ) Verlag C. J. E. Volckmann Nachf. E. Wette, Berlin-Charlottenburg 



Figure 90. Protractor for determining the gliding angle. On a 
short flight stretch, the flight path is sighted on from a point 
located to the side. The protractor and the flight path are seen 
as in Figure 89, One pupil does the sighting and a second makes 
sure that the plumb bob hangs free and still and reads tne angle 
on command of the first pupil. In contrast to this, on a longer 
flight stretch the protractor is set at the starting point and the 
model airplane is sighted on through the crosshairs by pupil 
No. 1 from the instant of start until it lands. Pupil No. 2 reads 
the angle, for example, every 5 seconds. Every 5 seconds, 
pupil No. 3 gives the command to read. Pupil No. 4 notes in a 
previously prepared table the angle. Pupil No. 5 observes the 
model airplane through a telescope and upon landing gives the 
command to stop so that total elapsed time can be read. The 
circumference of a circle of 57.3 cm radius is 360 cm; there- 
fore, 1 cm on the perimeter corresponds to 1°. The protractor 
can easily be made as a homemade instrument. 


♦ 
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Figure 91. A simple model sailplane with prostretched wire 
as an indueer of turbulence. 

VII. SUMMARY 


The measurements study the behavior of wing sections in the low- 
turbulence stream (turbulence factor = 1.06) in a measuring range from 
Re = 20,000 to 170,000 in which the profiles pass through a critical region. 
This measuring region corresponds approximately to the Reynolds number 
range of model airplane wings. 


h Behavior of Airfoil Profiles in the Critical Range of Reynolds Numbers 

a) Although for the flat friction plate the critical transition from lami- 
nar boundary flow to turbulent begins at Re. < 1 million, and at Re. < 405,000 

for the sphere, the measurements on airfoil profiles show a beginning of transi- 
tion at Re < 100,000, and for flat and cambered thin plates at the lowest angle 
of attack at Reynolds numbers even below 20, 000; on the intermediate airfoil 
section N 60 (d ■ 12 percent, f = 4 percent, r - 1.4 percent of chord t) at 
Re k = f° r the thick airfoil section G 625 (d = 20, f = 6, r = 3.4) it is 

at Re, = 105,000. 
k 

b) Although for the flat friction plate the "transition" to fully turbulent 
boundary layer flow develops over a range up to Re » 10 T , the fully turbulent 
condition on the two above mentioned airfoil profiles is. finished after a Reynolds 
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i' 1 " 1 " i I »f 50,000 above Re^ and at Re = 155,000 for the thick profile 

; measurements of two thin plates (1. flat, 2. curved) , the transi- 

*' h 1,1 " considered as finished at the lower measuring limit attainable 

"Hi l be i. istallation, Re = 20,000. 


e beginning of the transition at the eritioal Reynolds number Re 

IV 

by a sudden increase in lift coefficient c and a corresponding 

a 

i' '«•**!« * i" ag coefficient c^, because the laminar separated flow previously 

!1 "'" oal region overcomes the pressure increase at the airfoil upper 
" i i . " \ I • iason of the boundary layer becoming turbulent or by the now 

■ 1 1< < l i vc « . jy transport from the outside flow, and attaches itself; at low 

•i.i- lcf: ni k this occurs near the trailing edge. 


t | )(>).» it i. 

; :u f;u‘P ;i i;tl 
i o *' Ll 1 1 ; 1 1 
. m * I s «tN :i iii 
r lit hii' 

• m » n point 


Itical transition, that is, the attachment to the airfoil upper suf- 
; a lower Reynolds number the small is the camber of the upper 
) thinner the airfoil nose, and thus at very small angles of attack 
te with a sharp leading edge. The sharp nose of a thin profile 
lence edge in connection with the great acceleration effective in 
le of attack and circulation because of the movement of the stagna- 
>w around the nose. 
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the transition region, a small increase in Reynolds number 
rapidly increasing turbulence produces a relatively large move- 
nsition point toward the leading edge and a movement of the 
it toward the trailing edge, so that now the attainable separation 
and the pertinent c & max value rapidly become larger, for 

•ofile N 60 from 4° to 16. 5° , or c = 0. 5 to 1. 3, iin the range 

) to 113,000. A measurement with a turbulence wire on the thick 

where Re, and the c increase is displaced from 
tv a max 

o Re = 20, 000, shows that die intensity of turbulence and not the 
it or the energy content of the external flow is decisive. 
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he range of the transition, the c o line in the plot of c versus 

a max a 

ith the separation line. The end of the transition evidences itself 


a break in the c & line, which now ascends only slightly 
( corresponding to the entraining effect, which increases with Re) ; 








r ^ < r . 






2. from the break on, the separation line is below the c 

line. a max 

With increasing angle of attack, the turbulent separation point moves 
toward the leading edge on the profile; if it passes the transition point, the flow 
separates, and then the separation point moves away to the pressure minimum, 
that is, to the highest point of the upper side camber. 

f) Local laminar separation is superimposed on the incipient turbulent 
transition beginning at the airfoil trailing edge at low Reynolds numbers and low 
angle of attack. With increasing turbulence, the transition point moves toward 
the leading edge, where the accompanying local laminar separation collapses. 

The transition point reaches the leading edge more quickly, the smaller is the 
mean-line camber f / 1 and the thinner is the profile or the sharper the leading 
edge. When the transition point reaches the leading edge on a thin low-camber 
profile, it forms a transition vortex which broadens with further increasing 
angle of attack, since the flow separates in laminar fashion over the leading 
edge and then again attaches turbulently at the rear of the vortex. 

Because a laminar separation precedes the transition both at the trailing 
edge as well as now at the leading edge, the conclusion is evident that each turbu- 
lent transition on a wall is preceded by a laminar separation - even if it is of 
only molecular size. The origin of the transition vortex on the leading edge is 
evident in pressure measurement by a suction peak. With further increase in 
angle of attack, the transition vortex broadens flatly to the trailing edge, as the 
transition point moves again to the trailing edge. The separation surface of 
the initially small local separation then forms a "fluid” profile upper surface 
with enlarged camber and a pronounced suction peak, which shortly before the 
final separation oscillates , in intermittent separation and impulsive attachment 
of the flow, between the highest value of negative pressure and the pressure of 
the separated flow. The separation and attachment occur from the leading edge 
to the trailing edge. 

The development of the transition vortex is similar in character to the 
separation phenomena in cavitation. In addition, just as cavitation decreases 
with increasing total pressure in water, for example, on a diving submarine, 
so with increasing Reynolds number the transition vortex probably decreases, 
so that the transition vortex has a particular importance in the Reynolds number \ 
range considered here. 

g) Flat plate. The turbulent phenomena at the flow around the nose 
explain the supercritical behavior at a Reynolds number as low as 20,000 for 


the flat and cambered plate, because of the premature widening of the transition 





sharp leading edge acting as a turbulence edge becomes a separation edge at 
~ 5 to 7 6 angle of attack, and in spite of a steep slope in the lift line 
( dc /d 3^ « 2ir) , a c of only 0. 5 at 42,000 to 0. 7 at 420, 000 is reached. 

The flat plate, after separation with further increase in angle of attack, 

shows no decrease in the c value. The c value increases continuouslv 

a a max J 

with increasing Reynolds number because of increasing entraining effect and 

diminishing transition vortex, and this explains why, in spite of low Re , the 

K 

flat plate and symmetrical profiles are so unfavorable with respect to c at 

a max 

low Reynolds numbers, and why symmetrical profiles operate more favorably 
only at high Reynolds numbers. 

h) Cambered plate. This, too, works supercritically at Re = 20,000, 

but it reaches double the c value in comparison with the flat plate; at 

bl max 

“ 8 ° » a c a max = 1.06 occurs at Re = 42,000 and 1.3 at 420,000. The 

slope of the lift line is surprisingly steep; at Re = 42,000 the value of 
d^/da^ - 2. 77ir, and thus is even greater than for the flat plate. The superi- 
ority of the cambered plate - over profiles, too - extends to Re » 100, 000. There 
are three reasons for the superiority: 

1. Favorable interaction of tangential approach flow to the leading 
edge at a high angle of attack with nose turbulence. 

2. A maximum concave camber of the lower surface so that it 
participates as far as possible in the generation of lift. 

3. The relatively small camber of the upper surface. 

The attempt to make the lower surface as heavily cambered as possible 
and the upper surface as little as possible to avoid premature separation th us 
leads to the thin cambered plate as the most favorable profile for model airplanes 
and model airplane propellers to generate high lift and low drag at low Reynolds 
numbers. 


2. Conclusions for the Model Airplane 

a) Choice of Profile. A model airplane achieves optimum flight 
performance in its size class only when its flight state is supercritical. This 
flight state can always be realized in the practically required limits by a suitable 
profile choice following the basic principle: the lower the Reynolds number the 
thinner must the profile be and the sharper the leading edge. 
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The model airplane has a low wing loading ( up to 3 kg/m 2 ) , It flies 
slowly (v (=» 3 to » 13 m/sec) , and accordingly has a high c (0.5 to 1) . The 

9 . 

cambered plate 417 a at Re = 42,000 and = 2. 5“ has the high value e -0.7 
and the minimum of profile drag c^ ^ = 0. 027. This gives a profile maxi- 
mum drag-lift ratio of c Woo /° a = 1:26 which, converted to the aspect ratio 

A = 10, is 1:16.5. This drag-lift ratio, unusually favorable at Re = 42,000 for 
a model-airplane airfoil, surpasses any profile in this Reynolds number range 
and thus explains the surprising flight performances of the room airplane model 
covered with film. A plot of the profile ceiling factor c 3 /c 2 versus the pro- 

file lift-drag ratio c / c for the measured profiles shows for the profile 

choice, with respect to camber f and thickness d, the following guiding values: 


Re 

f/t 

d/t 

r/t 

50,000 

0.06 to 0.09 

0.03 

0.004 

100,000 

0. 06 to 0. 08 

0.06 

0.007 

200,000 

0.05 to 0.07 

0.12 

0.014 


The most favorable location of the maximum camber appears to be at 
15 and 25 percent of the chord length. 

b) The Wing Plan. Wings with a tapered plan work unfavorably because 
of the correspondingly low Reynolds number at the wing tips: for the small model 
airplane, accordingly, the rectangular plan with slightly rounded tips is the most 
favorable form. This applies to the plan of the small model airplane propeller, 
too. The gain achieved in high Reynolds number by the rec tangular plan is 
greater for small model airplanes than the reduction in induced drag achieved by 
tapering of the plan. 

c) C onstant Moment Coefficient ( fixity of center of pressure) in the 
normal flying range is a characteristic of profiles with an S-shaped mean line 
only at high Reynolds number. As separation on the profile upper surface pro- 
gresses because of Reynolds number reduction or stalling, the S-shape becomes 
ineffective. Only the flat plate and the plate with the S-curve work with a fixed 
center of pressure in the region of Reynolds numbers around 20,000. 

d) Artificial Turbulence. By a turbulence wire or twine stretched along 
the airfoil nose and making turbulent the boundary layer of the airfoil upper sur- 
face the critical transition on profiles can be reduced to about Re = 20,000. 

This makes it possible to use favorable construction possibilities even for small 
model airplanes: 
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1. Profiles with a thickness of up to 20 percent of the chord 
length. 

2. Profiles with an S-shaped mean line for flying-wing model 
airplanes. 

3. Airfoils with tapered plan. 


3. Comparison: Large Airplane and Model Airplane 

a) The boundary layer should remain laminar as long as possible on the 
large airplane, to utilize the lower boundary friction as much as possible and to 
displace turbulent separation to as high an angle of attack as possible; on the 
other hand, on a model airplane, turbulence is desired as early as possible, to 
obtain high supercritical range of Reynolds numbers. For the large airplane a 
so-called laminar profile is desired, but for the model airplane a turbulence 
profile with a sharp leading edge is sought. The airplane wing should be 
especially smooth at the nose, but some roughness can sometimes be advanta- 
geous on the wing of the model airplane. In the large airplane, the elliptically 
tapered airfoil is advantageous, but on the model ariplane and on the propeller 
of the model airplane the rectangular plan is better. On the small model air- 
plane, the sharper wing nose, acting as a turbulence edge, produces the critical 
transition earlier, giving an increase in the c value; on the large airplane, 

it would act as a separation edge, so that here with an otherwise equal profile the 

wing with the higher nose radius must achieve the highest c value. All 

a max 

measures for artificial turbulence have a favorable effect on the model airplane 
but an unfavorable one on the large airplane. In any case, these are crude con- 
trasts, which require individual treatment, because the model airplane is not 
merely an airplane reduced in size. Thus model flying and aeronautical physics 
in our secondary schools are in urgent need of direction along these ways of 
thinking. 

b) The effect of force on the profile in the transition between the super- 
critical and the subcritical region is similar to the effect found during the 
approach to the velocity of sound c, as shown by a comparison of the measure- 
ment on profile G 625 with an American measurement on a similar profile of 
the same thickness ( d = 20 percent) ; the very unfavorable subcritical polar 
curve at Re = 21,000 for G 625 coincides with the polar curve at v/c = 0.85. 

The transitions are completed from the supercritical to the subcritical state of 
Re = 155,000 to 105,000, and then at v/c = 0. 5 to 0. 6 there is a new 
"hypercritical state. " Although the fluid mechanics causes are basically dif- 
ferent, the force effects are the same. The same practical result occurs for 
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the small model airplane as for the high-speed airplane: the lower the 
Reynolds number on one hand, and the more the flying speed approaches sonic 
velocity, on the other hand, the thinner must the profile be, and the smaller 
the camber and the more pointed the leading edge. 


4. Application of Results 

a) The polar curves of the measured profiles make possible for the 
first time the quantitative determination of flight performance and the stability 
of model airplanes in the Reynolds number range from 20, 000 to 170,000, New 
fundamentals are provided for the evaluation of experiments in the free-fall 
tunnel and in the vertical wind tunnel for spinning tests. More far reaching 
conclusions than were hitherto possible car be drawn for large airplanes from 
these experiments, through interposition of the polar curves measured at 
large Reynolds numbers for the same profile. The significance of the measure- 
ments for general aeronautics is involved here, along with the expansion of the 
physical picture. 

In addition, the test results find application in the design of small 
blowers, centrifugal pumps, and fans, as well as for evaluation of tests with 
model ship's propellers. 

b) The contrasts, hitherto considered as troublesome, between the 
flow laws of the large airplane and the facts of experience in model flight are 
explained in considerable part and reasons for the causes found. This has 
built a new bridge between aeronautical physics and model-airplane building in 
the intention of the "Ludwig Prandtl Prize. " 

c) It is shown that model flight is closer to Lilienthal's prototype, 

the flight of birds, than to the large airplane, because birds and the model air- 
plane fly in a region of similar Reynolds numbers. New and interesting connec- 
tions of model flight with biology instruction, in the sense of the biophysics of 
flight, result from these interrelationships. 

d) The desired impetus has been given to aeronautical joint efforts in 
the secondary schools for going beyond the usual model airplane construction 
and incorporating the constructive creative element of individual activity, so 
that the agreement between theory and practice can be proven in flight testing 
by measurement of flight in still air. Although the individual in the airplane 
factory or in the research laboratory can see only part of a large task, here 
the entire picture from design to flight testing is in the hands of teacher and 
pupil in the best form as a promising, vocationally directing incentive to the 
study of flying or aeronautics. 
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Figure 92. Summary of the five airfoil profiles measured. 
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VIII. NUMERICAL TABLES 


A correction for the measured values c and a measured on the 

w m m 

rectangular airfoil model (aspect ratio A « 5) to take into account the finite 
stream of the wind tunnel (cf. p. 90) : 

°wk6 = 01287 ^ “kb " 738 c a • 

For conversion of the test values of and a- to the elliptical wing 
plan and to infinite span ( A = <») • 


"Profile drag" c =c - c , = c - 0. 07897 c 

w«o wm wk» wm w • w ' 0£ ' , 

"True angle of attack" at the lower surface tangent a = 

= a - 4. 918 c . 00 

m a 
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a 
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TABLE 1: Polar Chart of Wing Section G 625 
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